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Die Temperaturabhängigkeit 
des Widerstandes des Nickels am Curiepunkt 


Von Bérje Svensson 
(Mit 2 Figuren) 


Die theoretisch wichtige Frage, ob die spontane Magnetisierung 
am Curiepunkt verschwindet oder ob sie nur in dem umgebenden 
Temperaturgebiet stark abnimmt, um sich dann mit steigender Tem- 
peratur dem Werte Null asymptotisch zu nähern, läßt sich durch 
die Messung der Magnetisierung kaum beantworten, weil das fiir die 
Messungen nötige Feld auch eine wahre Magnetisierung hervorruft, 
die sich der spontanen in solcher Weise überlagert, daß diese nicht 
für sich genau ermittelt werden kann. Die Beziehungen zwischen 
der spontanen Magnetisierung o und dem Widerstand o, die u. a. 
von Cabrera und Torroja!) sowie Gerlach und Schneiderhan?) 
näher untersucht worden sind, und die nach Johansson) und 

prelius®) auf einen ferromagnetischen Zusatzwiderstand von der 

:öße 
(1) do = c(o,? — 0%) 
zurückgeführt werden können, wo o, die spontane Magnetisierung 
im absoluten Nullpunkt und c eine Konstante bedeuten, öffnen aber 

1en indirekten Weg zur Prüfung der Frage. Geht do?/dt, wie die 
übliche Theorie fordert, am Curiepunkt sprungweise zu Null, muß 
sich nämlich nach (1) auch do/dt sprungweise ändern. Die Wider- 
standstemperaturkurve soll am Curiepunkt einen Knick aufweisen. 
Gerlach?) schließt aus gemeinschaftlich mit Schneiderhau vor- 
genommenen Messungen an einer Nickelprobe, daß kein solcher 
Knickpunkt vorhanden sein sollte, und scheint geneigt zu sein, an- 
zunehmen, daß die spontane Magnetisierung erst oberhalb des 
Curiepunktes allmählich verschwindet. Die Folgerung ist aber 


1) B. Cabrera u. J. Torroja, Ann. Soc. Esp. Fis. y Chim. 11. 8. 443, 
25. 1913. 
2) W. Gerlach u. K. Schneiderhan, Ann. d. Phys. [5] 6. S. 5, 772. 1930 
3) C. H. Johansson, Ann. d. Phys. [5] 84. S. 976. 1927. 
4) G. Borelius, Ann. d. Phys. [5] 8. S. 261. 1931. 
5) W. Gerlach, Phys. Ztschr. 33. 8. 953. 1932. 
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nicht bindend. Teils liegen die Messungen nicht dicht genug, um 
einen eventuell vorhandenen Knickpunkt deutlich hervortreten zu 
lassen. Teils enthielt die Probe 0,6°/, Mn, das stark auf den 
Curiepunkt einwirkt. Dieser Mangangehalt ist in Wirklichkeit 
groß genug, um den Curiepunkt um etwa 15° herabzudrücken und 
schon ganz kleine Konzentrationsschwankungen müssen genügen, 
um den Knickpunkt unscharf zu machen. Es ist deshalb hier im 
Zusammenhang mit einer Untersuchung über den Widerstand der 
Ni-Cu-Legierungen, die später veröffentlicht werden wird, auch eine 
neue Untersuchung über den Widerstand an reinstem Nickel in der 
Umgebung des Curiepunktes vorgenommen worden. 

Die Messungen sind an einer von der Firma Hilger in London 
bezogenen Nickelstange von 12cm Länge und 5 mm Dicke aus- 
geführt worden. Nach dem beigegebenen Analysenschein enthält 
dieses Nickel 0,102°/, Fe und insgesamt 0,036 Mg, Ca, Cu und Si. 
Für Stromzuführung und Spannungsabnahme dienten Nickeldrähte, 
die in Ausbohrungen an den Stabenden befestigt waren. Der Nickel- 
stab befand sich bei der Messung in einem evakuierten Quarzrohr, 
das durch einen mit Wechselstrom betriebenen Rohrofen von 45 cm 
Länge geheizt wurde. Um konstante Temperatur in der Umgebung 
der Probe zu erzielen, war das Quarzrohr noch von einem dick- 
wandigen Kupferrohr umgeben. Die vom Heizstrom erzeugte Feld- 
stärke an dem Ort der Probe war höchstens 2,5 Gauss. Die Tem- 
peraturgeschwindigkeit betrug etwa 0,4 °/Min. 

Die Widerstandsmessung geschah mit einer Thomsonschen 
Doppelbrücke und der Strom durch den Nickelstab betrug etwa 
60 mA. Zur Temperaturmessung diente ein Pt-PtRh-Element, dessen 
Schweißstelle elektrisch isoliert an der Mitte des Nickelstabs be- 
festigt war. Die elektromotorische Kraft des Thermoelements wurde 
mit einem Kompensationsapparat von Wolf gemessen. Das Thermo- 
element wurde an der Siedetemperatur des Wassers und den Schmelz- 
punkten von reinem Zinn, Cadmium und Zink kalibriert und die 
zwischenliegenden Temperaturen mit der Hilfe von Tabellen in den 
„Critical Tables“ interpoliert. 

: Die Resultate sind in der Fig. 1 zusammengestellt worden. Als 
_ Abszisse dient die Celsiustemperatur, als Ordinate die Größe 
0, R,/R, (wo o, den spezifischen Widerstand bei 0°C und R, und 
R, die absoluten Widerstände bei £ und 0° C bezeichnen) in Mikroohm 
x Zentimetern angegeben. Der Deutlichkeit halber sind die ver- 
schiedenen Kurven gegeneinander um 1 u2/cm parallel verschoben. 
Es sind Messungen sowohl bei steigender als auch bei abnehmender 
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In diesen Kurven tritt ein Knickpunkt am Curiepunkte recht 
deutlich zum Vorschein. Seine Lage kann aus den Kurven zu 
’ 2 „ 363,0 + 0,2° C bestimmt 
3 a 3° werden. Das im Nickel be- 
ae findliche Eisen diirfte nach 
Peschard') eine Erhöhung 
des Curiepunkts um 0,1° be- 
wirken, der Curiepunkt des 
Nickels ergibt sich also aus den 
Messungen zu 362,9 + 0,2°C. 
Die kleinen gegenseitigen Ver- 
schiebungen der Kurven 2—4 
in der Ordinatenrichtung sind 
wahrscheinlich Wirkungen 
der Entglühung und Ent- 
gasung zuzuschreiben. Für 
die Kurve 1 liegen nur rela- 
tive Messungen vor, die Ab- 
solutwerte sind den übrigen 
angepaßt. 
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Fig. 1. Temperaturwiderstandkurven _ Fig. 2. Temperaturwiderstand- 
fiir reines Nickel ae kurve fiir reines Nickel 
in der Nähe des Curiepunkts a zwischen 0—500° C 


1) M. Peschard, Théses Strassbourg 1925. 
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Fig. 2 zeigt die Widerstandskurve über einem mehr ausgedehnten 
Temperaturgebiete. Auch in diesem reinen Nickel ist, wie bei der 
Probe von Gerlach und Schneiderhan, die Kurve über dem 
Curiepunkt noch gebogen und erreicht erst einige Zehnergrad höher 
einen geradlinigen Verlauf. Es dürfte doch kaum richtig sein, aus 
diesem Beiunde auf das Vorhandensein einer spontanen Magneti- 
sierung auch oberhalb des Curiepunktes zu schließen. Möglicher- 
weise steht die Sache mit dem in diesem Gebiete vorhandenen 
starken Paramagnetismus in Beziehung. Es ist in diesem Zusammen- 
hange zu bemerken, daß auch die Metalle Pd und Pt, die einen 
starken temperaturabhängigen Paramagnetismus aufweisen, nach 
unten gebogene Widerstandskurven haben. 


Die vorliegenden Resultate sprechen sehr für die gewöhnliche 
Auffassung, daß die spontane Magnetisierung am Curiepunkt aufhört. 

Es sei mir gestattet, Herrn Prof. G. Borelius für viele er- 
giebige Besprechungen meinen herzlichsten Dank auszusprechen. 


BE Stockholm, Physikalisches Institut der Technischen Hochschule. 


(Eingegangen 9. Dezember 1934) 
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L. Posener. Theorie der freien elastischen Schwingungen usw. 


Ein Beitrag 
zur Theorie der freien elastischen Schwingungen 
von Zylindern und Rohren 


Von Ludwig Posener 


Inhalt: Einleitung. — 1. Dehnungsschwingungen einer dünnen zylin- 


drischen Schale: a) Ableitung der Bewegungsgleichungen; b) Widerlegung 


eines Einwandes; c) Axialsymmetrische Schwingungen; d) Grenzfille. — — 


2. Bemerkungen zur Theorie der symmetrischen Schwingungen eines Voll- 
zylinders: a) Die Pochhammersche Lösung; b) Beweis für die Nicht- 


existenz einer vollständig separierten Lösung, die alle Randbedingungen 
erfüllt; ce) Entwicklung der Lösung in Fouriersche Reihen nach 2; d) Er 


gänzung der periodischen Lösung. — 3. Zur kopplungstheoretischen Inter- 
pretation der elastischen Schwingungen: a) Kraftkopplung und Beschleunigungs- 


kopplung; b) Die Kopplungstheorie und die Schwingungen des Vollzylinders. 


Einleitung 


In einer Reihe von Arbeiten, die in den Jahren 1931-1983 


in den Annalen der Physik veröffentlicht worden sind, haben 


E. Giebe, A. Scheibe und E. Blechschmidt’) die freien Schwin- | 
gungen elastischer Stäbe und Rohre untersucht. Die Schwingungen 


wurden bei kristallinischem Material (Quarz) piezoelektrisch, bei 
Nickel- und Eisenstäben magnetostriktiv angeregt, und die Serien 
der Eigenfrequenzen bis hinauf zu hohen Ordnungszahlen aus- 
gemessen. Es ergab sich eine mit steigender Ordnungszahl zu- 


nehmende Abweichung vom Charakter der einfachen harmonischen 


Schwingungen, wie die Rayleighsche Theorie des „dünnen“ Stabes 


sie fordert. Auch die durch die ,,Dickenkorrektion“ verbesserte 


Rayleighsche Theorie vermag die Versuchsergebnisse nicht wesent- 
lich besser darzustellen. 


Giebe, Scheibe und Blechschmidt versuchten, ihre Er- _ 


gebnisse in folgender Weise zu deuten: Sie nahmen an, daß die 


gesamte Energie des schwingenden Gebildes sich auf die zwei 


bzw. drei Freiheitsgrade in ähnlicher Weise verteilt wie bei den 
Schwingungen „gekoppelter Systeme“. Als schwingende Grund- 
systeme wurden die Grenzgebilde betrachtet, die bei Vernach- 
lässigung einer bzw. zweier Dimensionen entstehen. So wurde z. B. 


1) E. Giebe u. A. Scheibe, Ann. d. Phys.[5]9. S. 93. 1931; E. Giebe 
u. E za Ann. 4. Phys. (5) 18. 8. =" u . 457. 1933. 
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die Schwingung eines zylindrischen Stabes zusammengesetzt aus den 
Schwingungen eines unendlich dünnen („Rayleighschen“) Stabes 
und denen einer unendlich dünnen Kreisscheibe. Die „Kopplung“ 
kann entweder in der kinetischen oder in der potentiellen Energie 
oder auch in beiden angenommen werden (Beschleunigungs-, Kraft- 
oder gemischte Kopplung). Aus experimentellen Gründen ent- 
schlossen sich die Verf. zur Annahme einer reinen Beschleuni- 
gungskopplung, was auch der modellmäßigen Vorstellung am besten 
entspricht. Die Theorie gestattete auch, über die auftretenden 
Kopplungskoeffizienten Aussagen zu machen. Sie werden im wesent- 
lichen von den Elastizitätskonstanten des Materials geliefert. 

Mit diesem einfachen Modell gelangt man theoretisch zu einer 
Doppelserie von Frequenzen, die in nullter Näherung (unendlich 
dünner Stab) natürlich mit den Rayleighschen übereinstimmen 
und in erster Näherung die Dickenkorrektion liefern. Im Rahmen 
der Meßgenauigkeit (etwa 0,1 °/,) stimmen sie ausgezeichnet mit 
den gemessenen überein, besonders gut bei dünnen zylindrischen 
Rohren aus isotropem Material (Eisen oder Nickel), für die Giebe 
aus seinem Modell die schon von Love’) für symmetrische Dehnungs- 
schwingungen angegebene Formel herleiten konnte. 

Die hier skizzierte modellmäßige Vorstellung bedarf einer Be- 
gründung aus den Grundgleichungen der Elastizitätstheorie. Diese 
gelingt ohne Schwierigkeiten im Fall des Hohlzylinders und soll 
im folgenden dargestellt werden. Dagegen gelingt sie nicht im Fall 
des Vollzylinders, dessen Theorie erheblich komplizierter ist. Be- 
kanntlich ist es bisher nicht gelungen, die strenge Lösung der 
Grundgleichungen zu finden, die zugleich die Bedingungen der 
Kräftefreiheit an der Oberfläche befriedigt. Außerdem scheint in 
diesem Falle die Giebesche Theorie nicht mit der allgemeinen 
Elastizitätstheorie vereinbar zu sein. Diese Fragen sollen im 
zweiten Teil der Arbeit untersucht werden. Im ersten Teil wird 
eine vereinfachte Theorie der Dehnungsschwingungen des Hohl- 
zylinders hergeleitet, die im Gegensatz zu der Darstellung bei Love 
nicht von flächentheoretischen Hilfsmitteln Gebrauch macht. Der 
dritte Teil ist einer genaueren Analyse der kopplungstheoretischen 
Zusammenhänge gewidmet. 


. Dehnungsschwingungen einer dünnen zylindrischen Schale 


a) Ableitung der Bewegungsgleichungen 
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Dehnungsschwingungen versteht man solche Schwingungen, deren ‘ 
potentielle Energie ausschlieBlich von den Verzerrungen der Mittele 
fläche der Schale herrührt, während bei den Biegungsschwingungen 

die Mittelfläche unverzerrt bleibt, d.h. kein Linienelement dieser 
Fläche bei der Bewegung seine Länge ändert. Die potentielle 
Energie dieser Schwingungen rührt also von den Verzerrungen der ER 
übrigen Schichten der Schale her, die linear mit der Entfernung eh 
von der Mittelfläche wachsen. Ist die Dicke 2h der Schale klein, — ee 
so wird die potentielle Energie der Biegungsschwingungen in erster 5 
Näherung quadratisch in h, während die kinetische Energie linear 
in h wird. Daher müssen die der Biegungsschwin- 


sind von derselben Galion in h, und die Frequenzen hängen = 
nicht von der Dicke ab’). 


Dehnungsschwingungen mathematisch formulieren. Wir rechnen in ¢ = 2 


Zylinderkoordinaten r, p, 2 und bezeichnen die Verschiebungen in 
radialer, azimutaler oni axialer Richtung bzw. mit u, v und w, pq a 
Komponenten des Deformationstensors mit ¢,,, € 
und die des Spannungstensors mit T,,, T,., T.,, T,,; 


az 
Dann gelten die Beziehungen ?): 


Or’ 

du dw 
Oz TR Or ’ 


rg? Cos ®: 


den 
ibes 
ing“ 
rgie 
aft- 
nt- 
« 
ten 
den 
nt- 
ner 
4 N 
1en 
mit 
1en 
be 
1 dw . öv [ev 
e e = — —+—, = = 
ler und 177 
rd | + r + r „rt Oz + ou r + 
, (du u 1 dv Ow dw 
T,. = — 3=+75-] + 20: 
v zz or r r Oz Oz as, 
2 
er 1) Vgl. auch Love „Dehnungsiose Deformationen einer dünnen Schale“, 
3 2) Love-Timpe, S. 66. 
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Man gewinnt sie durch Transformation aus den Ausdrücken für 
kartesische Koordinaten. 

Der Ausdruck für die elastische Energie der Schale lautet in 

Zylinderkoordinaten: 

ath 22 | 


Ein Linienelement der Mittelfläche =a verändert seine Länge bei 
der Bewegung um?): 


2 2 Aad, 
+ + 2e,,m nN, 
wo m und n die Richtungskosinus gegen die Parameterlinien 
gy =const und z= const bedeuten. Sollen reine Biegungsschwin- 
gungen vorliegen, d.h. sämtliche Längenänderungen verschwinden, 
so müssen die Deformationsgrößen e,, €,,, €,, einzeln verschwinden. 
Daher ist die Energiedichte der Biegungsschwingungen: 


Es liegt die Vermutung nahe, daß der zweite Summand des Inte- 
granden in (3) die potentielle Energiedichte der Dehnungsschwin- 
gungen darstellt. Diese Vermutung wird durch folgende Überlegung 
gestützt: Da die Frequenzen der Dehnungsschwingungen von der 
Dicke der Schale unabhängig sind, kann man diese beliebig dünn 
machen. An beiden Mänteln gilt nun die Bedingung der Spannungs- 
freiheit, die das Verschwinden derjenigen Komponenten des Span- 
nungstensors fordert, die einen Index r enthalten. Fordert man nun 
als „Dünnheitsbedingung“, daß diese Randbedingung noch im ganzen 
Innern der Schale erfüllt ist, so muß identisch gelten: 


(D) = = 9; 


e = 


und die Energiedichte der Dehnungsschwingungen wird: 
fo=} Coy +e, +¢,, 
Im folgenden werden wir uns ausschließlich mit Dehnungs- 
schwingungen beschäftigen. 
Aus der Bedingung T 
eliminieren 


(4) 


„=0 und Gl.(2) kann man Öu/ör 


4 


1) Love-Timpe, S. 64. 
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so daß wir schließlich erhalten 


1[ [u 1 dv , dw\2 


In diesem Ausdruck kann man die Laméschen Elastizititskonstanten A _ 
und u ersetzen durch den Youngschen Modul E und die Poissonsche © 
Konstante o, die eine mehr anschauliche Bedeutung haben. 
kinetische und potentielle Energie des Systems werden dann: 
ath I 2x 4 


u 1 dv dw\? 
a—h 
ath 1 22 


a-h 0 0 


r 


tet 


Nun muß berücksichtigt werden, daß die Schale sehr dünn sein 
soll. Beide Energien sind von der Ordnung h. Um die „Flächen- 
energien“ zu pi muß man bilden: 


lim - Vh)=V, und Th Tih) = T, - 
h—>0 
Nach dem Mittelwertsatz der ae kann man die Inte- 


gration nach r ausführen und erhält als Flächenenergien: 


wists 


“Oh 


8) 
dz\a ' a Og) 4\a Og Oz 
und entsprechend in T,. 


E 1 - 1 (4 dw öv\2 


jee 


Jetzt sind die Verschiebungen u, v, w nicht mehr als Funktionen "1 


von r zu betrachten. 


Mit diesen Energien bildet man das Hamiltonsche Integral ' xe 


des Problems, das zum Minimum zu machen ist. 


el + (ar) + Ge) | 


r 
in 
It. 
el 
n- rdrdgdz, 
n, | 
+ 
j 
e- 
ar 
n 
| 
| * 
r 
(9) | (“+ 1 dv . dw 
+5 _ 1/1 dw dv)? 
+0 |dz \a a dp 4 \a 59 + |} 
-Mn 
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Die Eulerschen Gleichungen für dieses Variationsproblem lauten: 


ov + 


1-0? öz 

é ou 

l+o | dz 2 a Og\dz 
dw 


E — 
2(l+n |a dgdz 


Dieses Verfahren könnte dem Bedenken begegnen, daß hier ohne 
Berücksichtigung der Nebenbedingungen T,, = T,, =T,,= 0 variiert 
worden ist, denen u, v, w unterliegen. Es ist aber klar, daß u, v, w 
nach den vorgenommenen Eliminationen und nach der Tntegrakien 
nach r andere Funktionen geworden sind als die, die in den Neben- 
bedingungen auftreten. Der Übergang zur „Flächenenergie“ bedeutet 
auch den Übergang zu einem neuen Variationsproblem, das nun un- 
abhängig von Nebenbedingungen behandelt werden muß. 

Man kann unser Verfahren noch auf eine andere Weise recht- 
fertigen, die in mancher Hinsicht befriedigender ist, aber die Zu- 
sammenhänge nicht so deutlich hervortreten läßt. Legt man näm- 
lich den allgemeinen Ausdruck (3) für die elastische Energie zugrunde 
und betrachtet einen Hohlzylinder endlicher Dicke, so erhält man 
als Hamilton-Integral 


-rdrdgdzat. 

Die Eulerschen Gleichungen dieses Integrals sind die be- 
kannten Bewegungsgleichungen in Zylinderkoordinaten: 
OT, 
ot 
w 


(11) 0 


| 
(I) Q i-# +o 4a 
‘ oy FL 4 
(I) | dw 
dell’ 
| 1 69 
| 0 4 
= 
| 
"3 
R 
— 


en: 
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Als Randbedingungen muß man das Verschwinden sämtlicher 
Spannungen an beiden Mänteln und an den Endflächen verlangen, also: 


(14) T,,=0, T,,=0, T,,=0 für r=(a—h), (a+h). 


(15) T,,=0, T,,=0, T,=0 fir z=0l. 
Läßt man jetzt die Dicke 2h immer kleiner werden und ver- a 
langt, daß mit abnehmendem h auch T,,, T,,, T,, gegen Null a 
gehen, so vereinfachen sich die Gleichungen zu = 
° Fr do dz 
Die Verschiebungen und Spannungen, die in diesen Gleichungen fa 
auftreten, sind nur noch als Funktionen von ¢ und z zu betrachten. deg 


Ausdrücke wie 0u/Or oder Öv/Ör, die noch formal in den Span- 
nungsgréBen auftreten, sind durch die Nebenbedingungen (14) zu oie 
eliminieren. Dann erhielten wir aus (11’), (12’), (13’) die Bewegungs- “4 
gleichungen (I), (ID), (IID) von S. 106. er 


ce) Axialsymmetrische Schwingungen 
Zu diesen Gleichungen kehren wir nun zurück. Betrachtet ig 
man nur axialsymmetrische Schwingungen, setzt also alle Ab- bs. 
leitungen nach @ gleich Null, so lassen sich Dehnungsschwingungen Be 
und Rotationsschwingungen trennen; denn die Gl. (IT) und (IT) ent- ig 
halten © nicht mehr, und für das letztere liefert Gl. (III): 2 

2(l+o) dz? 


also die Frequenzen: 


= = u 


in Ubereinstimmung mit Love-Timpe S. 624. 
nehmen die Gestalt an: 


v= A-sinpzett, 


Eo 


E 


Mit der Rar 


4 
— 
4 
q a 
| 
4 
ne 
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on 
n- - 
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‘Die eT der Gl. (I) und (II) können in der Form 
u= A-sinpzeirt, 
w= — Ücospzeir! 
angesetzt werden, wo p durch die Bedingung (16) bestimmt wird. 
Es ergibt sich: 
ne sinpz | = =0 oder p= > 
z= 0,1 


“- 4 undC sind durch die folgenden linearen Gleichungen verbunden: 


a el-0) a 
Eo P E 2 
Die Frequenz v bestimmt sich durch die Determinantengleichung: 
E 1 | 
o(l — 
~ — 0%) 


— 


d) Grenzfälle 


Geht man zur Grenze } 0 über, so erhält man die radialen 
Man 


u und die zweite Spalte durch p nn und geht dann zur Grenze 


— © über. 


1 
— . 
a @ 
; Die Longitudinalschwingungen eines sehr diinnen Stabes stellen 
_ den zweiten Grenzfall dar. Multipliziert man GI. (I) mit a und 
läßt a gegen Null BER 


3 
5 
of 
er 
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= 
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Dies ist die bekannte Formel für die Kopplung der longi- 
tudinalen Schwingungen eines sehr dünnen Stabes mit den trans- 
versalen. Die zugehörigen Frequenzen bestimmt man entweder aus 
der 

E 


oder man geht mit (21) in die Gl. (II) ein und erhält die bekannte 
Rayleighsche Differentialgleichung: 
d? w 0? w 
Um die zwischen den Komponenten u und w der symmetri- 
schen Schwingung bestehenden Beziehungen besser übersehen zu — 


können, gehen wir mit dem Ansatz 

u = U(f)sin pz, Re 

in die Gl. (I) und (II) ein: 


@®U E 1 
(I) dt? + (1 — a 


ew 
(I) + p? W it) — 


Dies sind „Kopplungsgleichungen“ im Wienschen Sinne. col = 
besteht ng mit den Kopplungskoeffizienten 
pe 


Soll die Kopplung ae werden, muß o = - 0 gesetzt 
werden. Nun steht aber der Kopplungskoeffizient o auch in den 
Koeffizienten von U in Gl. (I) und W in Gl. (II), die eigentlich m 4 
Frequenzen der Teilsysteme darstellen. Es ist nicht eindeutig 
entscheidbar, ob auch in diesen Koeffizienten o = 
oder nicht, d.h. ob 


die Teilfrequenzen sind. 
Eine genauere Untersuchung der Kopplungstheorie wird lehren, Sun 
wie diese Frage zu entscheiden ist. ‚Giebe RN daß ww 


re 
4 


d 
d. 
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Kopplung reine Beschleunigungskopplung sei und daß eindeutig der 
Fall a) vorliege. Nun kann man aus den Differentialgleichungen (I) 
~ und (II) auf 8.107 ganz einfach zwei äquivalente gewinnen, die 
durch „Beschleunigungskopplung“ miteinander zusammenhängen. 
Man setze: 


aw 

dt? 

Setzt man in diesen Gleichungen o = 0, so erhält man die 
Differentialgleichungen der Radialschwingungen des Kreisringes und 
der Longitudinalschwingungen des dünnen Stabes mit den Fre- 


U 


Bow a? 
+ W +oap 


E 


(22) 


man dasselbe, wenn man in Gl. (I) / oder in Gl. (II) a —> 0 gehen läßt. 

Damit ist Giebes Behauptung bewiesen: Die symmetrischen 
Schwingungen eines dünnen Rohres können aufgefaßt werden als 
Schwingungen eines gekoppelten Systems, dessen Grundgebilde ein 
schwingender Stab und ein symmetrisch schwingender Kreisring 
sind und bei dem Beschleunigungskopplung mit dem Kopplungs- 
faktor o herrscht. 

Die Interpretation der unsymmetrischen Schwingungen auf 
Grund der Giebeschen Theorie scheint unmöglich zu sein. Denn 
die Integration der Gl. (I), (II) und (III) würde auf eine Säkular- 
gleichung dritter Ordnung für »? führen. Nach Giebe stehen aber 
als Ausgangssysteme nur der unendlich dünne Stab und der Kreisring 
zur Verfügung, und man sieht nicht, wie aus der Kopplung dieser 
beiden Systeme eine Gleichung dritter Ordnung für »? gewonnen 
werden kann. Selbst wenn man die (unsymmetrischen) Oberschwin- 
gungen des Kreisringes hinzunimmt, erhält man zwar mehr Fre- 
_ quenzen, aber immer nur Gleichungen zweiten Grades zu ihrer Be- 
stimmung. 


1) Vgl. E. Giebe u E. Blechschmidt, 
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2. Bemerkungen zur Theorie der symmetrischen Schwingungen 
eines Vollzylinders 


a) Die Pochhammersche Lösung‘) 
Das Problem der symmetrischen Schwingungen von Voll- 
zylindern ist zum erstenmal von L. Pochhammer im Jahre 1876 
in Angriff genommen worden. Die Differentialgleichungen des — 
Problems lassen sich unmittelbar aus den allgemeinen Gleichungen _ 
für Zylinderkoordinaten entnehmen [GL (11), (12), (13), 8.106), 
wenn man darin alle Ableitungen nach gleich Null setzt. Sie 


u ou 


te 


Ou +52 1 
Ov v 

Da v nur in der Gl. (C) auftritt, lassen sich auch hier Torsions- __ 
schwingungen und Dehnungsschwingungen trennen. Fiir die ersteren me je 
findet man die vollständige Lösung bei Love S. 334. Das mathe- 
matische Problem der Dehnungsschwingungen ist schwieriger, weil 
es sich um das System der beiden partiellen Difleventialgleichungen (A) 
und (B) handelt mit den Randbedingungen 


; identisch in 
(24) = 0 z und t. 


für z= 0,1 


identisch in r, t. 


Pochhammer integriert die Gleichungen durch lineare Kombination Sk 
einer divergenzfreien und einer rotationsfreien Lösung, und es 
lingt ihm, drei der Randbedingungen zu erfüllen. 

lautet: 


(27) ws [Ahd, (hr) + (pd, (kr)] sin pzeirt. 
(28) w =[— ApJ, (hr) + BkJ,(kr)] cos pzert, 
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bedeuten. Die Bedingungen am Mantel ergeben zwei homogene 
lineare Gleichungen für (A) und (B), die sich nur befriedigen lassen, 
wenn die Determinante 


2h aJ, (kR) 
30) | — 
— 2hp J, (hR), (k?— p) J, (kR) 


verschwindet. Das Verschwinden der Normalspannung an den 


Enden [Bedingung (26)] erfordert p = u. Damit sind alle Kon- 


stanten bestimmt, und es gelingt nicht, gleichzeitig die Schub- 
spannung an den Enden zum Verschwinden zu bringen [Bedingung (25)). 
Eigentümlicherweise ist die Frequenz aus den Gl. (29), (30) bereits 
bestimmt, ohne daß alle Randbedingungen erfüllt sind. Für sehr 
dünne, lange Stäbe jedoch geht die Pochhammersche Lösung in 
die richtige über, denn da die Schubspannung am ganzen Mantel 
verschwindet, kann sie an den Endflächen keine beträchtlichen 
Werte mehr annehmen. Auch analytisch sieht man das ein, wenn 
man die in den Lösungen (27), (28) auftretenden Besselschen Funk- 
tionen nach Potenzen von r entwickelt und bei aor ersten (bzw. 
zweiten) Ordnung abbricht. Es wird dann 

oz’ 
und iie Frequenzengleichung geht in die für den unendlich dünnen 
Stab über. 


b) Beweis fiir die Nichtexistenz einer vollstindig 
separierten Lösung, die alle Randbedingungen erfüllt 


Will man untersuchen, warum das Pochhammersche Ver- 
fahren nicht zum Ziele führt, so muß man zunächst den Ansatz be- 
trachten: Alles weitere ist dann zwangsläufig. Pochhammer setzt 
die Funktionen u und w in der Form an: i rie 


(31) u= Alnett-r9, 


Dieser Ansatz erscheint willkürlich. Unter allen Umständen muß 
indessen die Lösung den Zeitfaktor periodisch enthalten; aber in 2 
braucht sie a priori nicht periodisch zu sein. Es läßt sich aber 
zeigen, daß die Pochhammersche Lösung die einzige ist, die zu- 
gleich in allen drei Variablen separiert, die Zeit periodisch enthält 
und drei der vier Randbedingungen befriedigt. Dies beweisen wir so: 
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lassen sich die Differentialgleichungen (A) und (B) folgendermaßen 
schreiben, wenn man die Dichte mit 1 setzt: 


ö 
(A 


Wir differenzieren Gl. (A) z, Gl. (B) nach r und 
dem wir 


ou Ow 
(82) @z u Ww, 
als neue Variable eingeführt haben: 
ey 


Subtraktion dieser beiden 
a 
33) —p(Poy+ —p (Pow 4+ 
Durch Anwendung des ee 
o? 
auf Gl. (A’) kann w’ mit Berücksichtigung der Relation (33) elimi. ert 
man erhält eine Gleichung vierter für 


0? 
Diese Gleichung vierter Ordnung kann man ee schreiben: 
a? a 
(84a) — @ + 2m) (Po + - + )|y=®. 
Es wird also auf die unbekannte Funktion y das Produkt zweier 


vertauschbarer Differentialoperatoren zweiter Ordnung angewandt. 
Die Lösungen = beiden Gleichungen a Ordnung 


(36) az? 


sind sicher auch Lösungen von (34). Außerdem könnte noch eine 
Lösung auftreten, oo. dem System 


genügt (oder umgekehrt). 
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a) Untersuchen wir zuerst die Gl. (35) und (36). Verlangt man, 
daß die Lösung separiert ist, d.h. sich in der Form ansetzen läßt: 
(37) 
ergibt Gl. (35): 

vfir)- pat 2wfir 
und nach durch f(r)-p(2(a 
A+2u f(r) (2) 

Die Summe aus einer Funktion von r und einer Funktion von z 
kann aber nur konstant sein, wenn die Summanden konstant sind: 
Pofr) ___ pe. at 

Integration dieser beiden gewöhnlichen Differentialgleichungen 
liefert: 

fi = A-d, (hr); (2) + sin (yz + 9) 
mit 
h?+ 72 = 

Analog ergibt (36): 

y, (r) = B- J, (kr) sin (#2 + ©) 
mit 

ß) Der Fall, daß y dem System (Q) genügt, liefert nichts Neues. 
Denn aus der ersten Gleichung des Systems folgt: 
yY=R-.J,(hr). sin (yz + 0). 
Setzt man dies in die zweite ein, so ergibt sich die inhomogene 
Differentialgleichung: 
(38) + uf (2) = Red, (hr) sin yz. 
Division durch f(r) - p(2)- u liefert: 
+ _ J, (hr) sin (yz +9) 


Diese Gleichung ist von der Struktur: 
const + x(r) + 7(z) =a(r)- 
Differentiation nach r und z gibt: 


(7) = 
a’ (r) = 0 


Also entweder 


oder 
(2) 0, 
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wir eine gewöhnliche inhomogene Differentialgleichung für f, deren 
vollständiges Integral (unter Fortlassung der durch die Natur des 
Problems ausgeschlossenen Besselfunktionen zweiter Art) lautet: 


f(r) ~ J, (hr) +C-J, (kr). 


Und dasselbe ergibt sich, wenn wir die andere Alternative (r) = 0) 
wählen. 
Wir haben also: 


(40) w=[E-J, (hr) sin (yz + 6) + Jd, (kr) sin (Bt + eirt. 

Durch Integration nach z bzw. r ergeben sich w und w: 

39) u=[C-J, (hr) cos(yz + 0) + D-d, (kr) cos(@z + 8) + «(rj]eir' | 

40) w=[G-J, (hr) sin(yz + 0) + H-J,(kr)sin(d2 + 8) + B(aje”’. 


Die Funktionen «(r) und #(z) werden durch die ursprünglichen x A 
Differentialgleichungen (A) und (B) bestimmt zu 


ae 


a(r)=M-J,(§r), = 


Die Randbedingungen am Mantel müssen identisch in z gelten. 
Dies läßt sich nur erreichen, wenn M=N=0, B=y, und das 
führt auf die Pochhammerschen Lösungen. 

Damit ist bewiesen, daß die Pochhammerschen Lösungen die | 
einzigen vollständig separierten Lösungen der Differentialgleichungen _ 
sind, die drei Randbedingungen erfüllen und daß es somit keine 
vollständig separierte Lösung gibt, die alle vier Randbedingungen 
erfüllt. 


by 


7 


ec) Entwicklung der Lösung in Fouriersche Reihen nach z. 


Die Lösung muß jedenfalls periodisch in ¢ sein, also von der 
Form 
u = Alr,.).eirt 
w= Bir,2)-eirt, 


Die Funktionen A (r, z) und B(r,z) dürfen in den Variabeln r und z © 
nicht separiert sein. Der nächste Weg, der sich zu einem allgemei- — 
neren Ansatz für A und B bietet, ist die Entwicklung nach einem 
vollständigen Orthogonalsystem in z. Da die Differentialgleichungen ai 
homogen von der zweiten Ordnung in z sind, bietet sich von selbst die = 
Entwicklung in Fourierreihen dar. Die ersten und zweiten Ableitungen = 
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der Funktionen nach z lassen sich dann sicher auch als Fourierreihen 
darstellen; es ist aber fraglich, ob gliedweises Differenzieren erlaubt 
ist, da man wohl die Funktionen über ihren Definitionsbereich (Stab- 
länge) hinaus periodisch fortsetzen kann, aber nicht sicher ist, daß 
diese Fortsetzung stetig und mit stetigen Ableitungen ausgeführt 
werden kann. Die Bedingung für gliedweises Differenzieren ist die 
gleichmäßige Konvergenz der Reihe selbst und die Konvergenz der 
Reihe der Ableitungen. A priori kann man nicht sagen, ob die 
gesuchten Funktionen diese Bedingung erfüllen. 

Die physikalische Bedeutung der stetigen, differenzierbaren, 
periodischen Fortsetzbarkeit ist diese: Setzt man an den Zylinder 
einen zweiten gleicher Länge an, so bleiben beide während der 
Bewegung immer in Berührung, und der Charakter der Schwingung 
wird nicht gestört. Oder, anders ausgedrückt, die n-te Oberschwin- 
gung des Zylinders hat genau dieselben Eigenschaften wie die Grund- 
schwingung eines Zylinders von der Länge !/n. Diese physikalische 
Vorstellung liegt sowohl der Theorie von Pochhammer als auch 
der von Giebe zugrunde. Daß die mathematischen Voraussetzungen 
dafür nicht vorliegen, soll im folgenden gezeigt werden. 

Der Einfachheit wegen machen wir die Länge des Zylinders 
gleich 2a (Anfangspunkt — x, Endpunkt + 2). Dann lassen sich u 


und w in der Form ansetzen: 
| 4 


w=»>ß,(r) sinnz-eirt, 

Der reelle Ansatz ist einfacher als der komplexe und bedeutet keine 
Einschränkung der Allgemeinheit, zumal für u aus den Gleichungen 
das Verschwinden der sin-Glieder folgt, wenn w antisymmetrisch in 
bezug auf z = 0 angesetzt wird. Die Differentialgleichungen werden 
dann: 

Die — a, O+A+2u)P ta, (nj cosnz=0. 
2p)n*) PB, (r)—A+p)n Fa, sin nz =0, 
Sie haben auch die Gestalt von Fourierreihen; also muß jeder 
Fourierkoeffizient einzeln verschwinden, d. h. es muß gelten: 


)+ (A+ 2p) Pda, 
—(A+2u)n*)B = 0. 


Dies sind gewöhnliche für «,(r) und /,(r) 
Ihre Lösungen sind dieselben Besselschen Funktionen, die auch 
in (27), (28) auftreten. Auch diese Methode führt also nicht über 
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die Pochhammerschen Lösungen hinaus. Daraus ergibt sich, d aS Er. 
u und w nicht als periodische Funktionen in 2 angesetzt werden ee = i 
können. 


d) Ergänzung der periodischen Lösung. 

Nun kann man aber jede Funktion, von der Existenz und Stetig- = “a 

keit bis zur n-ten Ableitung gefordert wird, durch Hinzufügen eines 
Polynoms von höchstens (n + 1)tem Grade stetig und periodisch 
über ein gegebenes Grundintervall hinaus fortsetzen. Man erkennt 
dies leicht so: Sei y die gegebene Funktion. Dann entsteht durch 
Hinzufügen einer linearen Funktion die neue Funktion z(2), de 
stetig sei: 


Die Periodizitätsforderung lautet: 


Daraus bestimmt sich b. Das Verfahren kann nach ein- oder mehr- a 
fachen Differentiationen fortgesetzt werden. Da wir Ableitungen bis — by : 
zur zweiten Ordnung brauchen, müssen wir zu u und w Polynome © 

dritten Grades in z hinzufügen. Die so ergänzten Funktionen können 
dann in Fourierreihen entwickelt und in die Differentialgleichungen 
eingesetzt werden. Vorher können wir aber noch einige Einschrän- 
kungen aus Symmetriebetrachtungen und aus den Randbedingungen 
für z =+ n herleiten. Wes 


(53) (r)cosnz+ A(r)+ Bir)z + + D(r)2?]. 
(54) w (r) sin nz + a(r) + + e(r)2? + dir)z*). 


Aus der Konfiguration des Systems geht jedenfalls hervor, dB u _ 
symmetrisch in z sein wird. Daraus folgt, daß sich das Poly- 
nom A+Bz+Cz?+Dz? auf das quadratische Glied reduziert. _ 
Für w wollen wir nichts voraussetzen. Die Bedingungen an den End- _ 
flächen müssen für z=— 2a und z2=+n erfüllt sein. Das ergibt 
zwei Gleichungspaare: 


(55) 


06) | +b(r) + 2c(r)- -a+d(na?) = 0. 


Aus diesen Gleichungen folgt, daß w eine ungerade Funktion vonz_ fi 
ist. Zur Bestimmung der fünf Funktionen «,, f,, b, d, C haben wir x 
jetzt vier Gleichungen, die beiden Differentialgleichungen und die == 
beiden für 2 N. Um das zu einem 
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bestimmten zu machen, setzen wir d = 0, so daB der endgiiltige An- 
satz für u und w lautet: 
(53a) + C(r)2?). 
(54a) w= (r)sin nz + 
Die Randbedingungen an den Endflächen fordern: 
(55 a) 2C(r)+ d'(r) =0 
(56a) D(-1P (Ada, + (24+ 2u)nB,) 
Und die Differentialgleichungen (A) und (B) nehmen die Form an: 
(A,) 9?(«,(r) cosnz + v?C(r)z? + (4 + 2u)(P ta, (r) cosnz + P&C(r)2?) 

+u(—n?e,(r) cos nz + 2C(r)) 

+ (2 + u) (np, (r)cosnz + b'(r)) = 
B,) v8, (r)sinnz + u(d PB, (r)sinnz + P b(r)2z) 

+ (2 + 2y)-(— (r) sin nz) 
+(A+ p)(—n Fa, (r)sinnz —22C(r)z) =0. 

Um diese Gleichungen zu lösen, entwickeln wir ihre linken Seiten 
in Fourierreihen und setzen jeden Fourierkoeffizienten einzeln gleich 
Null. Dadurch gewinnt man gewöhnliche Differentialgleichungen von 
der Form: 


a, (r)) + (4+ =0. 
v2 ß,(r) n+ — (A + 2u)t? (r) 
+ (# + Poe) 1 
++ = 


1) Physikalisch läßt sich diese Ergänzung so interpretieren: Die stehenden 
Wellen in dem Zylinder entstehen nicht durch einfache Reflexion an den End- 
flächen, sondern durch einen komplizierteren Vorgang, der durch die elasti- 
schen Eigenschaften des Materials und durch die Spannungsfreiheit am Rande 
bedingt ist. 

Die Knotenflächen der longitudinalen (w-)Schwingung sind nicht äqui- 
distant. Ihre Lagen werden durch die transzendente Gleichung 

o, sin nz2+C-z=0 

gegeben, also durch den Schnitt der;Sinuskurve mit einer schwach geneigten 
Geraden durch den Ursprung. Ihre Abstände werden nach den Enden des 
Stabes zu immer kleiner. In analoger Weise werden die Knotenflächen der 
transversalen (u-)Schwingungen durch die Schnitte der Kosinuskurve mit einer 
Parabel gegeben. Da b und C fiir den unendlich dünnen Stab (R = 0) verschwinden 
müssen, können sie auch für Zylinder endlicher Dicke als sehr klein angesehen 
werden, so daß die Neigung der Geraden und die Krümmung der Parabel im 
Bereich von — rn bis +2 ebenfalls sehr klein sind. Unsere Schwingungen 
unterscheiden sich nur von stehenden Wellen. 
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Gl. (B,) differenzieren wir nach r 

| +@+m(- tPHa,(r) —(— +P = 
und gehen mit den Abkürzungen = 


= 0. 
; 

0 in die Differentialgleichungen ein. Dann erhalten wir schließlich Ba 
c unter Berücksichtigung der Gl. (55a) 
3 a) e(r)+ (A+ Pe r)—prs 

sd Dies System kann gelöst werden unter der Voraussetzung, daB C ~ 


eine reguläre, gegebene Funktion von r ist. Obgleich dies nicht 2 ae 
zutrifft, wollen wir vorläufig C bekannt annehmen. Um =: ob 


A) Pe Pr) + (A+ 2u)P (r)— (A + 
(B,°) »? — + = 
die lautet: 


=0. (59) = AJ, (h,r) + BLY, (h, J, (k,r) + (k,"), 

den (60) = —A_J,(h,r) —B_Y,(k,r) J, + D,Y,(k,r) 

nd- mit 

jui- Das inhomogene Problem wird nach der Methode der Variation der 
Konstanten behandelt, indem A, B, E, D als Funktionen von r 
betrachtet werden mit den Nebenbedingungen: 

t 

(61) (hr) + Br) Y, (br) + 

(62) — A’(r)J, (hr) — Y (hr) +  (E(nJ, (kr) + D’(n)Y,(kr)) = 0. 
len Gehen wir mit den Ausdrücken Pa 
re (59) = (hr) + (hr) + Bes (kr) + DY, (kr), 
en 


(60’) (r) = — A(r) J, (hr) — Br) Y, (hr) + = (r) J (kr) + DY 
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in die Differentialgleichungen A,, B, ein, so erhalten wir Ei 5 

A { A’ (r) J, (hr) + Bn) + E (r) J,’ (kr) + D'(r) Y,' (kr) 

(Ay) = (- 

A’ (r)J (hr) + Bin) (hr) — = (E' (r) J,’ (kr) + D'(r) Y,’ (kr) 
= 2p) Cin). 


Diese Gleichungen bestimmen zusammen mit den Nebenbedingungen 
(61), (62) die Funktionen A, B, E, D. So wird z.B. 


(— +1.4.(4+2u)C(r) Y,'hr) - en - 5 Y,'(kr) | 
Y, (ar) J, (kr) Y, (kr) 
ke 
Y,hr) - - = Y, kr) 


J, (hr) Y, (hr) J,’ (kr) Y, (kr) 
J, (hr) Br u Y, (hr) J, (kr) Y, (kr) | 
| J, (hr) - - 


Die darin auftretenden Determinanten lassen sich wesentlich ver- 
einfachen. Man erhält: 
2 
. 24) Yan. Cr) 

A’ (r) = (— ti 


(& + (hr) Y’ (hr) - J, ar) 


Der J,’ (hr) Y, (hr) — J, (hr) (hr) im Nenner hat den Wert 
u = ; also bleibt der Bruch für r = 0 endlich, unter der Voraus- 


sl daß C regulär ist. Daher existiert sicher — übrigens auch, 
wenn C logarithmisch unendlich wird — 


Ebenso bestimmen sich die übrigen Koeffizienten. Es wird dann 


Y, hr)}ode 


+L, f C(e){¥, ke) J, (kr) + J, ko) Y, (kn)}ode 


0 
+ M, J, (hr) + M, Y, (hr) + +N, Y, (er). 
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Dabei miissen M, und N, verschwinden, um die Endlichkeit fir a 
r= 0 zu garantieren. Die Funktionen @ und # werden dann aus - Ber - 
den GI. (57), (58) bestimmt. Diese Ausdrücke haben nur dann einen Er 
Sinn, wenn C bekannt ist. Dies ist a priori nicht der Fall, Indessen > 
haben wir in der Gl. (56a) eine Beziehung zwischen a, f, b und C. Er: 
Um b zu eliminieren, differenzieren wir die Gleichung nach r un 
erhalten unter Berücksichtigung von (55) 
(56 b) 
Nun können die hier auftretenden Differentialausdrücke P ta, P 
durch wiederholte Anwendung der Gleichungen A,, B, linear durch “a 
&, und C ausgedrückt werden: 
Pita, = +(— PsC(n. 


= b, PF (a, + a, 4, + a, C(r)) 
Dann nimmt Gl. (56b) folgende Form an: a “i 
(56¢) 


0 


Gleichungen dieser Art heißen Volterrasche 
Der Kern dieser Gleichungen ist im vorliegenden Falle ausgeartet et 
und symmetrisch, also ist die Gleichung im Prinzip lésbar. ae o 
Der hier angegebene Weg zur Lösung des Problems ist AR" 
nicht befriedigend, und die rechnerische Durchführung stößt auf 
große Schwierigkeiten. Doch scheint dadurch wenigstens die Existenz 
der Lösung bewiesen. 
Es bleibt noch zu zeigen, da8 die Randbedingungen am Mantel 


64) 8 (R)cosnz + Ab(R) 
+ 2ua,'(R)cosnz + 2uC(R)z?*=0. 

(65) —nae,(R)sinnz + 2z2C(R) + £,'(R)sinnz + b'(R)-z=0. 
Die angewandte Methode wird auch auf diesen Fall übertragen: — a 
die linken Seiten werden nach Fourier entwickelt und die Koeffi- ; 
zienten einzeln Null gesetzt. Unter Berücksichtigung der Ab- — rs 
kürzungen (57), (58) erhalten wir N 
8, (B) + 0B,(R)) + 0. 


— «,(R) + 7,(R) = 0. 
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Mit (59), (60) und den Nebenbedingungen (61), (62) reduziert sich 
dies auf eine reine Bestimmung von Konstanten, die ganz ähnlich 
erfolgt wie bei Pochhammer. 


. Zur kopplungstheoretischen Interpretation 
‘der elastischen Schwingungen 


a) Kraftkopplung und Beschleunigungskopplung 


Im folgenden soll untersucht werden, welche Schlüsse aus 
unseren Überlegungen im Hinblick auf die Giebesche Interpreta- 
tion der elastischen Schwingungen als Schwingungen „gekoppelter 
Systeme“ gezogen werden können. Es konnte (S. 109) bewiesen 
werden, daß die Schwingungsgleichungen eines Hohlzylinders auf- 
gefaßt werden können als Differentialgleichungen gekoppelter Sy- 
steme; es konnte aber nicht eindeutig bestimmt werden, welches die 
schwingenden Grundgebilde sind und ob es sich um „Kraft“- oder 
„Beschleunigungskopplung“ handelt. Es gelingt vielmehr, aus den 
Gl. (1), (II) auf S. 107 zwei neue Gl. (I), (IT) zu machen, so daß 
beide Systeme äquivalent sind — sie entsprechen übrigens beide der 
Gleichung vierter Ordnung: 

Po E E? 

Das erste System entspricht einer Kraftkopplung, das zweite einer 
Beschleunigungskopplung. 

Es liegt die Frage nahe, ob man, wenn die Grundgebilde nicht 
eindeutig gegeben sind, Kraft- und Beschleunigungskopplung ein- 
deutig unterscheiden kann. 

Gegeben sei ein System von zwei Freiheitsgraden mit der kine- 
tischen Energie 


T = 2a,, 
und der potentiellen 
7 1 
= Cn +c 


Es liegt also Beschleunigungskopplung vor. Die Bewegungs- 
gleichungen für das System lauten: 


ax 
m | 2 +¢,%,+4 

| xy 


V 
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ae _ Ge _ 4, 
11 22 


die Kopplungskoeffizienten. Ihr Produkt k,,-k,, = q? bezeichnen we 
wir als Kopplungsfaktor. uns 


Durch Auflösung nach om und 3 können die Glei- % 


chungen (IV) umgeformt werden in das System: 


, d? 
a1, + %1 + Cig = 0, 


(V) 
| + Cr + 0% = 
Diese ‘allie ein System mit Kraftkopplung dar. Es ist dabei, ue 
wenn die Determinante a,, a,, — mit D bezeichnet wird: 
D Cy @ 
a= —, Cio = — An, 
Coe Cee 
Coq A 
- el, Ap, 
Po," @ 


Der Zusammenhang zwischen den alten Teilfrequenzen und ; 
Kopplungskoeffizienten und den neuen ist gegeben durch: x? 


(69) @, PP= Pio Py = Ky ky = 


Der Kopplungsfaktor ist also eine Invariante, während die 
Teilfrequenzen andere sind als früher. Theoretisch bedeutet der 
Ubergang vom System (IV) zum System (V) eine kanonische Trans- _ 
formatior. Dies sieht man leicht ein, wenn man an Stelle der Ge- 
schwindigkeiten die Impulse in die Energiefunktion einführt. Sie 
sind für das a (IV) gegeben durch: 


= a,, 4, 


Damit a die Hamiltonfanktion H: 


= 
| 


11 2 
00) 2H= p24 Ze pm C, 


Die kanonischen bleiben invariant, wenn man 
durch eine Berührungstransformation neue Koordinaten und Impulse _ 
Xu einführt. 
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Cy 
oder nach (68): 


Die Hamiltonschen Gleichungen lauten: 
aH 


aH 
“OP, = X,; 


Das erste System liefert: 


(71) 


das zweite: 


ay, Xı + + X2 = 0, 


9X + Xo + = 
und das sind die Differentialgleichungen (V). i 
Um die allgemeine Theorie auf unseren Fall anwenden zu 
können, bemerken wir, daß die Systeme IV und V derselben Gleichung 
vierter Ordnung äquivalent sind: co» 
+o, daz @,? 
: wobei jetzt unter w, und w, die „ungestörten“ Teilfrequenzen bei 
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Man wähle als Transformationsfunktion: 4 
wird: 
amd die Hamiltonfunktion (70) geht über in: 
i 
ay, 
=— P N 
oder [unter Berücksichtigung von (71)]: 
. 
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Mit dem Ansatz x = A.eirt rer wir die Frequenzgleichung: 


@,? + 


v? 


(72) — 
mit der Lésung: 
(73) = 


+ + Vio? — + 49° @,” 
2(1 — 9°) 


(Longitudinalschwingungen mit der Frequenz w, = 


den Kreisring (‘Transversalschwingungen mitder Frequenz, = V =) 


und setzt Beschleunigungskopplung mit dem Koeffizienten q an. Als 
Wert von q ergibt sich die Poissonsche Konstante o, wie es den 
Differentialgleichungen (I’), entspricht. 

Den ursprünglichen Gl. (T), (II) entsprechen aber zwei andere 
Grundsysteme mit den Frequenzen 


e(l — 0%) o(l - 


(S. 109). Können diese Systeme auch physikalisch interpretiert werden ? u 
Setzt man die Eigenfrequenzen ®,’ und ®,’ in die Frequenzen- 5 
gleichung ein, so lautet diese: 
— (0, ?+ 0, 0, ?0,?(1 — = 0 
(p: Kraftkopplungskoeffizient). 
Ihre Lösungen sind, wenn ®,’ harmonisch in der Form o,’ = ER = 
angesetzt wird: 


Wir erhalten zwei Serien von Frequenzen. 
wachsendem k einem Grenzwert zu 

Y_ (max) = yl —p*, 
die zweite » , hat bei 1 = 00 ein Minimum 


Die erste strebt mit 


Y+(min) = 


Dieses Minimum ist also die Frequenz der Balielheingeng eines 
unendlich langen Rohres. Sie ist: 


Einen anderen Grenzfall erhält man für »,, wenn man sien Radius 
des ins wachsen läßt. Es wird dann: 


2 ar 
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ka,’ bedeutet also die longitudinale Eigenfrequenz einer an zwei Seiten 
begrenzten, unendlich dünnen ebenen Platte. 


1 LV g(l 0%) 


Damit ist die Bedeutung von w,’ und w,’ aufgeklärt und die Frage 
von S. 109 beantwortet. Zusammenfassend kann man sagen: Die 
Schwingungen eines sehr dünnen Hohlzylinders können als Kopplungs- 
erscheinung aufgefaßt werden; und zwar kann man als schwingende 
Grundgebilde ansehen: entweder den unendlich dünnen Stab und 
den symmetrisch schwingenden Kreisring. Dann herrscht Be- 
schleunigungskopplung (a = 0 und / = 0); oder das unendlich lange 
Rohr und die zweiseitig begrenzte ebene Platte (I—> oo und 
a —> oo). Dann herrscht Kraftkopplung. 


b) DieKopplungstheorie und die Schwingungen des Vollzylinders 

Wendet man ähnliche Überlegungen auf den Vollzylinder an, 

so gelangt man zu einem Widerspruch. Giebe und Blechschmidt 

betrachten auch die Schwingungen des Vollzylinders als Kopplungs- 

erscheinung. Als Grundschwingungen werden die radialen, symmetri- 

schen Schwingungen einer dünnen Kreisscheibe und die Longitudinal- 

schwingungen eines unendlich dünnen Stabes betrachtet. Die 
Eigenfrequenzen sind bzw. 

oO, = as E und V- E 

1 l 22 R e(1 — 


Die Frequenzgleichung des Gesamtsystems stimmt formal mit der 
Gl. (73) überein. Auch hier kann man die Grenzfälle untersuchen, 
wenn man R oder | unbegrenzt wachsen läßt. Es ergeben sich als 
Grenzfrequenzen: 
. @, - 
muß mit der Eigenfrequenz der stehenden Schwingungen 
einer unendlich ausgedehnten, zweiseitig begrenzten Platte von der 
Dicke / übereinstimmen, @ mit der Eigenfrequenz der Radial- 
chwingungen eines unendlich langen Rohres. Aus den Grund- 
gleichungen kann man beide leicht ausrechnen. Sie sind: 


/ E(l-o) ak E(1 — 
e(1 + a)(1 — 20) ’ e(1 + «)(1 — 20) ’ 
wobei der transzendenten Gleichung £J,'(&) + J, (§ = 0 
genügt. 
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Daraus folgen die Gleichungen: 
2a) (1 +0) = 


die nur dann gleichzeitig befriedigt werden können, wenn 
yi 


ist. Diese Beziehung geben auch Giebe und Blechschmidt an, 
aber nur als Näherungswert, während sie eigentlich streng erfüllt 
sein müßte, wenn Giebes Theorie richtig wäre. Daß sie streng 
erfüllt ist, ist unmöglich, da die Wurzeln zweier verschiedener 
transzendenter Gleichungen nicht in einer konstanten Beziehung 
zueinander stehen können. Die angegebene zahlenmäßige Über- 
einstimmung muß als zufällig angesehen werden. 

Die Frage, in welchem Umfang die kopplungstheoretischen 
Vorstellungen trotzdem aufrechterhalten werden können, läßt sich 
ohne explizite Kenntnis der Lösungen der elastischen Gl. (A) und (B) 
nicht beantworten. Sie scheinen mit diesen Gleichungen nicht ver- 
einbar zu sein, wie aus folgender Überlegung hervorgeht. 

Die Kopplungstheorie verlangt, daß die Gl. (A) und (B) etwa auf 
Form werden kénnen: 


Pout+k, = = (0, 


+ ky, =0. 
Denn hebt man in diesen Gleichungen die Kopplung auf, so er- 
geben sich die Differentialgleichungen fiir die Schwingungen der 
Kreisscheibe bzw. des diinnen Stabes. 

Nun kann man die Gl. (A) und (B) nicht auf diese Form 
bringen, sondern nur — in ziemlich kiinstlicher Weise — auf die 
Form: 


62 o(l+o) OT, _ 
at z 


bw 
Or’ 
Es bleiben also noch Differentialausdriicke in den nicht ver- 
schwindenden Spannungen T,, und T,, übrig. 

Setzt man übrigens in diesen Gleichungen die Spannungen T.,, 
und T,, gleich Null, so ergeben sich die Differentialgleichungen des 
Rayleighschen bzw. der Kreisscheibe. 
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Das wichtigste Argument gegen die Deutung von Giebe ist 
jedoch die früher bewiesene Tatsache, daß stetige, periodische und 
differenzierbare Fortsetzung der Lösung der Grundgleichungen über 
den Stab hinaus nicht möglich ist, daß also die n-te Oberschwingung 
des Stabes nicht identisch ist mit der Grundschwingung eines Stabes 
von der Länge //n. 

Zusammenfassend muß gesagt werden, daß die Anwendung 
kopplungstheoretischer Vorstellungen auf den schwingenden Voll- 
zylinder auf einen Widerspruch führt und außerdem in den all- 


gemeinen Gleichungen der Elastizitätstheorie keine Rechtfertigung 
findet. 


Meinem hochverehrten Lehrer, Herrn Prof. v. Laue, möchte 
ich auch an dieser Stelle für die Anregung zu dieser Untersuchung 
und für sein förderndes Interesse und seinen Rat herzlichst danken. 
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Widerstandsänderung von Metallen im Magnetfeld 129 


Über die von Metallen 
im Magnetfeld 


Serban 


§ 1. Einleitung 
auf den Widerstand hun Leiter hat und daß die un ; 
des Widerstandes, mindestens für kleine Werte des Magnetfeldes, _ 
quadratisch vom Magnetfeld abhängt. Qualitativ kann man dieses __ 
Ergebnis ganz leicht erklären: das Magnetfeld verursacht eine 
Krümmung der Bahn des Elektrons und dadurch eine Vergrößerung 
des Weges, den das Elektron zurücklegen muß, um von einem be- 
stimmten Punkte zu einem anderen durch das Gitter hindurch- 
zuwandern. Eine Vergrößerung des Weges bedeutet aber eine Zu- 
nahme der Anzahl der Zusammenstöße mit dem Gitter, also eine 
Zunahme des Widerstandes. Da der Effekt unabhängig vom Vor- _ 
zeichen des Magnetfeldes sein muß, enthält eine Entwicklung nach 
Potenzen der Feldstärke nur gerade Potenzen und für schwache = 
Magnetfelder kann man die Reihe nach dem quadratischen Glied _ 


abbrechen. 

Es ist nicht schwer, die Theorie in dieser Näherung quantitativ IE 
durchzuführen und den Koeffizienten zu berechnen. Unter Berück- SR 
sichtigung der Entartung des Elektronengases ist das Problem von A ze 
Sommerfeld!) behandelt worden. Diese Theorie führt aber zu er 
großen Schwierigkeiten, wenn man sie mit dem Experiment ver- AS 


gleicht; denn nach ihr sollte erstens ein zum Strom paralleles Magnet- _ 
feld keinen Einfluß auf den Widerstand haben; experimentell ist 
aber der Effekt im longitudinalen Feld durchaus von derselben 
Größenordnung wie im transversalen Feld. Zweitens ist der be- 
rechnete Koeffizient viel kleiner als der gemessene. Allerdings hängt 
der experimentelle Koeffizient, nach Kapitza?), stark von der Vor- — 
behandlung des Materials ab. Drittens haben die Untersuchungen _ 
von Kapitza?) gezeigt, daß für stärkere Magnetfelder statt des 

quadratischen Gesetzes ein lineares gilt. Die Neigung der Geraden, 3 


1) A. Sommerfeld, Ztschr. f. Phys. 48. 8.1.18. 
2) P. Kapitza, Proc. Roy. Soc. S. 292. 
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die dieses Gesetz darstellt in einem Diagramm, in dem die relative 
Widerstandsänderung als Funktion des Magnetfeldes aufgetragen 
wird, hängt nicht von der Vorbehandlung des Materials ab; ihre 
Lage dagegen ist stark davon abhängig. 

Um die Theorie mit dem Experiment in Einklang zu bringen, 
hat Frank!) die Sommerfeldsche Theorie unter strenger Berück- 
sichtigung des Magnetfeldes durchgeführt. Die ersten beiden von 
den oben erwähnten Schwierigkeiten bleiben dabei bestehen, aber 
für die Abhängigkeit vom Magnetfeld bekommt er ein komplizierteres 
Gesetz, das für kleine Felder in das quadratische übergeht und 
für große einer Sättigung entspricht. In der Nähe des Wendepunktes 
verläuft die Kurve ungefähr linear, und damit glaubte Frank das 
Ergebnis von Kapitza erklären zu können. 

Zur Lösung der ersten Schwierigkeit hat Bloch?) versucht, 
die Wirkung des Feldes auf den Spin heranzuziehen. Im Magnetfeld 
ist die zum Feld parallele Stellung des Spins bevorzugt; nach dem 
Pauliprinzip hat diese Spinorientierung einen Einfluß auf die Ge- 
schwindigkeitsverteilung der Elektronen, der sich in einer Änderung 
des Widerstandes widerspiegelt. Es ist klar, daß diese Änderung 
unabhängig von der Richtung des Magnetfeldes ist. Die Durch- 
rechnung ergibt das quadratische Gesetz, aber ebenfalls mit einem 
zu kleinen Koeffizienten. 

Das Problem wird auch in einer Arbeit von Peierls*) be- 
handelt. Er verfeinert das Sommerfeldsche Modell durch Berück- 
sichtigung der Gitterkräfte. Bekanntlich ist die Wellenfunktion eines 
Elektrons im Gitter [Bloch®)] eine modulierte ebene Welle, aber 
die Energie ist eine kompliziertere Funktion der Wellenzahlen als 
die eines freien Elektrons. Wesentlich für die Peierlssche Er- 
klärung der Widerstandsänderung ist die Richtungsabhängigkeit der 
Stoßzeit (Zeit zwischen Zusammenstößen des Elektrons mit dem 
Gitter. Aus dieser Anisotropie folgt nämlich die Ähnlichkeit des 
Effektes im longitudinalen und im transversalen Feld. Da Genaueres 
über diese Anisotropie nicht bekannt ist, lassen sich aus der 
Peierlsschen Theorie keine quantitativen Schlüsse ziehen. Um die 
richtige Größenordnung der Widerstandsänderung zu bekommen, 
ist man gezwungen, eine erstaunlich große Anisotropie zu postulieren. 
Rein qualitativ kann man die Abhängigkeit vom Feld ermitteln: 
für kleine Felder bekommt man das erwartete quadratische Gesetz, 


1) N. H. Frank, Ztschr. f. Phys. 64. S. 650. 1930. 
2) F. Bloch, Ztschr. f. Phys. 53. S. 216. 1929. 
3) R. Peierls, Ann. d. Phys. [5] 10. 8. 37. 1931. 
sehr. f. Phys. 52. S. 555. 1928. 


fi 
L 
n 
h 
i 
P - u 
I 
R 
8 
€ 
N 
( 
A he a 
| 
| 
| 
| 
| 


für große eine Sättigung. Nun hat man experimentell nur bei $ es 
Halbleitern eine Sättigung gefunden; das bei den gewöhnlichen 


Leitern gefundene lineare Gesetz folgt also nicht aus der Theorie 
von Peierls. 
Es gibt noch eine Wirkung des Magnetfeldes, die man bis jetzt 


noch nicht zur Erklärung der Widerstandsänderung herangezogen en 


hat: die Quantelung der Elektronenbahnen im Magnetfeld. Diese 
Quantelung spielt eine wesentliche Rolle in der Theorie des Dia- 


magnetismus [Landau'), Teller?)], da nach der Quantentheorie, 


im Gegensatz zur klassischen Theorie, die Suszeptibilität des Elek- 
tronengases von Null verschieden ist. Es ist daher zu vermuten, 


daß die Quantelung auch auf den Widerstand einen Einfluß hat _ 


und es ist der Zweck der vorliegenden Arbeit, diesen Einfluß zu 
untersuchen. 


Die Wirkung eines Magnetfeldes auf im Gitter gebundene 
Elektronen ist sehr schwer zu behandeln; Peierls®) hat aber ge- 


zeigt, daß man die Bindung, genau wie im feldfreien Fall, ganz grob 


durch Einführung einer scheinbaren Masse für das Elektron berück- Sr ae 


sichtigen kann. Für unsere Rechnung werden wir uns mit dieser 


Näherung begnügen; das Termspektrum und die Kigenfunktionen 


werden als identisch mit denjenigen freier Elektronen mit dieser 


scheinbaren Masse angesetzt — diese Vereinfachung der Eigen- 


funktionen ist zwar sehr weitgehend, sie läßt sich aber nicht leicht 


vermeiden. Um den Widerstand zu berechnen, brauchen wir noch 


einen Ansatz für die Wechselwirkung der Elektronen mit den 
elastischen Schwingungen des Gitters; für diesen werden wir eine 
potentielle Zusatzenergie der Elektronen annehmen, die an jeder 
Stelle im Gitter proportional ist der lokalen relativen Volumänderung, 
die von den Gitterschwingungen verursacht wird. Dieser einfache 
Ansatz ist nicht wesentlich von dem, der in der gewöhnlichen Leit- 
fähigkeitstheorie abgeleitet wird, verschieden. Durch Vergleich der 


beiden könnte man auch den Proportionalitätsfaktor ermitteln, aber 7 


fiir das folgende brauchen wir ihn nicht. - 
Ferner werden wir mit der Blochschen Annahme, daß die _ 
Gitterschwingungen sich im thermischen Gleichgewicht befinden, 


die Rechnung durchführen. Die unter dieser Annahme gewonnenen _ 


Ergebnisse stimmen mit dem Experiment gut überein und es 
sprechen auch theoretische Argumente dafür, daß der Peierlssche 


1) L. Landau, Ztschr. f. Phys. 64. S. 629. 1930. 
2) E. Teller, Ztschr. f. Phys. 67. 8.311. 1931. 
3) R. Peierls, Zt ‚8. 763. 1933. 
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Einwand gegen diese Annahme für ein reales Gitter nicht so be- 
rechtigt ist wie für ein ideales. 


Es ist endlich noch ausdrücklich zu betonen, daß in der vor- 
liegenden Untersuchung das Magnetfeld nicht als kleine Störung 
behandelt wird, sondern schon in nullter Näherung streng berück- 
sichtigt. Als Störung werden nur die elastischen Schwingungen des 
Gitters und das elektrische Feld betrachtet. 


§ 2. Die Widerstandsänderung im transversalen Magnetfeld 


Es ist sehr leicht, im Rahmen der klassischen Mechanik, die 
Bewegungsgleichungen eines Elektrons, das sich unter dem Einfluß 
eines elektrischen und eines, zum elektrischen senkrecht stehenden, 
magnetischen Feldes befindet, zu integrieren. Die Bewegung parallel 
zum Magnetfeld ist gleichförmig. Senkrecht dazu ist sie im ein- 
fachsten Fall eine gleichförmige Translation längs einer Geraden, 
die senkrecht zu beiden Feldern steht; im allgemeinen Fall über- 


a ‚lagert sich ihr noch eine Schwingung um diese Gerade als Achse. 


Jedenfalls ist der Zeitmittelwert der zum elektrischen Feld parallelen 
Geschwindigkeitskomponente gleich Null. Wenn wir dieses Ergebnis 
auf unser Problem anwenden, so folgt, daß das Elektronengas keine 
zum elektrischen Feld parallele Stromkomponente hat. Ist aber 
dieses Gas dem Einfluß der Wärmebewegung des Gitters unter- 
worfen, so ist es möglich, daß die Bahnachse jedes Elektrons eine 
säkulare Verschiebung erfährt und dadurch wird ein Strom parallel 
zum elektrischen Feld hervorgerufen. Daß das tatsächlich der Fall 
ist, kann man leicht einsehen, wenn man die störende Wirkung 
der Gitterschwingungen als eine Art Reibung betrachtet und die 
Reibungskraft als proportional und entgegengesetzt der Geschwin- 
digkeit annimmt. 

Wenn man den Proportionalitätsfaktor mit —a@ («>0) be- 
zeichnet, so heißen die Bewegungsgleichungen in demjenigen Koordi- 
natensystem, dessen z-Achse parallel zur elektrischen Feldstärke F 
und dessen z-Achse parallel zur magnetischen Feldstärke H steht 


, =— eF H—«0,, 


e 
=‘v,H—av,, 


Im folgenden werden wir von der Bewegung in der z-Richtung 
_absehen 
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S. Tifeica. 


durch +, = 6,=0. Man bekommt dann sofort 


eF-<H 
c 


eFa 
2 er 2 
a +H 
c 


’ y 


a +H 
c 


Um die allgemeine Lösung zu finden, machen wir den Ansatz 


und bekommen fiir A und B die Gleichungen 


4 = — — . — = — 
mo, -H—av,,, m 6, H-av), 


oder für die Größe u = v, +40,’ die Gleichung 


Die Lésung dieser letzten Gleichung ist en, 
a 


nl. (mit komplexem C\. 


Wegen des Dämpfungsfaktors sind die Zusatzglieder v,’ und v, für 
unsere Diskussion uninteressant, so daß wir uns auf die einfache 
Lösung beschränken können. 
Man sieht zuerst, daß jetzt die Geschwindigkeitskomponente v,, a 
also auch die ihr proportionale Stromkomponente s, von null ver- 
schieden sind. Der Widerstand o ist definiert durch 
wenn die Reibung als klein angenommen wird, so ist s,, weil durch © = ; 
die Reibung verursacht, viel kleiner als s,, und kann im Nenner — 
von (1) vernachlassigt werden. Dann wird dl 


die RN. die allgemein durch 


definiert gleich R =— ist. 


8, 


Nach der Quantentheorie läßt sich das Problem ohne Reibung Hs 
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_ gewonnenen, sehr ähnlich. Wenn man für das Vektorpotential den 


a Ansatz 4,=4,=0, A, = H.x macht, so wird die Hamiltonfunktion 
des Problems: 


1 H _\2 1 
W= + (p, + 2) + +eFx 


wo m die scheinbare Masse und — e die Ladung des Elektrons ist. 
Man sieht sofort, daß p, und p, mit W vertauschbar sind, so daß 
Fa diese Größen als Zahlen betrachtet werden können, wenn man die 


Eigenfunktionen in der Form f(z)e * schreibt. Man kann 
dann W folgendermaßen umformen: 


2, Mm Py 
(=+ 


e? Ft 
mo, M Wo” m 2m a,” 


2p." mM @,2 ler 
+ 
_ dabei ist w, die doppelte Larmor-(Kreis-)frequenz und 


eF 

0 Mo, 
AR Aus dieser Form von W sieht man sofort, daß derjenige Faktor 
der Eigenfunktion, der nur von & abhängt, die bekannte Eigen- 
_ funktion g(x — x,) des Oszillators mit der Masse m, der Frequenz ®, 


Piss und der Gleichgewichtslage z, ist; wenn man jetzt p, durch , und p, 
fe: durch die entsprechende Wellenzahl / = fe ausdriickt, so werden 
Eigenfunktionen 


‚eF 


und die entsprechenden Energiewerte 
Ey, 2,.1 = hw, (n + 3) + +eF. 
_ wenn man von der additiven Konstante = 
0 

&, ist der Mittelwert von 2; der Mittelwert der A: 
der Geschwindigkeit ist daher der Zeitdifferentialquotient von z,. 
ia Wenn keine Störung vorhanden ist, so ist v,= 0. Es ist auch in 
diesem Fall zu vermuten, und das werden wir später genau verfolgen, 
er unter dem Einfluß der Wärmebewegung des Gitters x, nicht 
3 konstant bleibt, sondern sich irgendwie mit der Zeit ändert. 
> ae Den Strom in der z-Richtung kann man dann direkt aus seiner 
a berechnen: man lege senkrecht zur z-Achse eine Wand; 
dann ist diese Komponente des Stromes gleich dem Überschuß der 
Be Anzahl der Elektronen, die pro Zeiteinheit vom negativen zum posi- 
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positiven zum negativen Halbraum übergehen, multipliziert mit der 
Ladung des Elektrons. Die Stromdichte bekommt man durch Division 
durch den Flächeninhalt der Wand. 

Zur Berechnung des Widerstandes braucht man auch den Strom 


in der y-Richtung. Da v, = = (p, + <7) ist, so bekommt man 
durch Mittelwertbildung 
1 eF 


ein Resultat, das genau dem klassischen Wert der mittleren Geschwin- 
digkeit senkrecht zu beiden Feldern entspricht. Man sieht also, daß 
aus unserer summarischen Berücksichtigung der Bindung nur der 
normale Hall-Effekt folgt. 

Die Stromdichte in der y-Richtung ist 


(6) 


MC Wy 


wobei » die Anzahl der Elektronen pro Volumeinheit bedeutet. 


§ 3. Wechselwirkung mit den Gitterschwingungen 


Wie in der Einleitung geschildert wurde, werden wir jetzt den 
Einfluß einer Störung untersuchen, deren Energie proportional der 
lokalen kubischen Dilatation ist 

> 
wobei u die Verschiebung eines Gitterpunktes und D der Proportio- 
nalitätsfaktor ist. 

Um die Berechnung der statistischen Gewichte zu erleichtern 
ist es vorteilhaft, einen Wiirfel mit zu den Koordinatenachsen 
parallelen Kanten und mit der Kantenlänge L als Grundbereich 
einzuführen und von allen räumlichen Eigenschaften des Systems 
Periodizität mit der Periode’L in allen Koordinatenrichtungen zu 
fordern. Für die Eigenfunktionen des Elektrons läßt sich die 
Forderung der Periodizität in der z-Richtung nicht streng erfüllen; 
wenn man aber bedenkt, daß die Oszillatoreigenfunktionen im Un- 
endlichen exponentiell verschwinden, kann man immer L so groß 
wählen, daß es größer wird als das Intervall, in dem diese Eigen- 
funktionen wesentlich von Null verschieden sind, und dann kann 
man ganz leicht die Eigenfunktionen periodisch fortsetzen. Wegen 
des exponentiellen Verschwindens kann man auch die Integration 
der Eigenfunktion über alle Werte von x im Grundintervall durch 
ein Integration von — oo bis + 00 ersetzen. 
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Aus der Periodizitätsforderung folgt, daß die (Guastonnahlen 
und x, diskreter Werte fähig sind: 


Die Wirmebewegung des Gitters zerlegen wir in fortlaufende 
ebene Wellen. Da nur die longitudinalen Wellen eine Dilatation 
verursachen, so werden wir uns auf diese beschränken: 


Er; 


> 


wobei 5 f ein ganzzahliges Vektor, und f = |f| ist. Die b's, wie 


aus der Theorie des festen Körpers wohlbekannt, sind folgende 
Operatoren: 


b> = YN> +1 4%, br sire 4>, 
wobei M die Masse des Kristallwürfels mit der Kante L bedeutet, 
N> die Besetzungszahl der betreffenden Eigenschwingung und A> 
bzw. 4> diejenigen Operatoren sind, die N> in Nr +1 bzw. I-95 
transformieren. 


Im folgenden werden wir die Dispersion der elastischen Welle 
vernachlässigen und mit einer konstanten Schallgeschwindigkeit w 


die Rechnung durchführen, so daß DIN 
(8) wf. 


Unter dem Einfluß der Wärmebewegung des Gitters werden die 
Elektronen Übergänge von Zustand zu Zustand ausführen. Da die 
Störung linear in den b’s ist, so ist ein Übergang des Elektrons mit 
einem Sprung der Quantenzahl irgendeiner Eigenschwingung des 
Gitters um eine Einheit verknüpft. Wegen der Analogie mit der 
Strahlungstheorie werden wir diese Elementarprozesse „Absorption“ 
und „Emission“ eines Schallquants nennen. 

Wir werden noch für die Quantenzahlen des Elektrons die 
kürzere Bezeichnung gq einführen. Dann lassen sich in bekannter 
Weise, durch Anwendung der Methode der Variation der Kon- 
stanten, die Wahrscheinlichkeiten der Elementarprozesse leicht be- 
rechnen. Die Übergangswahrscheinlichkeit pro Zeiteinheit für den 
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(9b) 


Das Zeichen (®),-, bedeutet das Matrixelement der ‚Funktion, die 7 


in der Klammer steht: eles 


integriert über den Grundbereich. 
Daraus bekommt man sofort die Auswahlregeln 


m 
für eine und 
h 
(10b) =I—f, Lo = + 
fiir eine Emission. 
Wenn diese Auswahlregeln erfüllt sind, dann ist das 
element von Null verschieden und gleich 


+00 


(2 — Lo) Pn’ (4 — 2’) )dz 


+00 
bzw. (@— 2,) Pu — 2, )dx 

fiir eine Absorption bzw. Emission. Das Quadrat des absoluten Be- Fs 
trages ist in beiden Fällen dasselbe und wird mit Wy(f,, %, %) 
bezeichnet. Diese Funktionen werden später eine große Rolle spielen, 
darum seien hier die wichtigsten nents derselben näher 
untersucht. 

Wir gehen vom 


jua 
Sorte f ist gegeben durch = ees 
1 — cos (=" 0) t 
’ tMw E,-E 2 f 
h 
de fiir denselben Ubergang und die Emission eines Schallquants: ra a 
— cos | —— + @ 4 
1e { 
. 
atrix- ‘ 
ie 
ie 
it 
r 
“ € 
r 
n aus, so daß Way ist. Wir führen die neue Integrations- _ 
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4 
variable = vr x und die Funktion ®, (tl) = ein. 


Die ®, (t) sind aay rg normierten Eigenfunktionen des Problems 


Für sie gibt es bekanntlich eine erzeugende F vn 


1? 


Mit den neuen Bezeichnungen wird N 


_ miteinander, dann das Produkt mite " “> und integrieren über t 
—oo bis +00, so bekommen wir eine erzeugende Funktion für 


iA 
1 ma 


Daraus folgt, daß die Jan gleich der TE 


2)” 


1 fx (20 + %') 


multipliziert mit einem Polynom in 


(n 
in 
zu 
; | D, i— i 0 7 h ( 
vy = \ 
ti 
| Vx 2" m 
= n! n / 
0 
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sind, und daB W,, „= |J„ „|? nur von der Größe 


(nach den Auswahlregeln) patel Es genügt also, um die W, 
in voller Allgemeinheit zu haben, die J„„ nur in dem Spezialfall 
M Gp 


zu untersuchen, in dem f, = 0 und /* I — %) =a ist. Man 


hat dann nur zu quadrieren und «? durch — (f,’+f,?) zu ersetzen, 


um die W„„ zu bekommen. In aE Spezialfall wird die er- 
zeugende Funktion: 


Jn n(@)U" 0”. 


Man verifiziert leicht, daß F Differential- 
te ae t(1-)F 


a 


or or or 1 oF OF 
2 m»? : . 
dv? 2u0 ou | + | 1) (u Ou Ov 


a? 
Setzen wir die Potenzreihe für F in diese Gleichung ein und 
vergleichen wir die Koeffizienten ‘gleicher Potenzen, so bekommen 
wir für J, „- folgende Differentialgleichung: 


da? da a? 


Später werden wir Summen von Ausdrücken, in denen die W’s 
vorkommen, durch Integrale ersetzen. Dazu ist es zweckmäßig, 
die W’s mit Hilfe der Kramerschen Methode aus dieser Diffe- 
rentialgleichung zu berechnen. Das Quadrat des oszillierenden 
Faktors werden wir durch seinen Mittelwert ersetzen; man be- — 


kommt dann 


(12) 


7 a* 
(n +n’ + a? — (n — mM)’ — 


Man kann die Normierung feststellen, wenn man von der Darstellung 


+00 


V = Wi (a? + y’) = | fan + y) ®,, (s)ds 


+00 +00 


f + (t + y) (8) 
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ausgeht. Es folgt: 
+00 +00 


f (a? + y*)dax dy =2 (02) «de 


= ®,(+Y) ()dsdtdaxdy. 


Die Integration über x ergibt, mit der Diracschen Symbolik, 
2nö(s —t); dann kann man die Integration über ¢ leicht ausführen 
und man bekommt fiir die rechte Seite von (12a): 


Wegen der Normierung der ®’s ergibt die Integration über y 
und dann die letzte Integration auch eins; es bleibt also: 


[Wan =1. 
0 


Nun ist der Zahlenfaktor von (12) auch so eingerichtet, daß 
1 ade 


at 
(n +n’ + 1) a?— (n — — 


wobei die Integrationsgrenzen die positiven Nullstellen des Nenners 
sind. 


= 1, 


$4. Die stationäre Verteilungsfunktion 


Um den Strom zu berechnen, muß man zuerst die stationäre 
Verteilungsfunktion ermitteln. Es läßt sich leicht zeigen, daß die 


1 


Fermifunktion 7, = „5, —— , wobei «, die Energie ohne das Glied 


mit dem elektrischen Feld ist: 


1 


dieser Forderung genügt in der Näherung, in der man nur lineare 
Glieder in der Entwicklung nach Potenzen der Feldstärke F be- 
rücksichtigt. 

Um das zu beweisen, muß man die Zeitableitung von 7, be- 
rechnen. /, ist gleich der Anzahl der Elektronen, die pro Zeit- 
einheit in den Zustand q springen minus der Anzahl derjenigen, 
die aus diesem Zustand springen. 


- 
(12a 
= 
- 
2a (8 + ¥) B, (8 + y)@,, 
— 
D 
f- 
B 
| 
Bes, 
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Bet: 
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q 
f 
1 — cos (=; o)t 
E.- E 2 
h 
av 1 - cos o)t] 
| h o) 


Daß die rechte Seite für F = 0 verschwindet, folgt aus der 
bekannten Stationarität der Fermifunktion für diesen Spezialfall. 
Die rechte Seite hängt aber nur durch die Ausdrücke E,— E, von 
F ab, also kommt F nur in der Kombination eF(z,'— %,) vor. 
Das Glied erster Ordnung bei der Entwicklung nach Potenzen von 
F ist also das Produkt dieser Kombination mit einer geraden 
Funktion von (2, — z,. Da das Ganze noch über x,’ summiert 
wird, so verschwindet dieses Glied der Entwicklung und unsere 
Behauptung ist bewiesen. Wichtig ist, daß über die Temperatur 
keine Annahme gemacht worden ist. 


§ 5. Der Strom 


Jetzt gehen wir zur Berechnung der z-Komponente des Stromes 
über, nach dem Verfahren, das oben (§ 2) besprochen wurde. Als 
„Wand“ wählen wir die Ebene z = 0. 


| 


- —2,) 2,1 — x,)(N 


1 — cos o)t 


h 


| Wir eliminieren zuerst f, und f, mit Hilfe der Auswahlregeln, 
und dann ersetzen wir die Summationen über /,/' 2,2,’ und f, 
durch Integrationen. Das Gewicht jedes der Integrationselemente 
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dl, d ist gleich L/2r, das Gew von dz, und dz,’ ist 
Mo, 
[nach (7)]. 
e D? 
Mhw 


falar —z,)N 


1<0,2'>0 
1 — cos (75% - o)t 
z,( —2Z)(N + 1)] 


bekommt man 


1 e D? 
2 


dvda,da, f W,,,, (1 — x4) N 


Wir führen zuerst die Integration über u durch. Dabei be- 
1 — cos 2t 
$2? 

_ wir in die anderen Faktoren des Integranden für w die Lösung der 
Gleichung einsetzen und die Brüche allein integrieren. Beide 


nutzen wir den ö-Funktionscharakter der Brüche ‚indem 


Glieder ergeben dasselbe Resultat: = Soe Wenn man es 


in die Formel (15) einführt und noch die Differentiation nach t 
durchführt, bekommt man 
e D? (3x) f dvdz (A + B), 


u h Mhw \2n h 
n,n’ 


x 
(1 
SE 
6 
(1 
st 
E,- E, 
P 1 — cos i +olt 
4 —-z\N1- : . rt 
Durch Einführung der Integrationsvariablen wu=1+U und § (1 
ü 
5 / 2 2 
— %,(1 — x,)(N + 1) - 
SC 
al 
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ie fi 
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wobei 


aw —%_) N — x, (1 — x) (N +1], 


Ww lt, N + 1) —z,(1 — x,) 
In A ist für u die Lösung der Gleichung 
(17a) E„-E,-ho=0 
einzusetzen; in B diejenige der Gleichung 
(17b) 
Jetzt entwickeln wir S, nach Potenzen der elektrischen F eld- 
stärke F. 


Das Glied nullter Ordnung verschwindet, da in dieser Nähe- 
rung A=B=0. Für F=0 gehen nämlich (17a) und (17b) in 


(18a) &,—&,—ho=0, 
(18b) —& tha =0 


über; de z, die Fermifunktion der Variablen «, int aa N die 
Plancksche Funktion in Aw/kT ist, so kann man das Verschwinden 
von 4 und B mit Hilfe von (18a) und (18b) leicht verifizieren. 
Um das Glied erster Ordnung zu bekommen, ersetzen wir in 
N ho durch ¢,,— ¢,+ eF (2, — bzw. — + eF — x,) und 
entwickeln es nach eF'(z,’— x,). Die ersten Glieder von A bzw. B 
werden 


(18a) nn kT 


bzw. 


(18h) B, = FINN + 


Die Faktoren N (N + 1)(y,,—z,) hängen, durch u, noch von F ab, 
wenn man sich aber auf die Glieder erster Ordnung in F be- 
schränkt, braucht man nur in diesen Faktoren F = 0 zu setzen, 
also u als Lösung der Gl. (18a) bzw. (18b) zu betrachten, was wir 
von jetzt an tun werden. 

Setzt man jetzt A, und B, in die Formel für S, ein und 
führt man noch die neuen Variablen r= a2, — x, und s = 2,'+ 2, 
ein, so bekommt man: 


ne e D? L\5/mo 
Die Integrationsgrenzen der ER über s sind —r und +r 
und diejenigen der Integration über r sind 0 und +00. Da der 
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Integrand nicht von s abhängt, so kann man leicht die Integration 
über s durchführen: 

eF eD L M Wy 


_ Der Integrand hängt nur von r? ab; da nur für positive Werte 
von r integriert wird, so kann man überall, nach der Auswahl- 


h 
regel (10), r = 2,’ — x, durch win f, ersetzen und über f, von 0 
bis + 00 integrieren: 


(a) 2 äj,dof (A —B). 


kT Mhw 2a 


Führen wir noch in der f,, f,-Ebene RETTEN ein: 


f,=f cst, f,=f snd 


und integrieren über # von 0 bis a, so wird 


> far dv f'3(4,— B,). 


x 2 kT Mhw \2n M ©, 


Der Integrand ist eine gerade Funktion von v, man kann also als 
Integrationsgrenzen 0 und +00 wählen und das Ganze mit 2 
multiplizieren. Dann kann man für |v| einfach f, schreiben (nach 
den Auswahlregeln): i 


kT \?2n) mao, 


Fiir die weitere Rechnung brauchen wir noch die Ausdriicke 
von e, und &,, aus denen / und l’ eliminiert sind mit Hilfe der Aus- 
wahlregeln und des Energiesatzes. 

Für A, bekommt man 


und für B, 

-n+n+1)+ + 57, + (n’ — n)]? 


Führen wir die Abkürzungen 


1 ho, 


ein 


(21 
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und 


1 


D er ihe 
> far af, Wow N(N+1) 


ein, so wird der Strom 


9 n,n’ 
(21) 1 1 1 
e “+1 e rae e “+1 e 
wobei ö die Dichte des Materials bedeutet: ö= nz . 
eckige Klammer im Integranden werden wir auch die Abkürzung K 
einführen: 
1 
x-3 X+5 2-7 X+5 sr 
e “+1 e +1 e +1 e 


Jetzt können wir den Widerstand ganz leicht berechnen. Definitions- 
gemäß ist der Widerstand gleich dem Verhältnis zwischen der 
Komponente der elektrischen Feldstärke in der Richtung des 
Stromes und dem Betrag der Stromdichte. 
Wenn man mit s, und s, die Komponenten der Stromdichte 
bezeichnet, so ist der Widerstand 
(23) 
s, rührt nur vom Magnetfeld her und ist von Null verschieden 
auch wenn keine Wärmebewegung des Gitters vorhanden ist'); 
s, dagegen ist von dieser Wärmebewegung verursacht. Wenn die 
störende Wirkung der Gitterschwingungen als klein betrachtet 
werden kann, dann ist s, wesentlich größer als s, und man kann 
diese letzte Größe im Nenner von (23) vernachlässigen. Man be- 
kommt die Komponente s, der Stromdichte durch Division der 
Stromkomponente S, aus der Formel (21) durch L?, den Querschnitt 
des Metallstiickes. s, haben wir schon früher berechnet (6); wenn 
wir jetzt diese Werte in (23) einsetzen und die oben besprochene 
Vernachlässigung nce so bekommen wir: 


‘ m?’ 
24) = df. Wee N(N+1)K. 


§ 6. Diskussion des Widerstandes 


Um den Ausdruck (24) fiir den Widerstand durchsichtiger zu 
gestalten, muß man ein Näherungsverfahren anwenden, da die 


1) Vgl. § 2. 
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Summation über n und n’ sich nicht streng durchführen läßt. Für 
kleine Magnetfelder kann man in gröbster Näherung diese Summation 
durch eine Integration ersetzen. Die W,,, sind aber nur für ganz- 
zahlige n und n’ definiert, daher muß man sie zuerst, um die 
Integration durchführen zu können, durch die angenäherten Aus- 
drücke (12) ersetzen. Wir schicken noch die Bemerkung voraus, 
daß das Polynom unter dem Wurzelzeichen von (12) nur dann 
positiv zwischen den Nullstellen ist, wenn n +4 und n’+ } positiv 
sind. Man kann daher ruhig als Integrationsgrenzen diese Null- 
stellen wählen, weil sie immer im Quadranten n +4>0, n’+1>0 
liegen. 

Wir führen zuerst die neuen Variablen p=n+n’+1 und 
q=n’—n ein. Dann haben wir ae 


K. 


wi ir integrieren zuerst über p. Die BER ist ganz leicht 
durchzuführen, wenn man die Fermifunktion durch eine Treppen- 
funktion ersetzt. Dann ist die Funktion 


(26) 


5 = leich 1 im Intervall — = .< X< re und sonst Null; da aber X 
7 

eine lineare Funktion von p ist, so ist (26) als Funktion von p be- 
_ trachtet gleich 1 in einem Intervall p, <p<p, und sonst Null. 
Dann ist das Den. über den ersten Summanden von K dr den 


und xX-4=0, 


wobei P eine für weitere, uninteressante Glieder ist. 
Der Koeffizient von unter der rechten Seite 


h A 


| 
Pr 
; 
= 
k 
: u 
k 
AT ma, 1 
0 
z 
b 
& 


4 


§. Titeica. Widerstandsänderung von Metallen im Magnetfeld 147 


von (27) ist — (1+ =— 


Wert (f,>+f,7) einsetzt. Die Integration über q jedes der 


m 


wenn man für @ seinen 


Summanden der rechten Seite von (27) ist äquivalent der Quadratur — 


einer halben Ellipse, deren eine Achse gleich der Wurzel des von gq __ 


freien Gliedes unter dem Wurzelzeichen des Summanden und deren 


andere um den Faktor f,/f größer ist; diese Integration ergibt also | 


das von q freie Glied multipliziert mit > ra - Die Differenz der © 


beiden Integrale über q ist dann gleich der Differenz der q-freien _ 
Glieder multipliziert mit demselben Faktor, also + 


aS 


Der zweite Summand von K ergibt nochmals dasselbe, so dab 


AT gh, 


Wenn man dieses Ergebnis in den Ausdruck (24) einsetzt, so be- ke a 


kommt man fiir den Widerstand 


D: 2 1 


Nach Einführung von Polarkoordinaten f, = f cos 9, f’ =f sin 3 


kann man ganz leicht die Integration iiber den Winkel ausfiihren, $ 


ho hw f 


und wenn man noch f durch §= 77, = —y ersetzt, so bee 


kT 
kommt man 


9/1 
1 5 
(29) = ren 


24 n® Oh? ws = e* 


Die obere Integrationsgrenze ist, in bekannter Weise, durch die 


charakteristische Temperatur @ des Metalls ausgedrückt. In der _ 

Näherung, in der die Quantelung zu vernachlässigen ist, hängt also = 
der Widerstand nicht vom Magnetfeld ab, so daß nach unserem © 
Modell die Feldabhängigkeit als typischer Quanteneffekt aufzufassen 
ist. Es ist noch zu bemerken, daß man genau die Formel (29) be- 


kommt, wenn man die Blochsche Rechnung für die Leitfähigkeit 
ohne Magnetfeld auf unser Modell überträgt- 


§ 7. Widerstandsänderung für kleine Magnetfelder 
Um nun als nächste Näherung die Feldabhängigkeit für kleine 
Magnetfelder zu bekommen, muß man die Summe über n und n’ 


besser aproximieren als wir es bis jetzt getan haben. Wenn man 
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vom Integral als nullter Näherung ausgeht, kann man durch Addition 
von Korrekturgliedern die Approximation verbessern. Bekanntlich 
gibt die Euler-Mac Laurinsche Summenformel eine Reihen- 
entwicklung für die Korrekturglieder. Wir werden hier nur das 
erste Glied dieser Reihe berechnen. Die Summenformel lautet dann 


a y) dady 
0 0 


n=0 n’=0 


1 1 
+ £200, + at f,@, O)dz. 
0 0 


In dieser Form läßt sich aber die Formel nicht auf unser Problem 
anwenden, da der Ausdruck (12) eine gute Näherung fir W,,,, ist, 
aber seine Ableitungen nach n und n’ keine guten Näherungen für 
die Differenzen der W mehr sind. Deswegen müssen wir mit den 
Differenzen der strengen Ausdrücke für die W rechnen, und statt 
längs des Randes des ersten Quadranten der n, n’-Ebene zu inte- 
grieren, müssen wir jetzt für ganzzahlige n und nm’ summieren. 
Wir brauchen nur die Werte am Rande und die Differenzen nach 
der inneren Normale, also die Ausdrücke von W,,9 und W,ı — Wyo 
und diese lassen sich ganz leicht aus unserer erzeugenden Funktion (11’) 
berechnen 


1 [a®\n-1/a! 
-W,o=(n+ 1)(Wa+1,0 — Wa) — n—1,0)+ 


ae > Wenn « einigermaßen groß ist, hat W,„,.. als Funktion von n 
betrachtet, ein scharfes Maximum für diejenige ganze Zahl, die am 
nächsten an «?/2 liegt. Es ist also eine erlaubte Näherung, um 
eine Summe 


zu berechnen, wobei g nicht allzu schnell veränderlich ist, diese 
Funktion in der Umgebung der Stelle n = = in eine Potenzreihe 
zu entwickeln und nur das erste Glied beizubehalten: 


| 
= = 
5 
ic 
i 
7 
| 
| us 
IM 
i 
at 
] 
DS 
h 
; ‘ 
( 
| 
( 
i oo 
h 
2 ‘ 
N 


§. Titeica. Widerstandsänderung von Metallen im Magnetfeld 149 Re 


In derselben Weise folgt 
> War — Wao) pn) = 


~- > — Wn -1,0)' (n) 


‘ ‘ 1 
symmetrisch in n und m’ ist, so geben die Ränder n +, =0 


und n’ + 5 denselben Beitrag: 


oo 


Ke [War K dndn’ + 


—1/2 -ı/2 n 


co 


109 
+ 2 > — W,,0)(Kw =0- 


n=( 


‘ 


Das Integral ist im vorigen Paragraphen berechnet worden. Wenden 
wir unsere obige Rechenregel auf die Summation über n an, so be- 
kommen wir für das Korrekturglied folgenden Ausdruck: “8 


- 


n’=0 ni =0 n’ =0 


Pet 


Wenn man mit 7 und 7” die Ableitungen der Fermifunktion © ® 
1=— - : nach dem Argument u bezeichnet und mit X bzw. X 
e 


die Werte von X bzw. X’ fürn = 


so ist der obige Ausdruck (30) gleich 
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Um das Korrekturglied des Widerstandes zu bekommen, muß man 
zuerst diesen Ausdruck mit . [Nn (N + 1) multiplizieren und dann 
über f’ und f, von 0 bis oo FRE Es ist aber zweckmäßig 


als Integrationsvariable f und tof einzuführen. Man hat dann 

mit = f*N(N + 1) 

und über ¢ von 1 bis + 00 zu integrieren. Wir führen zuerst die 
Integration über ¢ durch. Die Argumente der Ableitungen der Fermi- 
funktion sind linear in ¢ Bei der Integration benutzen wir die 
Eigenschaft der negativen Ableitung der Fermifunktion, daß sie ein 
scharfes Maximum für den Wert 0 des Arguments hat. Man kann 
daher die anderen Faktoren als langsam veränderlich betrachten 
und in Potenzreihen um diese Stelle entwickeln. Die Integration 
der Glieder, die die zweite Ableitung als Faktor enthalten, kann 
durch eine partielle Integration anf den ersten Fall zurückgeführt 
werden. Das Ergebnis der Rechnung ist 


zu multiplizieren, über f von 0 bis + oo 


- 


mw\? \t,? 
2 1 
| 


wobei t, bis t, die Lösungen der vier Gleichungen X-; = (0 


X + 4 = 0 sind, also 


der Maximalwert von Awf gleich k@ ist, so ist 2 — sehr klein 
0 


fiir alle in Frage kommenden Werte von f. Man kann daher (31) 
nach Potenzen dieser Größe entwickeln und man bekommt 


~ 
+4 


und über f zu integrieren. Führen wir = §= 


_ grationsvariable ein und multiplizieren wir mit den Vorfaktoren 
Ne von (24), so bekommen wir die Widerstandsinderung 
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O/T 
nan 
D: mot 

Big (32) oT 

foe) 0 

di und durch Division mit dem Widerstand für H = 0 ch die relative my = 

Widerstandsänderung 

li- 

= wobei 

ten 

hrt § 8. Der Grenzfall starker Magnetfelder 

Die Bedingung für die Anwendbarkeit der im vorigen Para- — 
graphen diskutierten Näherung ist ha,<kT. Man kann leicht einen — 

b anderen Fall diskutieren: es ist der Fall, in dem das Magnetfeld so 
stark ist, daß alle Elektronen sich in den Zuständen befinden, deren 
Quantenzahl n gleich Null ist. Man kann sich sofort überzeugen, 

* daß in diesem Fall unser Ausdruck (24) für den Widerstand linear 


vom Magnetfeld abhängt. Von der Summe über n und n’ hat man ~~ 

nämlich hier nur das erste Glied, mit n = n’ = 0, zu berücksichtigen. 

Die Fermifunktionen hängen nicht mehr vom Magnetfeld ab, so daß 
A 


ur der Faktor vor dem Integral von (24) und Woy = e tna,” av 
von H abhängen. Der Faktor ist proportional H und Woo kann 
man in eine Potenzreihe entwickeln; wenn der Maximalwert des 
Exponenten klein gegen 1 ist, wenn also die Bedingung ; 
__(kTyP 
(34) Qmwho, <l 
erfüllt ist, kann man die Reihe nach dem zweiten Glied abbrechen 
und man bekommt aus (24) für den Widerstand in diesem Grenzfall: 


wobei 2K der Grenzwert von K ist 
1 


am) 2 


(35) 
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Um den Ausdruck (34) diskutieren zu können, muß man die 


Integrale abschätzen. Dazu bemerken wir zuerst, daß die Argumente 
der Fermifunktionen von K sich folgendermaßen umformen lassen 


+5 


so daß beide keine negativen Werte annehmen können. Die Fermi- 
funktion ist aber nur für kleine positive Werte des Arguments 
wesentlich von Null verschieden. Da das zweite Argument nur für 
f=0, also £= 0, verschwindet, kann man die entsprechende 
Funktion nach Potenzen des Summanden £ dieses Arguments ent- 
wickeln und man bekommt so für X 


A 2mw 

Wir führen jetzt wieder, wie im vorigen Paragraphen, f und 
und ¢t = - als Integrationsvariable ein. Für das erste Integral hat 


man zuerst folgende Integration über t von 1 bis + oo durchzu- 
führen: 


2kT 2 mw: 


und dann mit (ear) EN (N + 1) zu multiplizieren und über & zu 


integrieren. Für die folgende Rechnung werden wir annehmen, daß 
kT Sm w? ist, was keine wesentliche Einschränkung für die Tempe- 
ratur bedeutet, da mw?/k eine Temperatur von der Größenordnung 1° 
ist. Man kann jetzt das Integral (36), als Integral über die Ab- 
leitung der Fermifunktion, nach der gewöhnlichen Methode berech- 
nen: man entwickelt die Faktoren an der Stelle, wo das Argument 
der Ableitung der Fermifunktion verschwindet und man integriert 
dann die Reihe gliedweise; dabei kann man sich auf das erste Glied 


beschränken. Wenn man noch u = - (1 _ N. 5) als Integra- 
2 2muw? 


tionsvariable se so bekommt man näherungsweise für (36): 
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die wobei 7 = 0,708 ist. Die Integration über &£ ergibt dann für das 
erste Integral der Klammer von (35): 
sen ‘ & 5 
é 
Nach einer ähnlichen FOREN: bekommt man für das zweite Integral 
nil Wenn man diese Keen in (35) einsetzt, so bekommt man “ 
nts für den Widerstand 
ij m? 1 Ez’ a 
nt- und für die relative Widerstandsänderung 
I. ho, I 
§ 9. Der Effekt im longitudinalen Feld 
nd Wenn das elektrische und das magnetische Feld parallel sind, —__ 
al so ist die Berechnung des Widerstandes ganz ähnlich derjenigen, 
die Bloch für den feldfreien Fall durchgeführt hat. Die Bewegung 
ii des Elektrons parallel zu den beiden Feldern ist eine beschleunigte _ 
Bewegung. Es gibt also, im Gegensatz zum Fall eines transversalen 
Magnetfeldes, keine stationären Zustände in beiden Feldern, und man 
muß das elektrische Feld als Ursache einer zeitlichen Änderung der _ 
Verteilungsfunktion betrachten. In der üblichen Weise hat man ~~ 
zu dann diejenige Verteilungsfunktion zu ermitteln, die stationär bleibt 
B unter dem Einfluß des elektrischen Feldes und der Wärmebewegung _ 
ei des Gitters. Für diese stationäre Verteilungsfunktion machen wir _ 
0 den Ansatz x, + w,, wobei z, die Fermifunktion und yw, eine 
% Korrektur ist, die in erster Näherung proportional der elektrischen 
€ Feldstärke F sein soll. Bei der Aufstellung der Stationaritäts-- 
nt bedingung berücksichtigt man nur die linearen Glieder in F,d.h. © 
rt quadratische Glieder in oder das Produkt F werden vernach- 
d lissigt. Die gesamte zeitliche Änderung der Verteilungsfunktion, — 
die im stationären Fall verschwinden muß, setzt sich aus zwei Gliedern _ 
r zusammen. Das eine enthält den Einfluß des elektrischen Feldes und 


oder, nach der oben gemachten 


ist bekanntlich gleich — 7 it 


Voraussetzung, gleich er 


Elektrons fir die zum — bedeutet. 
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wobei / wieder die Wellenzahl des 
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sind. Mit ihrer Hilfe eliminieren wir f, und 2,’ und ersetzen dann 
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Das andere Glied enthält die Änderung, die von der thermischen 
Bewegung des Gitters verursacht wird. Um dieses Glied zu berech- 
nen, braucht man noch die Ubergangswahrscheinlichkeiten. Man 
bekommt sie aber ganz leicht, wenn man in unseren Formeln (8a 


und b) F =0 setzt. Da F nur in den Faktoren u . 
hat man bloß in dem 2 für die Energiedifferenzen ¢,,— &, zu setzen. 
Die Änderung der Verteilungsfunktion, die von den Stößen herrührt, 


wird dann: 


vorkommt, 


— y,)(N +1] 


1 — cos + o)t 


[77 


1-x)N - 21-2) N+D)+W 


(1 +2, N + 1) — (1 


1 — cos +o)t 


2 
o) 


+ w, +2%,N) %(N +1) + (1 N)] 


1 — cos Er = o)t 
wenn man die Glieder mit y, , vernachlässigt. Man kann jetzt 


diesen Ausdruck vereinfachen, wenn man von den Auswahlregeln 
[(10a) und (10b)] Gebrauch macht, die auch in diesem Fall gültig 


+ N 
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die Summationen über die anderen Wellenzahlen durch Integra- woe 
tionen. Es wird: hela 


wobei A eine Abkiirzung fiir die geschweifte Klammer von (38) be- ur 
deutet. 
Es ist zu bemerken, daß y, nicht von der Quantenzahl x, ab- a 
hängt. Für w, werden wir eine Lösung suchen, die dieselbe Eigen- a 
schaft hat. Die einzige Funktion der rechten Seite von (39), die er | 
noch von x, abhängt, ist W,,, da sie aber nur von der Differenz 
ty — z, abhängt, die wir schon durch f, ausgedrückt haben, so ver- 
schwindet die Abhängigkeit von x, vollständig. Ferner enthält der — 
Integrand von (39) f, und f, nur in der Kombination f,? + 6; so 
daß man eine Winkelintegration leicht durchführen kann: ie 


n’ —oo 


Wir werden jetzt in der üblichen Weise den Energiesatz aus- 
nützen, umf’ zu eliminieren. Dafür zerlegen wir das Integral von = im : 
(40) in zwei Integrale, von denen das eine sich auf die Emissions- 
prozesse und das andere auf die Absorptionsprozesse bezieht. Dann _ 
benutzt man den ö-Funktionscharakter der oszillierenden Faktoren 
1 — cos Lt 

$2? 
bekanntlich verschwinden dann die Glieder nullter Ordnung in F. 
Nach erfolgter Differentiation nach t bekommt man 


1 


wobei 


vom Typus , um die Integration über f’- durchzuführen; 


| 4=y, z)N+D+z, vw, [Lor ( 


(41) Er 


sungen der Gleic & + = 0, 
—ho=0,V—1 ay = 0 in die Klammern einzusetzen sind. 


Wir können jetzt die ge, aufstellen: 


eF ih 1 


(f? 4), = 0. 
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Es ist zu vermuten, daß die gestörte stationäre Verteilungsfunktion 
mittels einer kleinen Translation in / aus der ungestörten entsteht: 
(43) + win) = z(n,l— 


also 
y, =z (nl — (1 —z,): 
Daher ist es zweckmäßig, für w, folgenden Ansatz zu machen: 
wobei wir jetzt C, allgemein als Funktion der Quantenzahlen be- 


trachten. Für C, bekommt man aus (42) die Integrodifferenzen- 
gleichung 


| + (f? Waw N — xo) Jui C, »—1C,)=0. 

Für tiefere Temperaturen (T<®) kann man diejenige Methode 
zur Lösung dieser Gleichung anwenden, die Bloch!) für den feld- 
freien Fall angewandt hat, und schließen, daß in erster Näherung 
C = const der Gleichung genügt. Die Bedingung, die wesentlich 
für die Anwendbarkeit der Blochschen Methode war, dß’—I=+f, 
für tiefe Temperaturen sehr klein ist, ist auch in unserem Fall 
erfüllt. Für hohe Temperaturen ist aber unsere Gleichung viel 
schwieriger zu lösen als die Blochsche, da die Abhängigkeit 
von n und n’ sehr unübersichtlich ist. Daher werden wir uns hier 
auf den Fall tiefer Temperaturen beschränken. 

Der Wert der Konstanten, die wir jetzt wieder / nennen 
wollen, ist leicht aus der Gl. (44) zu ermitteln: man multipliziere 
die Gleichung mit /, summiere über n von 0 bis +00 und inte- 
griere über 1 von —oo bis +00. Man bekommt: ea 


> (l—z,) 
n 


II? Wan (1 — +[f°? Wan (1 — X) N),} = 0. 
Da die geschweifte Klammer unter dem Integralzeichen des zweiten 
Gliedes symmetrisch ist in den gestrichenen und ungestrichenen 
Variablen, so ergibt die Summation und Integration des Summanden 


(44) 


mit dem antisymmetrischen Faktor A : ('—1) Null und es bleibt: 


1) F. Bloch, Ztschr. f. Phys. 57. 8. 545. 1929. 
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tion 


(L — N), + Wan (1 — zo) 
Man kann die linke Seite durch eine partielle Integration umformen 


+00 

mk T kT 
er > "7er = ,eFr. 
Bei der letzten Umformung haben wir die Formel benutzt, die die 
Anzahl Z,dl von Zuständen, die eine bestimmte Quantenzahl n 


haben und deren Quantenzahl / in einem Intervall von der Länge di _ Br 
liegt, gibt: u 


Z dl = — mal). 


list aber, wie aus (43) sofort ersichtlich, die mittlere Wellenzahl ae 
in der Richtung der beiden Felder. Fiir die Stromdichte eh 2 a 


man sofort ss = — ev- = und für den W iderstand 


l 


Mit dem aus der Gl. (45) kauen Wert von / a ? 
kommt man 


D? mM 2 


Führt man jetzt f, und f’ statt 1 und I’ als Integrationsvariable ein = 
und benutzt man noch die Identitäten ag 


[Xo (1 —Z¢ = 1) — Za) 
[Xe (1 | (N + 1) (Ka — Xa) Ju; 
so bekommt man 


D? 2 


n,n 
1 
x-= X += 
e 741 


wobei X und X’ die Bedeutung von (20) haben. Die Integration ; 
ist nur über die positiven Werte von f’ und f, durchzuführen. 


1) L. Landau, a.a. 0. 
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Zur weiteren Umformung dieses Ausdrucks (46) können wir 
genau dieselben Methoden anwenden, die wir auf den Ausdruck (24) 
des Widerstandes in transversalem Feld angewendet haben. 

Selbstverständlich liefert die Integration über n und n’ unter 
Benutzung der Ausdrücke (12) für die W,,, genau dasselbe wie im 
transversalen Fall. Das kann man leicht auch rechnerisch be- 
stätigen. Der einzige Unterschied tritt bei der Integration über den 

Azimut in der f’, f,-Ebene auf: im jetzigen Fall hat man das Integra 
> aj2 a/2 
ff cos? sin td statt [ sin? td 
0 0 

als Faktor, und das kompensiert genau den Unterschied zwischen 
den Zahlenfaktoren von (46) und (24). 

Zum Zweck der Berechnung der Korrektur fiir kleine Magnet- 
felder- können wir unser Resultat (30’) unverändert übernehmen. 


Jetzt haben wir aber diesen Ausdruck mit + + zu multiplizieren 


und von 1 bis oo über ¢ zu integrieren. Nach einer Rechnung, die 
ganz ähnlich verläuft wie diejenige, die wir für den Effekt im trans- 
versalen Feld gemacht haben, bekommen wir 


Im Grenzfall starker Magnetfelder kann man auch dieselbe 
Methode anwenden wie im vorigen Paragraphen. Das Ergebnis ist 


4 ho 2 In 
BE th 82 2 mu 


§ 10. Diskussion der Ergebnisse 


Unsere Ergebnisse sind in den Formeln 
O/T 


m? 1 , | 
(k T) J dé limesho,—>0, 


[emus für ho, <KT, 
4-1 für ha, >E,, 


2 
- für ho,<kT, 


kT (I, kT ~ 2 mw 1, 


1 für ho, > E, 
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Im Grenzfall (29) muß der Widerstand identisch sein mit dem- A 
jenigen, den man bekommen würde, wenn man das Problem fürH=0 ° 


nach der Blochschen Methode behandelt. Unser Resultat geht 
aber über das Blochsche hinaus, da unsere Formel (24) für den | 
Widerstand auch für T ~ 9 Gültigkeit beansprucht, so daß auch (29) _ 


in diesem Temperaturbereich als gültig anzusehen ist. Griineisen!) _ 


hat gezeigt, daß die Formel (29) das empirische Material auch in 
dem Gebiet T ~ © gut darstellt. 

Für kleine Magnetfelder liefert unsere Rechnung, im longitudi- 
nalen wie im transversalen Fall, eine quadratische Abhängigkeit 
der Widerstandsänderung vom Magnetfeld. Es ist aber schwer, die 
theoretischen Resultate mit dem empirischen Material zu vergleichen, 
da der gemessene Koeffizient des quadratischen Gesetzes sehr stark 
von der Reinheit und der Bearbeitung des untersuchten Materials 
abhängt; nun führt jede Störung zu einem Zusatzwiderstand, der, 


aus Symmetriegründen, quadratisch vom Magnetfeld abhängt und oe 


man kann die verschiedenen Effekte nicht trennen. 


Für sehr starke Magnetfelder ergibt unsere Rechnung, im 
transversalen wie im longitudinalen Fall, eine Widerstandsänderung, pra 


die linear vom Magnetfeld abhängt. Die Neigung der Geraden, die 


dieses Gesetz in einem Diagramm darstellt, in dem 4o/o als a 


Ordinate und H als Abszisse aufgetragen wird, ist nach Kapitza 
praktisch unabhängig von der Reinheit und der Bearbeitung des 
Materials. Nach unserer Formel (37) hängt die Neigung von der 
scheinbaren Masse des Elektrons ab. Durch den Vergleich mit den 
experimentellen Ergebnissen kann man daher den Wert der schein- 
baren Masse bestimmen. Folgende Tabelle enthält das Ver- 


hältnis m/u der scheinbaren zur gewöhnlichen Masse u des Elektrons 


für einige Metalle: 


Metall | Cu | Ag | Mg | Cd | A Sb | Bi 
mu | 07 |. 10 | 01 | 04 | 05 |11-.10- | 3,3-10-¢ 
Für die Metalle der ersten, zweiten und dritten Spalte des 


periodischen Systems bekommt man also durchaus normale Werte. yey 7 
Im Fall des Wismuts bekommt man aber, als Folge der starken 


Feldabhängigkeit der Widerstandsänderung, einen sehr kleinen Wert 
für die scheinbare Masse. Zu einem ähnlichen Ergebnis ist auch 
Peierls?) gekommen durch den Vergleich zwischen seiner Theorie 


des Diamagnetismus bei tiefen Temperaturen und starken Magnet- __ 


1) E. Griineisen, Leipziger Vorträge 1930. 
2) R. Peierl . 192. 1933. 
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feldern und den Messungen von de Haas und van Alphen; er 
bekommt für das scheinbare Magneton der Bahnbewegung der 
Wismutelektronen einen Wert von 0,5-10-2%, der zu einem Ver- 


hältnis = = 1,8-10~* führt; die Übereinstimmung mit dem aus 


unserer Rechnung folgenden Wert ist also nur qualitativ. 

Über die Temperaturabhängigkeit der Neigung läßt sich wenig 
sagen, da keine systematischen Untersuchungen .darüber gemacht 
worden sind. Kapitza hat seine Messungen nur bei drei Tem- 
peraturen ausgeführt, die nicht sehr weit von der charakteristischen 
Temperatur liegen, für die also unsere Formel keine einfache Tem- 
peraturabhängigkeit darstellt. Es scheint jedoch, daß die Neigung, 
die aus der Theorie folgt, langsamer zunimmt mit sinkender Tem- 
peratur als die gemessene. 

Der Vergleich zwischen Experiment und Theorie für den 
anderen Parameter des linearen Gesetzes, der die Lage des Schnitt- 
punktes der Geraden mit der Abszissenachse bestimmt, ist schwerer 
durchzuführen als im Falle der Neigung, da dieser Parameter, nach 
Kapitza, stark von der Bearbeitung und Reinheit des Metalls ab- 
hängt. Nach unserer Rechnung beginnt der Gültigkeitsbereich des 
linearen Gesetzes erst für Feldstärken, die so groß sind, daß alle 
Elektronen sich in den Zuständen, deren Quantenzahl n gleich Null 
ist, befinden. Die untere Grenze der Feldstärken, die dieser Be- 
dingung genügen, ist gegeben durch 
1 2m a, \*/s 1 2eH \’h 
Es ist zweckmäßig, statt für die Schnittabszissen, den Vergleich für 
diese kritische Feldstärke durchzuführen, um die Anzahl » der 
Elektronen pro Kubikzentimeter zu ermitteln. Wenn man die ge- 
wöhnlichen Annahmen über die Elektronendichte macht, so bekommt 
man für diese kritische Feldstärke einen Wert von der Größen- 
ordnung 10° Gauss. Experimentell ist aber das lineare Gesetz für 
viele Metalle schon bei Feldern von der Größenordnung 10° Gauss 
gültig. Man muß also in diesen Fällen annehmen, daß » viel 
kleiner ist als der Wert, den man bekommt, wenn man die Anzahl 
der Leitungselektronen als gleich der Anzahl der Atome pro Kubik- 
zentimeter betrachtet. Für die Alkalien ist man zu einer solchen 
Annahme nicht gezwungen, da aus den Meßergebnissen von Kapitza, 
auch bei den größten magnetischen Feldstärken, nicht mit Sicherheit 
die Gültigkeit des linearen Gesetzes folgt. 

Für Metalle wie Cadmium und Wismut ist dagegen das lineare 


Gesetz schon bei Feldern, von der Größenordnung 3-10* Gauss 
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sicher gültig; man muß also in diesen Fällen schließen, daß » viel 
kleiner ist als gewöhnlich angenommen. Im Falle des Wismuts ist 
man auch auf anderen Wegen zum selben Schluß gekommen: 
Peierls’) bekommt aus den diamagnetischen Eigenschaften den 
Wertv=2.10'%; Bellia?) bekommt aus der Hallkonstante » = 1,2-10” 
und neuerdings Eucken und Förster?) bekommen aus Unter- 
suchungen über die freie Weglänge der Wismutelektronen Werte 
für » die zwischen 1,63-10!7 und 10,9.10!7 schwanken. Nach 
unserer Formel (49) berechnet man einen Wert » = 5-101". Die 
Übereinstimmung der Ergebnisse, die man nach den verschiedenen 
Methoden bekommt, ist zwar sehr grob, aber man kann doch daraus 
schließen, daß die Anzahl der Leitungselektronen des Wismuts 
wesentlich kleiner ist als die Anzahl der Atome. 


Zum Schluß möchte ich Herrn Prof. Heisenberg für die An- 
regung zu dieser Arbeit und die stete Hilfe bei ihrer Entstehung, 
und Herrn Dr. Bloch für zahlreiche interessante Diskussionen über 
die behandelten Probleme meinen herzlichsten Dank aussprechen. 


1) R. Peierls, a.a.O. 
2) O. Bellia, Ztschr. f. Phys 74.'S. 655. 1932. 
3) A. Eucken u. F. Fö+ ter, Gött. Nachr. Bd. 1. Nr. 3. S. 43. 1934. 


Leipzig, Institut f .eoretische Physik der Universität. 
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Uber den von einer ebenen Erde Teil 
der Strahlung eines vertikalen Dipolsenders 


Von K, F. Niessen 


(Natuurkundig Laboratorium der N. V. Philips’ Gloeilampenfabrieken 
Eindhoven/ Holland) 


Zusammenfassung 


Von der Gesamtenergie, welche von einer vertikalen Dipol- 
antenne emittiert wird (in einer Höhe h oberhalb der ebenen Erd- 
oberfläche aufgestellt) wird ein Teil T in der Erde absorbiert und 
der andere Teil 1 — T durch die (nichtleitende) Luft von der Erde 
weggestrahlt, wenn wir die Ionosphäre vorläufig außer acht lassen. 
Diese Teile T und 1 — T wurden berechnet und eine praktische 
Formel für T abgeleitet [vgl. (38)] für solche Werte von Wellen- 
länge (A) und Bodenkonstanten, bei denen der Verschiebungsstrom 
in bezug auf den Leitungsstrom vernachlässigbar ist und bei denen 
für den Absolutwert (n) des Brechungsindex gilt n >10. 

Zu gleicher Zeit wurde <t angenommen. 
a ity Graphisch findet man die Teile T und 1—T in Fig. 8 ab- 


gebildet als Funktion von n mit «= ni als Parameter. 


Formel (38) und Fig. 8 sind nur zu benutzen in den von den 
Figg. 5, 6 und 7 angegebenen Gebieten für h und A bei ent- 
sprechender Bodenart. 


§ 1. Allgemeine Formeln für die Strahlung nach oben und unten 


Wenn ein elektrischer Dipol entlang der 3-Achse eines ortho- 
gonalen Koordinatensystems x y 3 schwingt, sind die elektrischen und 
magnetischen Kräfte (E und H) bekanntlich von einer Hertzschen 
Funktion /7*(z, y, 3) ableitbar, die der Gleichung 


geniigt, wo 
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Lichtgeschwindigkeit c- (j = v-1) ; 
Im folgenden werden die von Sommerfeld in Riemann- 
Weber II benützten Einheiten angewandt. Man hat dann: 
—je 6 II* +je 6 II* 
dy ’ @ Ox ’ Hi, = 0. 
In Zylinderkoordinaten r, 3, g hat man unter Benutzung eines 
anderen Funktionszeichens 


I (r, d) II* Y; 8) 


IT 
+ Or 


2 6? IT 
E, = k? OF” 


Benutzt man in diesen Koordinaten neben 4 und @ nicht 


und führt als Funktionszeichen ein: 
so ergibt sich hierfür und für die in Fig. 2 angegebenen Kraft- = 
komponenten: | 


[2 
‚ aufgebaut ist aus der Dielektrizitätskonstante e, der elektrischen = ee 
= 4 
il 
ab- 
und für die in Fig. 1 angegebene 
IT 
n je 6 
3 H,= H,= & 
;ho- 
und 7, sondern a 
hen 


Annalen der Physik. 5. Folge. Band 22. 1935 


IT, 23 0°], 


9 


R? dR 
| on, 
Ee sın cos — sin 06 —- 


+ = =0, 


© ist hierbei der Winkel zwischen R und 3-Achse. Letztere Koor- 
dinaten werden meist benutzt, um die Gesamtstrahlung pro Sekunde ©) 
zu berechnen, wenn der Dipol sich in einem homogenen Medium be- 
findet, das den ganzen Raum ausfüllt. 

Hierfür gilt, wenn das Moment des Senders 1 ist: 
ei kk 


und die Strahlung berechnet sich nach der allgemeinen Vorschrift: 


= 


S [ 2rR sin ORAOcHoH,, 
0 


wo der Zeitmittelwert 
EoH, =i. (Ho- H , 
wenn {, bedeutet: Resller Teil von“, während @ den komplex 
konjugierten Wert von a angibt. Man bekommt dann mit den in (4) 
und (5) angegebenen Ausdrücken für Hg und H, die Strahlung: 
= R, {= mc? k? ke 


@ 


Ist das Medium nichtleitend (o = 0), so ist k reell: rs 


k=p 
(p positiv reell) und die Strahlung, die wir mit ©,” angeben wollen 
(Index f wegen der freien Aufstellung im Medium ohne Erde und 
Index p wegen des reellen positiven Wertes von k) wird | jetzt: 


(6) 


Auf als Maß wir immer die Strahlung. in kompli- 


rierteren Umständen beziehen. 
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Solch ein komplizierter Fall liegt vor, wenn wir die Strahlung 
eines vertikalen Dipolsenders betrachten, der in einer Höhe h ober- = 
halb einer ebenen Erde steht. 

Die Erdoberfläche sei 3=0 und die Luft („> 0) habe s, und Bi. 
die Erde (3 < 0) « und 2 als Werte für « und o und die zu- 
gehörigen Werte für k? seien k,? und k,?. Gr, 

Es gilt dann, wenn die Hertssche Fuaktion in der unmittel- 
baren Nihe des Senders diejenige eines freien Dipols mit Moment } ist, 
nach der Sommerfeldschen Rechenmethode fiir die Luft: 


oo 


kt ky? + ky? 3? 


wo J,(Ar) die Besselsche Funktion nullter Ordnung ist, während 
die Wurzeln immer mit einem positiven Realteil definiert sind. 
R, ist der Abstand vom Aufpunkt zum Sender und R, derjenige _ 
vom Aufpunkt zum Bildpunkt des Senders (gespiegelt in der Erd- 
oberfläche). 
Wir werden die ersten zwei Terme im rechten Glied in der 
bekannten Weise durch ein Integral mit Besselfunktionen ersetzen 
und in der Luft eine Ebene 3=h-+s oberhalb des Senders — 
(0<e< oo) und eine Ebene 3=h—e(<e<h) unterhalb ds 
Senders betrachten [e ist also immer eine positive Gréfe')]. Wir be- 
kommen dann: 


Jane” Jy an (th+e)y qq 


Jo une hy? (2h—e) 
> Adi 
vollen 
» und 
if kt ky? + ky? VP — 


1) Verwechslung mit den Brechungsindizes e, und s, wird nicht-zu be- 
fürchten sein. 
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Wenn wir nachher in beiden Fällen e->0 gehen a kann 
Gesamtstrahlung entstanden gedacht werden aus: > 


Srotaa = =n+0 + Gyan-o; 


wovon die erstere durch 3 =h-+« nach oben, die letztere durch 
4 =h —e nach unten ausgestrahlt wird («—> 0). In den Zylinder- 
koordinaten haben wir: 


(2) 


(8) f2are( (B,H |, 


worin die beziiglichen Kraftkomponenten aus 
bw 1,3 


zu bestimmen sind. 
Wir berechnen zuerst die nach oben gehende Strahlung ©, -;4,. 
Dafür ist nach (2) und (3) und bei Benutzung von J/, (3 =h + 8): 


aT, ~ (2%) 
a=h+e 


Orde 


r Or 03 


DP 


(Ar)e- 


0 
(ar) VP ke? 


Jy (ir) e” 1? dl 


’ - (2 
Jy’ (Ir) e 18 — (2h + €) dl 


Substitution von (9) und (10) in (7) und Vertauschung der 
Integrationsfolge führt zu: 


» . 
—— 
Wiis 
= — > | d, 
; 
| ky J 

\ 
‘ 


3 
A=0 
42 (2h + e) 


r=» 


Nun ist?): 


- 


0 Fi 
für Besselfunktionen der ne Ordnung und mit ö(x) als Diese: 


funktion von x, wofür bekanntlich 


wihrend 
(13) =0 fir «+0. 
Hieraus folgt also wegen J,’ (2) = J, (2): 


= (ar)J,(Ir)rdr 
=0 0 


1) Vgl. 
S. 560. 1932. 


q 
K.F. Niessen. Uber den von einer ebenen Erde absorbierten Teilusw. 167 
R, [22 ( H,);=n+.dr 
ch 2 Vie —k? 
. 
r 
Auch gilt nach (12) und (13): 
der Dun z.B. Balth. van der Pol u. K.F.Niessen, Phil. Mag. 13. a a 


4=0 


Nach Substitution von (14) in (11) kénnen wir (11) schreiben 


in der Form: 
l=co 


t=0 


F (4) = {-+-} von (11), 
G (l) = [---] von (11). 

Die Besselfunktionen sind bei dieser leteaneiion: der Strahlung 
über die ganze Ebene 3 = +. verschwunden. Ähnliches geschah 
bei L. J. Peters und J. Bardeen?) in einem anderen komplizierten 
Problem. 

Nach Anwendung von (15) kommt: 


l= 


wo | | den Absolutwert me 
Da F(l) auch den Faktor e~ V#-%*+ enthält, können wir auch 
schreiben: 


= 00 


l 
l= 


worin f(l) die Größe & nicht enthält: RE 


1 1 2h -Ve- 


18 l k,? . 
( )FO 2 Vi 2 + VP— + kt VR- k? 


Bevor wir diese Formel für &,-.+. für einen Sonderfall weiter 


vereinfachen, wollen wir zuerst auch die allgemeine Formel für die 
nach unten gerichtete Strahlung &,-„-. angeben. er 4 


1) L. J. Peters u. J. 2.8. 1088. 4 
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= 
In genau derselben Weise wie soeben finden wir dann aus (8), a 
aber jetzt bei Benutzung von J/, =h — 8), 


kB F, OG, (dl, 


i=0 


2 k,? ye- hy? + ky? 

1 
= + = ® 
hey? + Ve- k? 
Diese F,(l) und G,(l) sind nicht konjugiert komplex und sind 

nicht zerlegbar in einen Faktor ohne und einen KR mit «, wie 


2.B. 


§ 2. Strahlung bei nichtleitender Luft = 

Die Formeln (17) und (19) für Ganz. und ©,-ı-: wolle 

wir nun zuerst betrachten für den Sonderfall 
k?= p? (p positiv reell), 

wobei wir also von der Leitfähigkeit (o,) der Luft absehen, 
meistens erlaubt ist. 

Betrachten wir dabei zuerst: ©,-ı.+. 

Da nun jk?=jk?=jp? rein imaginär ist und wir in (17) 
nur den Realteil benötigen, braucht die Integration nur von / = 0 
bei l=p stattzufinden. 

Zu gleicher Zeit sehen wir, daß dann die s-Abhängigkeit, die 
ja nur im Faktor 


auftrat, gänzlich verschwindet, denn fiir — 
k’= pP, 
0 < 12 < p? 


1, 

_ Dieses Ergebnis war zu erwarten, denn bei k,? = reeil (nicht- 
leitende Luft) wird keine Energie von der Luftschicht zwischen den 
Flächen z=h+e® und 3=h-+" absorbiert, so daß alles, was 
durch die Ebene 4=h+<' nach oben geht, auch die Ebene 
§=h-+ :” in derselben Richtung passiert. 

Die der Annahme k,?= p? entsprechende Formel für &,-.-. 


(die wir andeuten mit ©’_,„_,) kann nur angegeben werden, wenn 


etzt mit: 
a 
{ 
| 
a 
I; 
2 
) daß - & 
= 
MY. 
2x 
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zuerst V1? — k,?= Yl?—p? in der richtigen Weise bestimmt wird. 
Dazu bedenken wir, daß k,? nach (1) eigentlich dem ersten Qua- 
dranten der komplexen Ebene 
zugehört (für einen Augenblick 
also 6, > 0 vorausgesetzt), so 
daß k,? im allgemeinen im 
dritten Quadranten liegt und des- 
halb YI?— k,? im vierten. Des- 
halb muß in Fig. 3 für das in 
Frage kommende Intervall 


genommen werden: 
Es ergibt sich dann nach (17) und (18): 
l=p 
l=0 
1 1 1 

18a) f ()=- - > — 
( fl) 2 2 Vp?— ky?— j ky*Vp?-P 
wie gesagt ohne irgendeine Abbangigkeit von & 

Aus demselben physikalischen Grunde wie bei ©? muß 


, von s unabhängig 


auch die nach unten gehende Strahlung ©... 
werden, sobald die Annahme k,*= p? gemacht wird, was wir jetzt 
zur Kontrolle an der Gl. (19) untersuchen wollen. Während die 
Annahme k,?= p? in (17) sofort die Integration von 0 bis oo auf 
eine von 0 bis p reduzierte, geschieht Ähnliches nicht in (19), denn 
wir bemerkten schon früher, daß F,(l)x G,(l) nicht reell war [im 
Gegensatz zu F(l) x G(l) = FO g Wi ir bekommen denn auch jetzt 
sowohl Integrale von 1 = 0 bis | = p wie einen von | = p bis l= 00, 
nämlich: 


l=p 
p? Ik _ FP dl 
t=0 


iJ 3 


= 0o 


L=p 


3 bedeutet N Teil von“, 


3 
‘ 
aX. 
SEA 
| 
= 
| 
4 
ve 
atte 
— 
> 
— 
. 
X 
— 
te 
| 


Wie man sieht, ist auch hier s verschwunden, genau so wie bei _ 
bei k’= p°. 


$3. Strahlung bei nichtleitender Luft und gutleitender Erde 
(0, > 8, 0) : 
Wir wollen nun unsere Einschränkung noch enger machen und — 
beide Strahlungen fiir den Sonderfall 
k?=p? und k?=jP:? 
(p und P positiv reell) unter- 
suchen, d. h. wir nehmen an: 
©_,,. und ©_,_, können 
dann noch mit nun bekanntem = 
Wert für yl®—k,? versehen 
werden, was die Ausrechnung 
der benötigten Moduli und Ima- 
ginärteile ermöglicht. 
3 Wir haben in Fig. 4 
VP—k,? = Al) —j BY) 
mit den Abkiirzungen: 


ee A(l) = OW cos o = y% + P* cos g, 
BO =O0Wsing = VF + Pising. 


Wir finden dann aus (17a) für die Strahlung durch 4 =h +e x 
nach oben im Falle k,*= p? und k,?=j P? (was mit Index p, P ar 
angegeben wird): 


l=p 
2m? 
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rd. 
= 
ne 
ick 
so a 
eS- 
| « 
eS- 
in 
12, 
3 
3 
| 
. 
18 
np As = an 
tat cos 2 = ——— 
4,17 
Ll, 
3 
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l 
+ 
Die nach unten gerichtete Strahlung , (also auch für 
k,?= p? und k,?= dagegen wird nach 


», _ 
oP 


l=p 
PB dl 


‚Vee- + Vit+ 


l=p 
Dieser Ausdruck, zusammen mit Cr nee LYgl (20)], geben in 


ar Limite « — 0 die Gesamtstrahlung des Senders im Falle 
= p*, k,?= 7 


§ 4. Strahlungsverteilung in Sommerfelds Spezialfall 


Die Ausdriicke (20) und (21) werden wir noch vereinfachen ; 
tels der weiteren Annahme }): 
Für den Absolutwert n des Brechungsindizes hat man dann: 
4 
—=n>1. 
> 
1) Bei diesen Annahmen bezüglich Bodenart und Wellenlänge leitete 
Sommerfeld bekanntlich die Hertzsche IZ-Funktion ab als Funktion einer 
von ihm definierten numerischen Distanz, und zwar im Falle eines auf der 


Erde stehenden Dipols (h = 0). 


fey 
BEN. 
Are 
3 
Spa 
Pr 
- 
Phe 
= - = 
v4 
(21) nc? p* 
© + 2 | 
l= 
| 
a 
eat) 
4 
* 
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Die Annahme PS > erlaubt in den Integralen von bu 
die Annäherung: it 
A —_—; B l 
BO) 


wonach sich ergibt: 


af2 

cos’ ycopt+d 

0 


1 
mit den Abkürzungen: 


= 


nc?p 
@ 


Der zweite Term im rechten Glied kann noch zerlegt werden ’ 
in einen Beitrag, der nicht, und in einen anderen, der wohl für 


P= (wo gutleitende Erde!) verschwindet. 
Man bekommt schlieBlich: 


@ j=hte 
mit der dimensionalen Zahl: 
af2 
yn>1 =f sin’ gdp 
s=hte J 


4 


TB (1 — cos 2h p 3 - sin 2h p34) +ö(l — cos 2hp 4) 


Behandeln wir 


n>1 


in ähnlicher Weise, so kommt aus (21): 


(23) 
x 


=h-e co?p+ycosp + 0 2 
1 


Die Gesamtstrahlung des Senders ergibt sich dann noch: 


Som = nc?p P* 
h-0 


0 bis 
- (20) se wir dann / = p sin g und vereinigen ‘erme ts a 
0 t 
IT 
| 
yr =_ 
4 
| 
iir 
i 
i 
| 
n (22) 
+ 
« 
e 
| 
r 
4 > 
XUN 
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(Da beide Bestandteile, dank k,?= p? von ¢ unabhängig waren, war es 
hier der Limitübergang « — 0 eigentlich überflüssig.) en 
Sei T der Bruchteil der Gesamtstrahlung, der in die Erde 
hineingeht und dort absorbiert wird und also 1 — T der Bruchteil 8 en 
der Gesamtenergie der nach oben ven der Erde hinweggeht, so be- 
kommen wir in unserem Fall mit den bis jetzt gemachten Be- 


gt 

schränkungen (k,?= p?, k,?= 7 P?, PS p): 
Rechtes Glied von (23) as Y 
= R. Gl. von (22) + R.Gl. von (23) p 


85 Über die ganze ausgestrahlte Energie bei den Zenneck-Sommer- 
feldschen Annahmen und beschränkter Höhe des Dipols 
Die ausgestrahlte Gesamtenergie interessiert uns nicht nur in 
bezug auf den Nenner im Ausdruck für 7, sondern hat auch an (2 
sich Interesse, da auch hierin die Anwesenheit der Erde sich be- 
merken läßt. 
Leicht zeigt sich bei Ausrechnung, daß in der Summe von (22) 


und (23) diejenigen Terme sich gegenseitig aufheben, die eine der ” 
kleinen Größen y oder ö als Faktor enthalten und überdies noch 
von h unabhängig sind. Man bedenke dabei, daß in (22) solch ein 
kleiner Term entsteht nach der Zerlegung: 
cos? sin® gdp yeosm +5 
cos’ + +d 3 ( 
sin® » (5 cos + 6) 
\2 do + | - 
cos? + ycosy + 0 P J +9) 
0 0 
l 
8 
1 
: 
t 


’ Die ersten zwei Terme rühren her von der Istrien vom | 
Bender und dessen PER in bezug auf die re denn ( 


4% 


4 
: > 
7 
a 
wen 
= 
e 
| J 
| 
4 
¥ 
ag 
= 
en 
: 
3 
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es sind nur diese Terme, welche übrig bleiben im Falle einer un- 
endlich leitenden Erde (|k,?| = 00, y=d= 0). 
Die letzten zwei Terme berücksichtigen die endliche (statt un- 
endliche) Leitfähigkeit der Erde. 
Der zweite Term berechnet sich sofort durch partielle Inte- 
p 
PY2 


5 transformiert, wodurch unter Benutzung von 


gration, der dritte wird mit 3 = — u und der vierte mit 


Pp 4 nc? p® . 1 
z<!i für die Gesamtstrahlung in —-— als Einheit: 


nV2 
yp (u + 2) cos xu + usin xu 
2u +2 du 


total 


(Absolutwert des Brechungsindizes) 


P 
p 
1 sin « cos «& r 


a sin (@ — x) — cos (a — x) 


cos (a — x) + sin (a — x) — 1 


x= 


Nehmen wir hierin zur Kontrolle den Sonderfall einer unend- 
lich gutleitenden Erde, also P = 00, wodurch n = 00 und x = 0, 
so kommt: 


= 00 


(26) s* = 00 


total, h=O 


Dieser Wert war zu erwarten, denn bei der Berechnung der 
Hertzschen Funktion nach dem Sommerfeldschen Ansatz war 
ein vertikaler Dipol vom Moment } betrachtet worden, d.h. einer 


var 
| 
aa “Pa 
7 
‘ 
1: 
an (24) 
4 
_ 
3 
hes 
KEITH 
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wofür die Hertzsche Funktion in der unmittelbaren — 
des Senders den Wert hat: 


1 

Ein freier Dipol vom Moment } würde, da die Gemein 
lung mit dem Quadrat des Momentes proportional geht, nach (6) 


ausstrahlen: 


3 

Ist dieser Dipol vom Moment 4 nicht frei, sondern steht er 

auf der Grenzebene zweier Halbräume, nämlich Luft und unendlich 
gutleitende Erde, so wird bekanntlich die Hertzsche Funktion in 
der Luft verdoppelt (das Moment auf 2 x 4 gebracht), während nach 
unten nichts ausgestrahlt wird, so daß S die Hälfte des S-Wertes 


eines freien Dipols mit Moment 1 wird, also 


- (also S= 3)‘ 


was gerade mit (26) übereinstimmt. 
Zur weiteren Durchführung der Berechnung der Gesamtenergie 
wollen wir die Reihe der sukzessiven Beschränkungen: 


noch : i einer vierten Annahme vermehren: 


| 
was nach (25) eigentlich er auf 4 
n | 


was öfters der Fall sein wird, da n sehr groß sein kann. : 
Im folgenden müssen wir immer h > 0 voraussetzen, was 
übrigens dem physikalischen Sinn von h als Höhe des Schwer- 
punktes einer wirklichen Antenne entspricht. (Da P>p, also 
eine gutleitende Erde vorausgesetzt ist, würde h = 0 sozusagen 
Kurzschluß geben, was mathematisch zu divergenten Integralen für 
die Strahlung in die Erde führen muß, wenn trotzdem noch eine 
Strahlung e/4 %/R, in die Luft verlangt wird.) 
a oo Der Annahme x <1 zufolge dürfen wir nun in (24) schreiben: 


: f *“udu -/[; e~ **udu fF 
w— 2u+ 2 —2u+2 
0 0 
Der erste Term im rechten Glied gibt bei Entwicklung der 


e-Potenz Terme vom Typus 
1/x 


rn 


q 
4 
Ber»! 
Umgebun d 
A 
A 
rt 
Ws 
Ag 
ag 
Bere 
| Az 
} ‘4 
IE 
: 
5 
TB 
> 
a 
# 
r 


was 
wer- 
also 


gen 
fiir 
eine 


ben: 


der 
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<1 wd 


deren hinzugefügter Näherungswert auf die Annahme x < 
m = 1,2, 3 basiert ist. ; 
Man bekommt also annäherungsweise 


1/x 
e *“udu =f udu 
u®—2u+2 uh — 2u + 2 
wo S (x) definiert ist durch 
1 
und folgendermaßen transformiert werden kann: 
xz S'(r)=1—e-?, 


ER 


Weiter ergibt sich nach PR mit nachheriger Berücksich- © 
tigung von x<l 
1/x 


u" —2u+2 


so daB schlieBlich: 


e~*“udu 


Nun ist bekanntlich !) 
1 1 
f 1 e ® ” =/7 Konstante) = 


(28) K= 


Weiter läßt sich auch das in ie gleichfalls benötigte Integral 


nV 
u +2u +2 


bei der Annahme x <1 berechnen. eta 
1) Whittaker and ee Modern Analysis 8. 243. am. 
Annalen der Physik. 5. Folge. 22. ; 


7 
Zz, = 
1 (6) 4 
- 
n 4 | 
artes 
2 4 4 
x 
R 
BER 
| ® 
ane 
..o ig 
* 
4 
d 
|| 
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Durch Entwicklung von cosxu und sin x u und name um 


u — bekommen wir Terme vom Typus: 


alx 

m m+1 m 

Ww. du m — a 

mid u +2u +2 m! m 

und andere vom Typus 
yp 


ajn AEN 
u™ d x™ 1 a\mn-1 
ui + 2u+2 

_ welche letztere fortgelassen werden können, während erstere den 
hinzugefiigten Näherungswert besitzen wegen also 


In solcher Weise ergibt sich 
nV: 
+ 2)cosxu + usin «u fw 
u +2u+2 + 2u +2 


0 
F (a) definiert ist mit 


=" 


13 


und folgendermaßen transformiert: werden kann: 


a F’ (a) = cosa + sina —1, 


cosa + sina — 1 ie 


~ Das erste Integral im rechten Glied von (29) ist sofort berechenbar 
und dessen Wert wird bei Berücksichtigung von « (wegen n> 1) 


u+2 

J ue + 2u+2 du 

so daß das gesuchte Integral wird: 

(u + 2) coszu + usin xu a 

+ 2u+2 


(32) 


0 
Die Differenz der Integrale (32) und (27), wie sie in (24) vorkommen, 
enthält also keinen Term lg mehr. 

Bei Substitution in (24) kommt (während x in C [a, x] [vgl. 25] 
 fortgelassen werden kann): 


n>1 


| 
| 118 K. 
. 
fl 
: 
Zu 
j 
+ F(a), 
\ 
r 
= 5 
at 
f 
> 
- 
3 
Br: 
| 
1: 
ig 
- 
3 
1: = 
2 
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8 mit: 
fi 
1 
1 
1 
1 — cosa — coB — 
= Sie) + Cie) — | da. 
0 
Nun ist?) 
1 
1 
1 — cosa — cos — ; 
0 
1) Denn: 
1 
cosa + 
2 —de= da + 
a 
) 
dt + 
t 
0 
1 1 
# t 
—dt+ ad dt + 
| t t 
Bei der Substitution ¢ = -- in den letzten zwei Integralen heben diese die oh oe My 


ersten zwei Integrale auf und es wird L = 0. 


. 
j : 
> 
a XUN 
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Wir bekommen also: 


mit: 
1 + (a — 1) cose — (a 1) sine 


(34) g(a) = Si(w) + Ci(a) + | 
Substitution in (83) gibt: 


awa. 4 “r 
und damit ist die sämtliche pro Sekunde ausgestrahlte Energie 
funden: 
5a) 
Hierbei wurde nichtleitende Luft und gutleitende Erde voraus- 
gesetzt, mathematisch festgelegt in k,*? = p?, k,?=jP?, P>p, 
P und p positiv. Weiter warn = -, “a= ., während D(e) und 
g («) in (25) und (34) definiert wurden. 
Wir erinnern daran, daß diese Berechnungen basiert waren 
auf einen Dipol, dessen Moment durch untenstehendes Verhalten 
der Hertzschen JJ-Funktion in der nächsten Nähe ua des | 
Senders festgelegt war: 


Als Energiemaß können wir auch die Gesamtenergie ©, wählen, 


welche von einem freien Dipol mit Moment 1 in einem homogenen 


Medium (Luft é=1, p=k,= = = En nach (6) pro Sekunde aus- 


gestrahlt wird: 
= $ 0 
Einführung dieser Größe in 
3 
—— 
4n V2 


mit den benötigten Größen: 


> 
(35 b) = 


sin « = 
+ ) V2. 


= Si( a) + Ci(a) + 1 + (a — 1) cosa — (a2 + l)sine 


= Absolutwert des Brechungsindiz 


ine 
5. Folge, Band 22. 195 | 
| 
Er, 
i 
a 
- 
. 
Pw. 
fy 
i 
1 
— | 
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§ 6. Berechnung des absorbierten Teiles 7’ der Gesamtenergie bei a 
Zenneck-Sommerfeldschen Annahmen über Bodenart und Wellen- — aa) 
länge und bei beschränkter Höhe des Dipols 8 


Zur Berechnung yon T brauchen wir neben der soeben be- 
stimmten Gesamtenergie noch die in der Erde absorbierte Energie 


@ g=h—e 


Im Ausdruck (23) fiir Kr ae setzen wir im ersten Term ° 


COs = (-—1-+1) 


= 


und im zweiten 


p 
Py2 


Das benötigte Integral wurde schon durch (27) und (28) bestimmt, 
so daß wir bei Substitution finden: 


n>1 


(36) [2lgn —lga +H} 


mit 


y = 0,494 . 


Der Bruchteil T, der die nach der Erde gehende 
in bezug auf die Gesamtstrahlung angibt, wird also nach (35) und (36) 
(37) T = 2 ign — Iga + 0,494 


n D(a) — g(a) + 7 


Die hierin auftretenden Größen D(a) und g(a) findet man in (25) 4 
und (34), während a = en. , n= tal. Alles war basiert auf die 

1 
Annahmen k,?= p*, k,’=jP?, P>p. P und p positiv reell, 


n 


h 


Im Grenzfall » =oo wird T = a —> 0, was der Fall 


sein sollte, da in eine unendlich gutleitende Erde kein Feld ein- 
ingen | kann, 


3 | 
| 
us- 
Bor: 
2 
he 
ren 
es 
| 
H=—, 
} 
XUM 


Annalen der Physik. 5. Folge. Band 22. 1935 


Für kleine @ bekommt man unter Anwendung der bekannten 
Entwicklungen für Si(«) und Ci(a): 


Die) — y2 (5 


g(a) —> (r- 3) +5 


Der Fall « = 0 würde also zu T = 1 führen in Übereinstimmung 
mit der in § 5 (S. 143) gemachten Bemerkung. Wie wir später 
sehen werden, brauchen wir diesen Grenzfall nicht zu betrachten. 

Führen wir Briggsche Logarithmen ein, die wir mit Log. 
(statt lg für die Neperianschen) andeuten, so wird: 


_ Log. n + Ua) 
(38) + Wü) 
mit 


U («) = 0,107 — + Log.@, 


3 


V (a) = 0,307 (3 + 


W (a) = 0341-0217 [Si(a)+Ci(a) + 


Für eine graphische Abbildung von T als Funktion von n und « 
müssen wir diese Größen in Übereinstimmung wählen mit den bei 
der Ableitung der Formel (38) unterwegs gemachten Annahmen, die 
gerade soeben releviert wurden. 

So brauchen wir wegen der reell, bzw. rein imaginär an- 
genommenen OR von 7 und k,? zu ae] 


wovon die erstere in nichtleitender Luft (o, = 0) sowieso erfüllt i 
Letztere erfordert wegen 


(o in C.G.S.-Einheiten), 
eine Beschränkung auf: 
(39) 


Auch hatten wir noch angenommen: 
LA 
k, 


| 
ai 
2 
ba 
( 
f 
N 
ag 
: 
1 ag rel 
An 
3 
= 
: 
: 
a 
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Wir werden eine Größenordnungsdifferenz durch einen Faktor 10 
ausdrücken, so daß wir den Bedingungen (39) und (40) genügen 
wollen durch die Wahl: 


(2 in Zentimetern), 


2co 


(42) i> oe (A in Zentimetern). 


Betrachten wir nun z.B. die 

für «, und o: 
für Seewasser . . . . = 80, 107" (Fall I), 
für Süßwasser . . . . ,=80, 107 (Fall ID, 
für nassen Grund. . . = 10, o=5-10714 (Fall II), 
für trockenen Grund. . 4, = 10725 (Fall IV). 


Wir bekommen dann: 


nach (39) nach (40) 


für Seewasser . . . . A> 13 m i> om 
(43) E Süßwaser . . . . A>13000m A> 1600 


für nassen Grund . . . A> 330m A> 330 
für trockenen ara . . A> 7000 m 4 > 17000 


werden, nehmen wir überall (d. h. für jede Bodenart) 
(44) 4 < 20000 m. 


Neben den jetzt schon berücksichtigten Annahmen o,=0, 
o,>«é,@ und n>1 wurde noch für die Höhe h des Dipols — 
halb der Erdoberfläche angenommen: 


= (hk und A in Zentimetern). 
Wir wählen deshalb: | : 
(45) 

was also bedeutet: 


fiir Seewasser 
für Süßwasser 
für nassen Grund 


für trockenen Grund . 


1) J. Zenneck, Ann. 


S. 721. 1909. 


2 
iter 3 
ten. 
108. 
| 
| 
a 
bei 
die 
un- 
st, 4 
(45a) 
. 
3 4s, | nun in Metern 
400 3 
1 
Phys. 23. S.859. 1907; A.Sommerfeld, Ann. = 
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wobei die Höhe h und die Wellenlänge A jetzt in Metern aus- 
zudrücken sind. 

Innerhalb der von (43) und (44) angegebenen Bereiche für die 
Wellenlänge darf die Höhe dann die von (45a) angegebene Grenze 
nicht übersteigen, damit Formel (38) anwendbar sei, und zwar in 


10 


wobei !/,, gegen 1 vernachlässigt wurde. 
Fig. 5 wurden diese (4, h) Bereiche für die erwähnten Bodenarten 
gezeichnet. Nach oben werden sie von der gezeichneten Kurve be- 
grenzt, nach unten von der Abszissenachse. 

Für Seewasser und nassen Grund gibt Fig. 6 ein vergrößertes 
Bild bei mittleren Wellenlängen und Fig. 7 ein noch stärker ver- 
größertes speziell für die sehr kurzen (oberhalb Seewasser) zu- 
lässigen Wellen. In der Ableitung von (38) wurde ausdrücklich h > 0 
angenommen, um divergente Integrale 


N» 


Dipolhöhe h in KM 


, 
. 
di 
M 
lie 
> 
Wellenlänge A in KM 2 
ad 
a 
I 
: ae 


aus- 


die 
enze 
r in 


liegen müssen. 


die Antenne ersetzen 


der Erde identifiziert werden. 
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nahme ist sehr naturgemäß, denn h ist die Höhe des Dipols, der 
muß und wird meistenfalls nahezu 
Mitte des vertikalen Antennenteiles, also oberhalb der Erdoberfläche, 


in der 


Bei Flugzeugen kann h mit der Höhe des Flugzeuges oberhalb 


Für eine Halbwellenlängeantenne, 


25 


20 


=~ 
DH 


Dipolhöhe h in KM 


Ss 
Ey, 2 


05 / 


- 
- 

- 


of 


05 


10 


Wellenlänge A in KM 
Fig. 6 


Entfern 


für eine langwellige Antenne mit einem vertikalen Teil von z. B. 
300 m h = 150 m sein, wenn die Stromamplitude als Funktion der 
Höhe nahezu konstant ist, sonst jedenfalls h ~ 100 m, dem Schwer- 
punkt der Stromamplitudenverteilung entsprechend. 
eine punktierte Kurve als Größenordnung der kleinsten h-Werte: 
h=12 (für 2A<600 m) und k = 150 m für A > 600 m. 

Man sieht in Fig. 6, wie dadurch noch der erlaubte erreich- 
bare (A,h)-Bereich beeinträchtigt wird. 
für h und 2 gibt Fig. 8 die fraglichen Teile T (in die Erde gehend) 
und (1—T) [von der Erde hinweggehend'!)] graphisch als Funktion 


ungen ermöglicht. 


deren unterer Teil auf 
einer Höhe v oberhalb 
der Erdoberfläche steht, 
würde h=v+yA sein 
and wenn er unmittel- 


bar auf der Erde steht: 
h=1i. Dagegen würde 
150 
= 
24 
SQ 
50 ye 
0 
ee, 50 100 
as Wellenlänge A inM 


Fig. 7 


In Fig. 6 gibt 


Innerhalb dieser Bereiche 


1) die in Wirklichkeit von der Ionosphäre reflektiert wird, und die Signal- 
übertragung auf größere 


> 
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\ 
PL: 
| 
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4nh 
vonn und mit e[= — als Parameter. Diese T-Kurven konnten 


nur für den Gültigkeitsbereich von (A, h) aus (38) berechnet werden, 
d.h. für A > Amin, welche untere Grenze nach (43) mit den Boden- 
arten sehr variierte. Das bedeutete wegen n = Y2co eine untere 
Grenze nin für n, und da wir nur A< 20 km betrachten wollten, 
ergab sich auch eine obere Grenze Mmax, beide abhängig von der 
Bodenart. Man bekommt so die Bereiche: 
für Seewasser . . . . . . 30 <n< 1100, 
| Süßwasser . ..38<n<B4, 
(#6) , nassen Grund. . . . . 10<n< 80, 
| trockenen Grund . . . 10¢n< ill. 


ruchteil 1- nach oben 
ausgestrahit 


| 


T 


25) Aude 
| Bruchteil T 
inder Erde absorbiert, 
| | 
it 


730 N=40 250 
Fig. 8 


Um nun in einem praktischen Fall das Verhältnis der ab- 
sorbierten und der nach oben gestrahlten Energie, d.h. T:(1—-7) 
zu bestimmen, muß man eine der vier erwähnten Bodenarten an- 
nehmen (oder in geeigneter Weise dazwischen interpolieren), sich 
dann überzeugen, ob man sich in einem erlaubten (A, h)- Bereich 
befindet und dann nach n = Y2coA (A in Zentimetern) den Modulus n 
des komplexen Brechungsindizes berechnen, d.h. 

für Seewasser . . mittels: n = Y60A, 

„ Süßwasser . . n = V0,06 2, | 

| „ nassen Grund . n = V0,34, 
trockenen Grund n = V 0,0064. 


berechnet, so gibt Fig. 


| 2 in Metern. 


Ist n so gefunden und wird«e = un 


sofort die erwünschten Teilbeträge T und 1—T. 
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Die Kurven werden gezeichnet nach Berechnung folgender 
Werte: 


a = 0,5 a=] a=2 


0,515 0,512 0,560 

0,445 0,409 0,443 
n = 20 0,349 0,342 0,370 a 
n = 30 0,267 0,261 ee 

n = 40 0,218 0,213 we 
n = 60 T = 0,170 0,162 0,158 

n = 80 T = 0,136 0130 | 017 


| 
uo 


Schließen wir vorläufig Kurzwellensender’) aus (die übrigens 
in unserer Annäherung nur oberhalb Seewasser zulässig sein würden) 
und werden gleichfalls Sender!) in Flugzeugen außer Betracht ge- 
lassen, so können wir nach den Bemerkungen über Dipolhöhe und 
Wellenlänge behaupten, daß für die auf der Erde stehenden Sender 
Anh 


etwa liegt im Bereich von 0000. 
> 


= a=0,1 bis amd 
z. B. z. B. {h = i}. 


die Größe & = 


\li= 20km 4 


Aus den in Fig. 8 gezeichneten Kurven, die T als Funktion 
von n bei « als Parameter angeben (und zwar bei « = 0,25, 0,5, 
1 und 2) geht hervor, daß für diese Kategorie von Sendern die 
Höhe des ersetzenden Dipols ziemlich gleichgültig ist, wenn es sich 
um den von der Erde absorbierten Energieteil handelt. 

Die durch Konstruktionsänderungen praktisch erreichbaren Ab- 
änderungen in der Höhe h des Dipols (etwa durch Änderung der 
Stromamplitudoverteilung über die Antenne oder durch Erhöhung 
der ganzen Antenne unter Beibehaltung ihrer Länge) wird also das 
Verhältnis T:(1— T) zwischen absorbierter und hinaufgestrahlter 
Energie nicht wesentlich ändern können. 

Unsere Behauptung über den geringen Einfluß einer Abänderung 
der Dipolhöhe gilt nur für diejenigen Kombinationen von Wellen- 
länge und Bodenart, die auch in Figg. 5, 6, 7 zugelassen wurden, 
nämlich 

100 m < 4 < 20000 m, 
Süßwasser. . . . . . 13000m<4< 20000 m, 
nassen Grund 

trockenen Grund 


1) Wir erinnern daran, daB alles basiert war auf vertikalen Antennen. 
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Wie es für die anderen A-Bereiche sein würde, kann mit der 
vereinfachten Formel (38) nicht untersucht werden. 

Diese Formel erlaubt aber wohl, die Frage nach der prozentualen 
Absorption zu beantworten für Kurzwellensender oberhalb Seewasser 


allerlei Werte 


bekommen, die das soeben betrachtete Gebiet 0,1 < « < 3 weit über- 
treffen. Natürlich muß man Wellenlänge und Dipolhöhe in Über- 
einstimmung mit den in Fig. 5 angegebenen (A, h)-Bereichen wählen, 
so daß h noch beschränkt bleibt, wiewohl z. B. «= 10 leicht zu er- 
reichen ist. 

Nach der allgemeinen Vorschrift kann man auch hier den 
Bruchteil T der absorbierten Energie bestimmen. Zwar kann man 
nach Berechnung von n mittels (40) oder (47) bei dem vorliegenden 
a-Wert nicht sofort die Bruchteile T und 1—T aus einer T-Kurve 
ablesen, weil diese nur für 0,25 < «<2 bestimmt worden sind, 
T ist aber sofort aus (38) zu berechnen. 


Eindhoven, den 1. Dezember 1934, 


und für Sender in Flugzeugen. Hier kann a@ = 


aa (Eingegangen 12. Dezember 1934) 
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Verte 
7 von Gasmischungen 
über- 
Jber- XXX. Die innere Reibung bei tiefen Temperaturen 
hlen, von Wasserstoff, Helium und Neon 
1 er- und binären Gemischen davon bis 90,0° abs. herab 
Von Max Trautz und Helmut Zimmermann 

den 

man (Mitteilung aus dem Physikalisch-Chemischen Institut der Universität Heidelberg) 
nden Inhalt: 1. Fragestellung, Meßverfahren, Berechnungsweise. — 2. Messungs- 
urve ergebnisse: a) 7 für He und Ne bis 90,0° abs. herab; b) 7 von Gemischen bei 
i 90,0° abs. — 3. Verwertung der Messungsergebnisse. — 4. Zusammenfassung. 
sind, g g 


— 5. Literaturverzeichnis. 


l. Fragestellung, Meßverfahren, Berechnungsweise 


Es fehlten Reibungsmessungen an binären Gemischen von H,, 
He, Ne für das Tieftemperaturgebiet und Neubestimmung von 7 für 
Neon (1, 2). 

Zur Durchströmungsmethode [Trautz und Weizel ( (3)] diente 
als Vergleichsgas bei 0° C Luft, bei 194,6° und 90,0° abs. H,. 


Die Versuchszeit betrug mehr als 200 Sek., damit die Hagenbachkorrektion 
nicht zu groß wurde, höchstens aber 900 Sek., um der Temperaturkonstanz inner- 
halb 0,1° sicher zu sein. 

Die senkrecht stehenden Kapillaren konnten wegen des Vakuumgefäßes 
höchstens 25 cm lang sein. Umkehrstrecken (3) verlängerten auf das Doppelte, 
bzw. Dreifache. Die oberen lagen mindestens 4 cm unter dem Badspiegel. 


Kapillarenbäder waren: 
0,00° C (273,1° abs.): Eis in Wasser, von schwachem Luftstrom gerührt. 
—78,5° C (194,6° abs.)+ 0,1°: Festes CO, mit Benzin, mit CO,-Gas gerührt. 
—183,1°C (90,0° abs.) + < 0,1°: Flüssiger O,, mit technischem O,-Gas gerührt. 


Ein Pt- Widerstandsthermometer, bei 0,00°C und — 183,1°C, wiederholt 
geeicht, auf 0,05° ablesbar, maß T bei — 78,5° C. Absolutgenauigkeit: +0,15°. 

Ein selbstgefertigtes Rein-Sauerstoffdampfdruck-Thermometer maß sehr 
genau Temperaturschwankungen wie auch den Absolutbetrag bei —183,1°. 

Die Dampfdruckkurve reinen Sauerstoffs ist sehr genau bekannt (6), (7), (8). 
Der Dampfdruck betrug gewöhnlich etwa 750 mm Hg; er schwankte höchstens 
zwischen 730 und 775 mm Hg, was einem Temperaturunterschied von 0,57° ent- 
spricht. Reiner Sauerstoff hat den Dampfdruck von 760,0 mm Hg bei — 182,97° C, 
Alle Messungen wurden aber zweckmäßigerweise auf —183,1°C, d.i. 90,0° abs. 
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Darstellung, Beschaffenheit und Reinigung der Gase: 


Luft: Außenluft durchstrich Glaswolle, festes CaCl,, konz. H,SO, und 
festes KOH. 


H,: Durch Elektrolyse von etwa 30prozentiger Kalilauge gewonnen ging 
durch CaCl,, durch auf 300°C erhitzten Pd-Asbest, konz. H,SO, 
und KOH. 


O,: Wie bei H,; ebenfalls über Pd-Asbest u. s. w. 


He und Ne stellte uns die „Gesellschaft für Lindes Eismaschinen A.-G., 
Höllriegelskreuth‘“ zur Verfügung, mit der Angabe bei Helium mindestens 
99,5 Proz. He. Einer weiteren Reinigung wurde es nicht unterworfen. Das 
Neon bezeichnete die Firma mit: Rein. Das 7 lag bei 0° C um 0,68 Proz. 
unter dem Wert von M. Trautz und H.E.Binkele (9) und um 0,44 Proz. 
unter dem von Rankine; das Ne war also wohl mit etwas He und N, ver- 
unreinigt. Ging das Roh-Neon bei der Temperatur des flüss. O, über aktive, 
zuvor im Vakuum ausgeheizte Kohle in langsamem Strome [Edwards (2)], so 
lag n nur noch wenig unter dem Wert von M. Trautz und H. E. Binkele, 
etwa um 0,2 Proz., und kam auf den Wert von Rankine (10). 

Der Millikan-Luftwert für 23,0° C liegt mittelbar den Messungen bei 0,00°C 
zugrunde. Die Interpolationsformel von Millikan hat einen Gültigkeitsbereich 
innerhalb 0,1 Proz. bis herab zu 12°C: 


ne 107 = 1824,0 — 4,93 (23 — 2). 


Fiir 000° C ergibt die Formel den Wert 1710,6. Dieser Wert ist sicher etwas 
zu hoch, denn mit sinkender Temperatur fällt die Reibung von Luft etwas 
mehr als linear. Nimmt man jedoch auch für 0,0°C eine Genauigkeit von 
+ 0,1 Proz. an, so wäre der obere Grenzwert 1712,3 und der untere 1708,9. 
Der wirkliche Luftwert wird nach dem eben Gesagten mehr dem unteren Wert 
zuneigen. Einer Arbeit von Sutherland und Maaß (11) entnimmt man einen 
Luftwert bei 0,0°C von 1708,3 (Methode der schwingenden Scheibe). 

Daher wurde hier für 0,0°C der Luftwert: 1709,0 festgesetzt und zu- 
grundegelegt. Bei —78,5°C (194,6° abs.) ist auf den Reibungswert 659,0 für 
H, bezogen, bei —183,1° C (90,0° abs.) auf 370,6 ebenfalls für H,. Die beiden 
H,-Werte entstammen der Interpolationsformel (12): 


k T, = 33,18; Nu = 0,975 . 


Diese Formel gab aüch noch die Reibungswerte von H, bei 20° abs. richtig 
wieder. 

Schwankungen der Temperatur des Ausflußapparates wurden durch Zurück- 
rechnen mit Hilfe des idealen Gasgesetzes auf eine feste Temperatur aus- 
geschieden, ebenso Schwankungen der Kapillarentemperatur, indem man für 
wenige Zehntelgrad die Durchströmungszeiten den absoluten Temperaturen 
proportional setzte. 

Korrektion auf Zweierstöße und Unidealität bei Mischung wurde nicht 
angebracht. 

Korrektion infolge veränderten Barometerstandes während der Messungs- 
reihe war nicht erforderlich, da die Eichmessungen mit Luft, bzw. mit H, 
jeweils kurz vor der Haupt 
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2. Messungsergebnisse 
a) Reingasreibungszahlen 


und - — 
Reingas T° abs. n.10° Meßreihen | Bezogen auf 
g g 
1,80, 273,1 833,4 2 | Luft: 1709,0 
0, Ä 282,1 1954 | 1 | yy : 1754 
273,1 1905 2 » : 1709,0 
He FR 90,0 884,1 1 H,: 370,6 
6 aaa 273,1 2960,5 4 Luft: 1709,0 
Wr. (Roh-Neon) . (273,1 2949 | 1 „ : 1709,0) 
— 194,6 2356 1 H,: 659,0 
a 90,0 1311 2 „: 370,6 
’rO2. 
= T = 90,0° abs. bezogen auf ng, = 370,6 PR 
tive, 
)], 80 
cele, 
0°C 
reich 
twas | 
twas 
08.9. Zusammenfassung der Misch-7: 
ia Mischg, | 0 | 10 | 20 30 40 | 50 | 60 | 70 | 80 | 90 | 100 
- H,-He . |370,6414,4| 458,2) 504) 552,6 602,1) 653,6 708 | 764 | 823 | 884,1 
für H,-Ne . 370,6)527,8| 677,4; 808| 915,5! 955,7.1080,4 1151 |1214,711268,011311 
‘in He-Ne . | 884,1/981,3 1064,2|1129 1179,4 1218,5 1250 1276,6 1295,7,1307,1 1311 
4 3. Verwertung der Messungsergebnisse 
Jetzt wird genauer als bisher [Berechnungsformel für 77 aus 
7, I’, und nr vgl. (12)] (9) für 
chtig 
Helium: 7,-107 = 98,0 (bisher 97,1); n, = 0,664 (wie bisher); 
k ’ 9 
rück- Neon: 7,107 = 708,5 (bisher 705); n, = 0,620 (bisher 0,613). 
aus- 
de Edelgase: Ty, me, Mg; Tgemess.’ ber.» Abweichung beider voneinander, n,,, 
im Temperaturbereich von 90,0° bis 550° abs. max. 
| Helium (9) Neon (9) Argon (9) Krypton Xenon (5) 
.i T, 519 | 44,74 | 1506 | 210,6 | 289,7 
Ä 710? 80 | 701,5 1209 ~1780—90 | 2235 
führt n, 0,664 “0,620 | “qos | — | 0m 
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Edelgase (Fortsetzung) 
Helium (9) Neon (9) Argon (9) Krypton Xenon 6) 
| gemess. ngemess. = n gemess. ‘7 gemess. gemess.| = 
"| » ber. n n ber. nm, | 7 ber. | x» ber. in ber. | 
"loo | "loo 00 
Abwehg | Abwehg. Abwehg. |Abwehg. Abwehg. 
gsa | 1311 | 
90,0 883,4 0,855 1317 0,987; — - |j—| — 
—0,8 "lee +4,5 | 
1480 2356 
194,6 — (08274 2345 [0,850 — — |—| — Te 
| 2961,5 | | 2314 | 
273,1 — | 2954 0,934 
— 2,2" loo | 
1941 3092, 3102 2217 2260 
2930 1944 0,816 | 3095 10,930 2206 0,9887 — |— | 2260 (0,99 
+1,5°/o9 —5,0%o9 | £.0,0°/o9 | 
2281 | 3623, 3635 2691 | 3036 
373,0 2282 0,810 | 3626 (0,919 -2669 0,9806 -| - | = 
+0,4%/,n -8,2%/ | 
3009 
1000 — - | — 4- | - 
+6,69 
| 
| 3351 
450,0 _ _ —_ — |—| 3355 10,994 
|+1,2°/, 
2672 | 4220, 4235 3217 | [+ 
4730 2673 0,804 | 4227° 0.90 3185 09716 — — | - 
+0,4 9/90 +1,7—2,0%/0 — 109/90 3652 
-| +, 1-1 — |—| 3861 0906 
| | 
2858 | 4501, 4518 464 | 
5230 2853 0,802 | 4507 (0,904). $422 083 — — |- 
= 1,7 */o0 +1,3—2,4°/o, — | 3954 | 
-| — |—| — |—| 8951. 
Edelgase (und Wasserstoff): 
Temperaturunabhängige Kenngrößen 
H, | He | Ne | Ar = | we 
T, 33,18 5,19 44,74 | 150,65 210,6 | 289,7 
+10” | 158,0 98,0 701,5 1209 ~ 1780 2235 
ng 0,670 0,664 | 0,620 0,602 _ 0,520 
10 
a 8,231 1,297 2,217 3,776 | 2514 | 2,206 
el 
ac 2169 1246 | ~ 1183 | 1296 
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Edelgase (und Wasserstoff) (Fortsetzung) 
H, He Ne 33 


i Me 4,611 | 18,88 15,68 8,025 | ~ 8,45 | 7,715 

k 

ne 26,56 43,02 | 104,9 98,50 | ~122,7 | 131,3 
107,8 48,99 156,2 191,4 » 195,5 | 195,1 
| 4,440 | 4142 | 3,975 | 4,785 5,634 

Gyr. AE| 2,08 08 1,12 1,54 1,68 1,90 


Temperaturabhängige Kenngrößen der Misch- und Unpaargasreibung: 


| q | 


H,—He: Try: 13,12; tug? 125,53 0,667; aye 10°: 956,6 
13,12 125,5 0,77 5 4,203] 4,291 0,002 
90,0 | 596 | 602 0,961 0,760 oo 
298 1339 1350 81 0,974 | 2,4 0,98 0,902 2 E 
373 1577/1587] 0,972 | ,55 0,90 2006, 501 0,7482/2,806 


473 |1864/ 1861 0,965 |1,17 ‚6156.0, 7467/2,748) 2,763 


| 
| 
| 
| 
| 
| 
| 
| 
| 


523 |2001|1991 0,961 | — 0.99 0,904 0,858 6,6006 0,7461 2,728 2, 736 5 


H,—Ne: Tey: 38, 53; Ma: 4045 0,645; 10%: 1049 


38,53; — | 404| —| 1,283 | — | — [1,55 4,332 9208 
| 


90,0 850,2862,5 1,220|— 4, 11,51810,980015,167 0,770 3,768 3,681 | 


293 | 1889/2031 70 | | 1,81, |1,43 1,31, 1,503/0,9270 5,684 (0,7617| 3,043 3,126 
| 


0,087 


0,08 083 
10 ‚7599 2,964 3,067 | 


0, 7568 2, 928) 3,008 | 


373 | 2406/2388 ia 1,54, \1,09/2,0), 1,504.0,9166 5,801 0.103 


173 |2848|2790 

523 3108 |2988 49,1) 464, |1,82)2,12 1,511 0,90315,962 0,7560 2,906 2,985 
He—Ne: Tyg: 15,24; „ug: 311,5; 0,642; 10°; 489,2 

15,24; — 115 — | 1,188 | — | — 2,27 |1,0000/5,0000 0,775 13,982 4,059 


1,58, |3,3 |2,06, 1,509 0,9065 5,919 10.080 


0,079 


4059 0,077 
90,0 |1314 1313 1,22, |1,8| — 11,68410,9251|5,705 0,762 |3,32, 3,357 0.034 
| | I | ’ 
293 ime (8/2877 35] | 1,2 1,28 |1,620)0,8709 6,392 (0,7496 2,77, 2,805 930 
373 | 3380/3365 1,17, | 2,0 6,528 |0,747612,72, 2,746 0.021 
473 |3930 3924|,7 1,16, | 1,5 | 1,13 |1,617\0,8536|6,656 |0,7458)2,68, 2,688|,. < 
523 | 4183/4187|95,| 1,16, | 2,3 | 1,18 0256316, 709 0, 7449|2,65, | 2,665 
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Unpaargas-, Mischgas- und Reingas-Konstante. 
(Temperaturunabhängige Kenngrößen) vgl. Ann. d. Phys. [5] 18. S. 856. 1933, 


| 


Im, VT,, l/a, 


| 
| 
VT,,| la, 


—3 |13,12|125,5|0,667|4,293 3,622) 1045 0,667|1,025, 
— (33,18) 153,0/0,670|4,440 5.760 217,0 | | — 
— 5,19 98,0 0,664 4, 142 2,278 | 53,0 | — 
+18 38,53|404 |0,645|4,332 6,207 953,7 11,396 | 
H,-Ne| —  33,18| 153,0 0,670|4,440 5,760, 217,0 “> 
— 44,74 | 701,5 | 0,620 | 3,975 | 6,689 | 63,8 | | 
+10 15,24|311,5 0,642|3,982 3,904, 489,2 | 1,91, 12,09 
He-Ne| 5,19 98,0 0,664|4,142 2,278| 53,0 | | — | 
| — 6,689} 638 | — | | — 12018 
4. Zusammenfassung 
Tieftemperaturmessungen von 7 bei 194,6° und 90,0° abs. wurden 
an den Reingasen He und Ne, bei 90,0° abs. an den 3 binären 
Mischungen zwischen H,, He und Ne angestellt. 
Die schon berechneten Kenngrößen (Ann. d. Phys. [5] 16. S. 881. 
1933 und Ann. d. Phys. [5] 18. S. 816. 1933) konnten bestätigt, bzw. 
durch um nur wenig verschiedene ersetzt werden. 


Die Durchführung dieser Arbeit ist dem Entgegenkommen der 
Gesellschaft für Lindes Eismaschinen A.-G., Höllriegelskreuth zu 
danken, die He und Ne in ausreichender Menge zur Verfügung stellte 
und der Notgemeinschaft der Deutschen Wissenschaft. Weiter danken 
wir der 1.G.Farbenindustrie, Werk Ludwigshafen, für kostenlose 
Überlassung von fester Kohlensäure, flüssigem O, und N,, sowie dem 
Chemischen Institut der Universität Heidelberg für kostenlose Bereit- 
stellung von ausreichenden Mengen des flüssigen 1.-G.-Sauerstoffs. 
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Der Widerstand von Stäben 


Pa Mit 9 Figuren) 


3 Von @. Sprick und A. Becker 


Die gegenwärtige Untersuchung erstrebt in Fortsetzung früherer 
Messungen!) eine Vervollständigung der bisherigen verhältnismäßig 
beschränkten Kenntnis über den hydrodynamischen Gesamtwiderstand _ 
von Kreisstäben, indem sie zunächst den Widerstand quer zur Achse 
(Querwiderstand) unter Variation der Versuchsbedingungen neu er- | 
mittelt und insbesondere nach dem Gebiet kleiner Reynoldsscher _ 
Zahlen (R< 10) weiterverfolgt und dann zum Vergleich die ent- 


niires sprechenden Verhältnisse des Widerstands in Richtung der Achse 
(Längswiderstand) der Beobachtung unterwirft. ao 

. 881. Die Mefmethode ist die schon früher benutzte: Beobachtet en : 

‚ bzw. wird die in einer Flüssigkeit sich einstellende konstante Fal- _ 
geschwindigkeit, bei welcher der hydrodynamische Widerstand mit & 

oa dem Gewicht des in der Flüssigkeit untergetaucht ruhenden Fall- 

hb ‘ae körpers identisch ist. Dieses Gewicht kann in weiten Grenzen 

stellte durch Einfüllen verschiedener kleiner Quecksilbermengen in den 


Fallkörper verändert und mit Hilfe einer über dem Fallraum auf- 
gestellten hochempfindlichen Torsionswaage mit großer Genauigkeit 

» dan bestimmt werden. Die zugehörigen Fallgeschwindigkeiten ergeben 
evel sich aus den gemessenen Durchgangszeiten durch einige in be- 

fis. kannten Abständen voneinander angebrachte Visiermarken. Zur er- 
schütterungsfreien Auslösung der Fallbewegung dient eine elektro- 
magnetische Greiferanordnung, welche den Fallkörper in bestimmter 
Ausgangslage festhält. Als Stäbe kommen kreisrunde Glasstäbe mit _ 
sorgfältig halbkugelig abgerundeten Enden zur Verwendung, deren 
Dicke von 0,2—1,6 cm und deren Länge von etwa 5—20cm — 
variiert wird. Als Widerstandsmedien werden die beiden in ihrer 
Zähigkeit wesentlich voneinander verschiedenen Flüssigkeiten Wasser 


und Glykol benutzt. : 
1. Der Querwiderstand 


anken 
nlose 


Da die Stäbe bei der Ermittlung ihres Querwiderstands den ae 
ganzen Fallweg in streng horizontaler Lage durchsetzen müssen, 
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wird es erforderlich, sie mit einem Tragkörper zu verbinden, der 
an der Bewegung teilnimmt und dessen Einfluß auf den Widerstand 
durch Differenzbildung an den Beobachtungsergebnissen in Abzug 
> HR gebracht werden muß. Um die hier- 

durch verursachte Beeinträchtigung 

der Genauigkeit der Endergebnisse 
möglichst gering zu halten, haben 
wir als Tragkörper Glashohlkugeln 
von so geringen Abmessungen be- 
nutzt, daß ihr Einzelwiderstand 
niemals einen sehr erheblichen 
Bruchteil des Gesamtwiderstands 
erreichte. Wir haben außerdem 
jetzt besonderen Wert darauf ge- 
legt, durch systematische Ver- 
änderung der durch die Bifilar- 
wo; aufhängung der Stäbe am Trag- 
körper (vgl. Fig. 1a) gegebenen geometrischen Verhältnisse eine 
etwaige Störung des Strömungsbildes zu erkennen und auszuschalten. 


Fig. 1. Fallkörper 


a) Fallversuche in Wasser 


Um einerseits eine größere Variation in den Abmessungen der 
Fallkörper vornehmen und anderseits etwaigen Einwänden gegen 
eine zu starke Begrenztheit des Fallraums begegnen zu können, 
sollten die früheren Fallversuche in einem möglichst erheblich ver- 
größerten Fallraum wiederholt werden. Dies wurde uns er- 
möglicht durch das freundliche Entgegenkommen der Direktion 
der I.G.Farbenindustrie in Ludwigshafen und insbesondere des 
Herrn Direktor Dr. O. Schmidt, durch welches uns ein großer Be- 
hälter zur Verfügung gestellt und das Arbeiten in dem Werke selbst 
gestattet worden ist. Der zylindrische Behälter hatte einen Durch- 
messer von 60 cm und eine Höhe von 140 cm und war mit 
6 großen Beobachtungsfenstern versehen, so daß bei Verwendung 
einer Anlaufstrecke bis zu 30 cm zwei aufeinanderfolgende Fall- 
strecken von je etwa 40 cm für die Messungen verfügbar waren. 
Das der Leitung entnommene Wasser wurde durch groben Sand 
und Holzkohle filtriert und von Luftblasen befreit. Zur Durch- 
mischung diente ein großes Drahtsieb, das auf und ab bewegt 
werden konnte und gleichzeitig das Heben der gesunkenen Fall- 
körper ermöglichte. Zur Verminderung der Temperaturveränder- 
lichkeit war der Behälter mit mehreren Lagen von Hanfseilen um- 
wickelt, und aus der  Ablesung in Höhe 
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befestigter Thermometer konnte der jeweilige Durchschnittswert der 
Versuchstemperatur hergeleitet werden. = 

Da die Beruhigung der verhältnismäßig großen Wassermenge _ 
immer eine erhebliche Zeit beanspruchte, so war die Durchführung 
längerer Meßreihen sehr zeitraubend. Wir haben trotzdem stets ; 
mehrere Geschwindigkeitsmessungen beim gleichen Gewicht des 
Fallkörpers gemacht, so daß alle Angaben des Folgenden Mittel- — 
werte aus mehreren, meist nur um wenige Hundertteile von- 
einander verschisdenen Einzelbeobachtungen darstellen. Wir haben 
auch jede Wägung sowohl vor wie nach dem einzelnen Fall- En 


uns aus dem Vergleich der zu den beiden Fallstrecken gehörigen x 
Durchschnittsgeschwindigkeiten zu versichern, daß wir tatsächlich 6% 
konstante Endgeschwindigkeiten messen. Es zeigte sich, daß eine Het 
Anlaufstrecke über 25 cm in allen Fällen ausreichte, um die End- I 
geschwindigkeit mit einer Genauigkeit von 1°/, zu erzielen, während Fi 
bei geringerer Anlaufstrecke größere Fehler merkbar waren. Um me 
festzustellen, wie weit eine Fälschung der Widerstandswerte durch 
die Verbindung der Stäbe mit der Tragkugel möglich sein könnte, a 
haben wir die Länge der Bifilaraufhängung, die entweder aus — 
Platindraht von 0,03 mm Dicke oder aus feinstem Kokonfaden 5 
bestand, von 5—20 cm verändert und in keinem Fall eine merk- = E 
liche Veränderung der Fallgeschwindigkeit beobachtet; es wurden | 
deshalb meist Abstände h (Fig. 1a) von 8—10 cm gewählt. 
Es war nicht zu erwarten, daß mit dieser Anordnung kleinste 
Widerstände von der Größe weniger Dyn mit der relativen Ge- 
nauigkeit der früheren sorgfältigen Laboratoriumsversuche meBbar : 
werden könnten. Wir haben daher meist dickere und zum Teil © 
auch etwas längere Stäbe benutzt, ohne außer Betracht zu lassen, 
daß insbesondere solche Abmessungen der Fallkörper verwendet 


werden sollten, die möglichst klein sind gegenüber den Abmessungen et 
des Wasserbehilters. Wir haben damit von vornherein darauf ver- 


zichtet, mit diesen Beobachtungen zu sehr kleinen Reynoldsschen © 
Zahlen vorzudringen; wir erwarten von ihnen lediglich eine — 
der vorhandenen 
Einen Auszug aus den MeBergebnissen gibt Fig. 2. Alle Beob- > DR 
achtungswerte sind in bekannter Weise auf die Temperatur von u: = 
korrigiert. Aus den ausgleichenden Kurven ermitteln wir durch er 
Differenzbildung den zur jeweils ausgewählten Geschwindigkeit v ae 
gehörigen Widerstand W der Längendifferenz | zweier Stäbe gleicher ee ie 
Dicke d= 2r und ‚ daraus die durch die übliche Be- 


ziehung W=y-od- 
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Funktion der Reynoldsschen Zahl R = 


Band 22. 1935 


d ‘ 
v-—, wo v die kinema- 


tische Zähigkeit und ö das spezifische Gewicht des Mediums ist. 


. Staboiche 0 82cm 


Tragkugel r=0:78 


Widerstand Win Dyn 


T 


l 


beschwindhgkeit v 


| 


3 
Widerstandsmessung in Wasser 


digend iibereinstimmt mit der in der 


4 Sek 


Wie aus Fig. 3 
ersichtlich ist, ist 
zwar eine gewisse 
Streuung der Ein- 

at zelwerte vorhanden, 
die aber auf keine 
systematischen Ab- 
weichungen der aus 
den verschiedenen 
Kombinationen fol- 
genden Daten von 
einem durchschnitt- 
lichen Verlauf hin- 
weist. Man ersieht 
ferner, daB dieser 
durchschnittliche 
Verlauf sehr befrie- 


friiheren Untersuchung an 


Stäben anderer Abmessungen gewonnenen Kurve. Wir sehen darin 


Kurve Vogt u. Becker 
Aus la » der fig.2 


eine Bestätigung der 
Allgemeingültigkeit 
des festgestellten 
w-R-Gangs für den 
untersuchten R-Be- 
reich. 


b) Fallversuche 
in Glykol 
Um den Gang 
der Widerstandszahl 
nach kleinen R- 
Werten hin zu ver- 


BE folgen, sind, wie 


früher schon hervor- 


gehoben, Beobachtungen der Fallbewegung in Wasser nicht geeignet, 
da es sich hier um Widerstandskräfte handeln würde, die nicht mehr 
mit genügender Sicherheit meßbar wären. Es war daher die Wahl eines 
Mediums mit wesentlich größerer kinematischer Zähigkeit erforderlich. 
Wir haben als solches Glykol benutzt, .das wir ebenfalls dem Ent- 
gegenkommen der Direktion der I. G. Farbenindustrie in Ludwigs- 
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hafen verdanken. Um keine allzu große Flüssigkeitsmenge zu be- 
nötigen, haben wir als Behälter den zu den älteren Versuchen 
verwandten, mit ebenen Wänden versehenen großen Glastrog mit 
quadratischem Querschnitt von 32,5x 32,5 cm? und 82cm Höhe 
gewählt, dessen Größe nach dem Ausfall der Wasserversuche als 
ausreichend erscheinen durfte. Zu weitgehender Prüfung der Er- 
gebnisse haben wir unter starker Variation der Versuchsbedin- 
gungen eine große Zahl sorgfältiger Messungen ausgeführt. Dies 
war deshalb mit 
nicht allzu großem 
Zeitaufwand mög- 
lich, weil das zähere 
Medium stets in 
kurzer Zeit zur Ruhe 
kam, wie deutlich 
aus dem Verhalten 
vorhandener suspen- 9 
70 


200 


ad 


dierter Teilchen zu 
erkennen war. Aus 
diesem Verhalten er- 
gab sich auch ein 
einwandfreies Bild 
von den Strémungs- 79.2990 
verhältnissen in der Traghugel r für 0-02 _| 
Umgebung der Fall- + 
körper, das mit 
Sicherheit auf prak- Geschwindigkeit v 
tische Störungsfrei- N l N 
in allen quantitativ = 
ausgewerteten Fällen 

schließen ließ. Da 

Veränderungen des Glykols während der Untersuchung nicht 
ausgeschlossen waren, so wurde sowohl das spezifische Gewicht 0 
wie die innere Reibung 7 für jede Versuchsreihe neu bestimmt. 
Für 5 fanden sich Werte, die, auf die Temperatur von 21° be- 
zogen, nur zwischen 1,121 und 1,123 schwankten, so daß der 
Wert 1,122 gem”? als Konstante festgehalten werden kann. Der 
Wert von 7 zeigte demgegenüber im Laufe der Versuche eine ständige 
Abnahme von 0,1855—0,1772 g em”!sec”! (bei 21°), was im 
einzelnen zu berücksichtigen war. Daneben war mit besonderer Sorg- 
falt auf die Temperaturveränderlichkeit von 4 zu achten. Zur Um- 
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rechnung des Wider- 
stands auf eine feste 
Temperatur — wir 
haben hierfür 21° 
gewählt — stand 
eine schon früher 
benutzte einfache 
Beziehung!) zur Ver- 
fügung. 

In Figg. 4 und5 
sind die MeBergeb- 
nisse fiir 4 zwischen 
02 und 08 cm 
liegende Stabdicken 
und mehrere Stab- 
längen graphisch 
wiedergegeben. Sie 
zeigen den Verlauf 
des Gesamtwider- 
stands von Stab und 
Tragkugel in Ab- 
hängigkeit von der 
Geschwindigkeit, je- 
weils gemittelt aus 
mehreren Einzel- 
messungen und auf 
die genannte feste 
Temperatur korri- 
giert. Man _ sieht, 
daß alle Beobach- 
tungen durch sehr 
befriedigend gesi- 
cherte Kurven zu- 
sammengefaßt wer- 
den können. Die- 
selben stützen sich 
auf ein so reiches 
Erfahrungsmaterial, 
daß von ihrer Aus- 


1) K. Vogt und 


A. Becker, Ztschr. f. 
Phys. 83. S. 168. 1933. 
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wertung neue einwandfreie Kenntnis des reinen Stabwiderstands 
erwartet werden muß. 

Ein Gesamtbild von dieser neuen Kenntnis gibt der zunächst 
für kleinere Reynoldssche Zahlen in Fig. 6 dargestellte w-R-Zu- 
sammenhang, wie er aus der geeigneten Kombination der zu ver- 
schiedenen Stablängen gehörigen Kurvenwerte der Figg. 4 und 5 
erhalten wird. Dieser Zusammenhang konnte, wie man sieht, mit 
befriedigender Genauigkeit bis zu etwa R= 0,5 herab festgelegt 
werden, und ein systematischer Gang der w-Werte mit den zu 
ihrer Herleitung benutzten Stabdimensionen macht sich nicht be- 


Aurve nach Vogt u. Becker 
I x  Fallbecbachtungen in 6lykol 
2 N 
8 » + /allbeobachlungen in Glykol Langs-Widerstand 
N Stabdicke d-0.625 cm 
. 7.0 
\ + ” 15 
47 
s+ 
- 
Oe 
+ 
0 DD DH 0 DU DD DW 
Fig. 7 


merkbar. Hervorzuheben ist, daß das angewandte Subtraktions- 
verfahren auch den Tragkugeleinfluß auf den gesuchten Stab- 
widerstand praktisch völlig ausschaltet, wie aus der Unabhängigkeit 
der w-Werte von der Größe der beteiligten Tragkugeln hervorgeht. 
Der starke Anstieg der Kurve zu relativ hohen w-Werten ent- 
spricht ganz der durch die auf kleinste R-Werte beschränkte hydro- 
dynamische Berechnung Lambs?) nahegelegten Erwartung. 
Vereinigen wir diese aus den Messungen in Glykol gewonnenen 
Ergebnisse mit den durch Fig. 3 veranschaulichten Ergebnissen der 
Wassermessungen, so erhält man das in Fig. 7 durch Kurven 
zusammengefaßte Bild, welches das Gebiet der Reynoldsschen 
Zahlen von etwa 2—130 umfaßt und einen anfänglich steilen, 


1) H. Lamb, Phil. Mag. 21. 8.10.11. 0000000000002 


a 

ider- 
feste 

. 
wir 
21° 

tand 
üher 
ache . 
Ver- 
- 
ndd 
tab- 
: 
isch 
sie 
der 
Ab. 
je- 2.7 
zel 
este 4 
eht, 
‚ch- 
h 
ehr © 
esi- 
)ie- 
ich 
hes 
al, | 
us- . 
ınd 

4 

XUN 


202 Annalen der 5. Folge. a: 


dann langsamen Abfall der w-Werte von etwa 8,5 auf 1 anzeigt. 
Man sieht, daß dieser Gang in dem schon früher untersuchten 
Bereich von R = 10—120 mit dem damals gefundenen (gestrichelte 
Kurve) völlig übereinstimmt, wie wir schon im vorhergehenden Ab- 
schnitt a) bemerkt haben, und daß auch dort, wo sowohl Wasser- 
wie Glykolmessungen für den festgestellten w-Wert in Betracht 
kommen, eine gute Übereinstimmung vorhanden ist. Die früher 
schon erkannte Abweichung der von Wieselsberger!) aus Wider- 
standsmessungen an Stäben im Luftstrom hergeleiteten w-R-Kurve 
(strichpunktiert) im Sinne eines schwächeren Abfalls nach größeren 
R-Werten hin bleibt bestehen. Die zusammengehörigen Zahlen- 
werte unserer Messungen enthält beistehende Tabelle. 


Tabelle 1 


80 100 
1,3 1,2 


Zur weiteren Veranschaulichung dieser Kenntnis kann die Be- 
trachtung des auf die Längeneinheit bezogenen Widerstands von 
Stäben in seiner Abhängigkeit von der Dicke und Geschwindigkeit 
dienen. Wir verweisen hierfür auf die Kurven der späteren Fig. 9, 
welche den Querwiderstand für die vier verwandten Stabdicken und 
Glykol mit der inneren Reibung 7 = 0,185 angeben. Beim Stab 
0,64cm sind die mit zwei verschiedenen Tragkugeln gefundenen 
Meßpunkte getrennt eingetragen, um ihre befriedigende Überein- 
stimmung zu zeigen. Bei den andren (gestrichelten) Kurven haben 
wir auf die Wiedergabe der Beobachtungspunkte verzichtet und 
ihre Berechnung mit Hilfe der ausgleichenden Kurve der Fig. 6 
bzw. der Angaben der Tab. 1 unter Benutzung der einfachen Be- 
ziehung W = .v-R- wy (W= Widerstand der Längeneinheit) vor- 
gezogen. Der Gang des Widerstands mit der Geschwindigkeit ist 
bei allen vier Kurven nahe der gleiche. Beachtenswert und prak- 
tisch von Bedeutung ist die Tatsache, daß der Widerstand mit zu- 
nehmender Stabdicke nur langsam wächst, so daß anderseits bei 
sehr kleinen Dicken noch recht nennenswerte Widerstände zu er- 
warten sind. Dies trifft ausgeprägt im Bereich kleiner Reynolds- 
scher Zahlen zu, wo die Abnahme von R mit einer erheblichen 


1) C. Wieselsberger, Phys. Ztschr. 22. S. 321. 1921; 23. S.219. 1922. 
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Zunahme von w verbunden ist, so daß das für W maßgebende Pro- 
dukt R- nur langsame Änderungen erfährt. 
Es tritt hier die Frage auf, ob nach dem von uns festgestellten 

— R-Gang der Fall eintreten könnte, daß abnehmendes r nicht 
immer mit abnehmendem W verbunden wäre. Die Betrachtung des 
4aW E AR 

R 
änderung AW/W im untersuchten R-Bereich stets das gleiche Vor- 
zeichen mit Ar/r hat, da die zu einer positiven relativen Änderung 
der Reynoldsschen Zahl gehörige negative relative Änderung von 
vw immer von kleinerem Absolutbetrag ist als die erstere. TR 

Es ist noch von Interesse, den Einfluß der Stabenden auf den Br 

Widerstand zu untersuchen. Wenn man auch eine noch weiter ge- 
steigerte Genauigkeit der Messungen wünschen möchte, um aus 
ihnen durch Differenzbildung zur genügend sicheren Beantwortung — 
dieser Frage zu gelangen, so kann man doch in erster Annäherung 
versuchen, durch Subtraktion des getrennt gemessenen Tragkugel- _ 
widerstands von dem zur jeweils gleichen Geschwindigkeit gehörigen 
beobachteten Gesamtwiderstand den reinen Widerstand desganzen Stabs 
herzuleiten und diesen mit demjenigen Wert zu vergleichen, den _ 
man durch a der Gesamtlänge mit dem auch zur vor 


Ausdrucks 


+ a zeigt, daB die relative Widerstands- 


in Tab. 2 einen solchen Vergleich { A die Stabdicke d = = 0,64 cm, 2 
4 


Tabelle 2 
d = 0,64 em 


eee 2 3 4 


Tragkugel r= 06 cm 
2,15 4,91 9,00, 13,5 
21,6 494 | 90,5 | 136 
20,2 46,2 84,8 | 127 
20,8 47,5 85 | 128 
Tragkugel r = 0,78 cm 
W. — W. ‘ 
= 488 a 2,15 4,82 8,92 13,3 | 
10,06 - W 21,6 48,5 90 134 
19,7 | 47,5 | 82 | 124 


W 10,06 ~ W eugel 


W = Weugel 
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die in drei Langen mit zwei verschiedenen Tragkugeln untersucht 
worden ist. 

Man erkennt, und dies findet sich auch fiir die anderen be- 
nutzten Stabdicken und Längen, daß die auf die Längeneinheit be- 
zogenen Widerstände W sowohl von den zu ihrer Ermittlung heran- 
gezogenen Absolutlängen wie von der gewählten Kugelgröße prak- 
tisch unabhängig sind. Man sieht ferner, daß gegen das Stabende 
hin meist eine Erniedrigung des Widerstands von etwa solchem 
Ausmaß angezeigt ist, daß der Halbkugelwiderstand der End- 
begrenzung vollständig zurücktritt. Ob diese Erscheinung quantitativ 
von der Geschwindigkeit abhängt, wie es bei den kleinsten Ge- 
schwindigkeiten angedeutet wird, dürfte mit Sicherheit wohl erst bei 
noch gesteigerter Genauigkeit der erforderlichen Beobachtungsdaten 
zu entscheiden sein. 


2. Der Längswiderstand 


Zur Untersuchung des Widerstands bei ihrer Bewegung in 
Achsenrichtung haben wir hohle Glasstäbe mit streng halbkugelig 
verschlossenen Enden benutzt. Dieselben besitzen, wie Fig. 1b zeigt, 
oben eine etwa 2 mm weite Bohrung, die gut abschließend mit einem 
kleinen Gummistöpsel mit feiner Drahtöse verschlossen ist, die zum 
Anhängen der Stäbe an die schon beschriebene Greifervorrichtung 
dient. Zur Erzielung verschiedener Fallgeschwindigkeiten wird wieder 
eine variable Quecksilbermenge Hg in den Fallkérper eingefiillt. 
Es ist von großem Vorteil, daß beim vertikalen Fall der Stäbe be- 
sondere Tragkörper und Aufhängevorrichtungen entbehrlich sind und 
daß infolgedessen unmittelbar die reinen Stabwiderstände gemessen 
werden können. Allerdings wird es auch hier erforderlich, durch 
Differenzbildung den Einfluß der Stabenden und der Mantelflächen 
getrennt herzuleiten. Es muß auch erwähnt werden, daß bei diesen 
Fallbewegungen, namentlich bei den dickeren Stäben und kleineren 
Geschwindigkeiten, bisweilen ein geringes Pendeln des oberen Stab- 
endes bemerkbar war, dessen Auftreten wahrscheinlich die in 
manchen Fällen sich zeigende geringere Sicherheit der Beob- 
achtungsdaten erklärt. 

Nach orientierenden Versuchen in Wasser, welche die relative 
Kleinheit des Längswiderstandes erkennen ließen, haben wir wieder 
für die ausführlichen Messungen Glykol als Widerstandsmedium 
benutzt. Zur Begrenzung des Fallraums diente jetzt ein 146 cm 
hoher Glaszylinder!) von 22 cm Durchmesser, in welchem wir eine 


1) Der uns ebenfalls von Herrn Direktor Dr. O. Schmidt freundlichst 
überlassen wurde. 
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Sprick u. Becker. Der hydrodynamische Widerstand von Stäben 205 


als ausreichend erkannte Anlaufstrecke von 40—50 cm und eine 
Meßstrecke von 60—70 cm verwenden konnten. Die Untersuchung 
erstreckt sich aus die drei Dicken 0,62, 1,0 und 1,59 cm mit je- 
weils vier Längen zwischen rund 6 und 20 cm. 

Einen Auszug aus den auf „= 0,185 korrigierten Versuchs- 
ergebnissen zeigt Fig. 8. Es liegt hier, wie man sieht, ein reiches 
Beobachtungsmaterial vor, das sich den ausgleichenden Kurven meist 
sehr befriedigend an- 
schließt. Von den + 
möglichen Folge- 
rungen aus diesen 
Kurven seien in den 


Tabellen 3 und 4 P™ 
die auf die Stab- “| & 
dicken d, = 0,62 und oe 
d, = 1,59 cm bezüg- 
lichen Verhältnisse 


im einzelnen an- 
gegeben. Dieselben 
verzeichnen für je- 
weils sechs Ge- 
schwindigkeiten die 
Werte des auf die # 
WO 
a 


Längeneinheit bezo- 
genen reinen Wider- 
stands W der Mantel- 
fläche, wie er sich 
durch einfache Sub- 
traktion aus der 
Kombination des 


Tabelle 3 + 
d, = 0,62cm 1, = 19,92 15,04 1,;,= 10,07 1, = 5,60 cm 


V em/see | 1 | 6 

(W,- Wp: 488 | 061 | 1,23 | 225 | 307 | 390 | 4,71 

(W, — W,): 9,85 0,66 | 1,32 | 223 | 304 3,96 | 4,67 

_ (W,- W):14,32 | 0,63 | 1,33 | 223 | 307 | 3,98 | 4,75 

W Mittel | 0,63 1,29 | 224 | 3,06 | 3,95 
ee W,-19,3 - W 2 | 6 8 13 19 
a: W,-144 -W 2 6 8 13 19 
Wa 9,45. 7 1,5 6 8 13 19 
W,- 498. W 2 6 18 
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Tabelle 4 
d,=159cm 1,= 20,55 ZL=1536 4=1056 1, = 6,19 cm 


V em/sec 


(W,- W,): 5,19 
(W,— W,): 10,0 
(W,—- W,): 14,36 
v Mittel | 1,36 | 2,59 
W, -18,96 - W 6 17 
W, —13,77 - W 6 16 
W,- 9,00. W 18 
W,- 4,60: W 7 18 
Ww Kugel d, 5 12 : | 42 


Widerstands von je zwei verschiedenen Stablängen ergibt. Man 
kann daraus 

Mantelfläche berechnen und 
nach einer Subtraktion vom 
Gesamtwiderstand auf den 
Beitrag der halbkugeligen End- 
flächen schließen. 

Man erhält für alle zur 
Berechnung benutzten Längen 
auffallend gut übereinstim- 
mende Werte von IV, woraus 
zu entnehmen ist, daß der 
0m Widerstand der Mantelfläche 

zum mindesten bis zu der 

kleinsten benutzten Stablänge 
9,:0625:m herab der Länge streng pro- 
portional ist. Wir haben 
seinen hier nachgewiesenen 
Gang mit der Geschwindigkeit 
in Fig. 9 graphisch einge- 
tragen. Vergleicht man diesen 
Gang für die gleiche Stabdicke 
beim Quer- und Längswider- 
stand, so zeigt sich ein ähnlicher, wenn auch nicht zahlenmäßig 
vollständig übereinstimmender Verlauf, wie dies die Fig. 9 für die 
Dicke 0,6 veranschaulicht. Reduziert man die den Querwiderstand 
angebenden Ordinatenwerte im Verhältnis 4,2:1, so kommt man zu 
den mit + bezeichneten Punkten, deren Abweichung von der ent- 
Kurve ‘Langswiderstands einen etwas schwi ächeren 


@verwicerstand 


W Stabwiderstand Dyn/cm 


T 


Längswiderstand 


Geschwindighert v 


‘ 
aver 
206 
— 
= 
70 | 8,60 
5 
55 | 8,43 
Le 
| 
; 
pe 
| 
~ 
er 
= 
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des Gangs von W mit der Stabdicke hinter der Proportionalität — 
haben wir auch beim Querwiderstand schon festgestellt. 

Ein Vergleich beider Widerstandserscheinungen ist auch auf 
dem Wege möglich, daß man versucht, den Längswiderstand durch 
eine in gleicher Weise wie beim Querwiderstand definierte Wider- 
standszahl als Funktion von R darzustellen. Man erhält so den 
in Fig. 7 eingezeichneten, für alle drei Stabdicken nahe gleichen 
Gang, der, von den erheblich tiefer liegenden Absolutwerten ab- 
gesehen, dem für den Querwiderstand gefundenen ähnlich ist. 

Besonders auffällig ist die große Übereinstimmung in den aus 
den verschiedenen Stablängen hergeleiteten Widerstandsbeiträgen 
der Stabenden. Sie zeigt für unsere Fälle eine gewisse Berechti- 
gung des groben einfachen Verfahrens der Addition der Teilwider- 
stände an, da sie keinen Einfluß der Länge der Verbindungsstücke 
auf die für den „Endwiderstand“ maßgebenden Strömungsformen 
an den Endbegrenzungen erkennen läßt. Daß indes ein solcher 
Einfluß bei kleineren Längen unterhalb 5 cm tatsächlich, wie zu 
erwarten, vorhanden ist, wird daraus ersichtlich, daß der festgestellte 
Endwiderstand nicht mit demjenigen einer aus den Endflächen ge- 
bildeten Kugel übereinstimmt, den wir am Schluß der beiden vor- 
stehenden Tabellen zum Vergleich angegeben haben. 

Man sieht, daß es möglich erscheint, auf dem hier eingeschla- 
genen Wege zu klaren Einblicken in die Widerstandsverhältnisse 
zu kommen. Indes wird man den Versuch, in der Frage des Längs- 
widerstands zu weiteren Einzelheiten quantitativ vorzudringen, 
zurückstellen müssen, bis ein noch reicheres, kritisch gewonnenes 
Beobachtungsmaterial vorliegt. Es wird eine besonders wichtige und 
vielleicht nicht ganz leichte Aufgabe der weiteren Versuche sein, 
die Sicherheit zu erlangen, daß die feststellbaren Widerstände völlig 
frei sind von einem Einfluß der oben bereits a 
Pendelerscheinungen. 


Zusammenfassung 


Zur Bestimmung ihres hydrodynamischen Widerstands wird 
die Fallgeschwindigkeit horizontaler und vertikaler Kreisstäbe mit 
halbkugeligen Enden in Wasser und Glykol beobachtet. 

Für den „Querwiderstand“ der zylindrischen Mantelflächen wird 
der durch frühere Untersuchungen gefundene Zusammenhang zwischen 
Widerstandszahl yw und Reyroldsscher Zahl R in großdimensio- 
niertem Raum geprüft und bestätigt und nach der Richtung kleiner 
EIERN Zahlen bis zu etwa 0,8 nat weiter verfolgt. Der 


4 Ps 
: Sprick u. Becker. Der hydrodynamische Widerstand von Stäben 207 oe 
q Anstieg der letzteren erkennen läßt. Das erhebliche Zuriickbleiben ae 
6 
3,29 4 
3,60 
3,43 
44 : 2 ‘ 
} 
2 
Man 
chen | = 
4 
vom 
den 
‚nd- } 
zur 
gen 3 a 
im- Bee 
che 
der 
nge 
ro- Be 
en 3 
1en 
eit 
sen 
er- 
tig 
lie 
nd 
zu 
en 2; 
2 \ 
XUN 
i 


208 Annalen der Physik. 5. Folge. Band 22. 1935 


festgestellte erhebliche Anstieg des w-Wertes in diesem Gebiet en 
spricht den auf kleinste R bezüglichen Aussagen der Lambsche 
Berechnung. Aus dem Gesamtwiderstand der begrenzten Stablänget 
wird außerdem der Widerstand der Endflächen hergeleitet und ala 
relativ geringfügig gefunden. 

Für den „Längswiderstand“ wird ebenfalls der Einfluß sowol 
der Mantelfliche wie der Endflächen in seiner Abhängigkeit va 
der Geschwindigkeit und den geometrischen Dimensionen unter: 
sucht. Die Ergebnisse sind von bemerkenswerter Einfachheit ur 
der praktischen Auswertung unmittelbar zugänglich. 


Der Direktion der I. G. Farbenindustrie A.-G. in Ludwigshafen 
a. Rh. und insbesondere Herrn Direktor Dr. h. c. Dr. O. Schmid 
schulden wir großen Dank dafür, daß sie uns die zur Durchfüh: 
rung der Arbeit erforderlichen Hilfsmittel sehr entgegenkommen 
zur Verfügung gestellt haben. 


Heidelberg, Physikal.-Radiol. Institut der Universität. 


(Eingegangen 11. Dezember 1934) 


Verantwortlich: für die Redaktion: Prof. Dr. E. Griineisen, Marburg/L.; für Anzeigen? 
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Zeitschrift 
fürtechnische Physik 


Im Auftrage der deutschen Gesellschaft für technische os 
Physik E. V. herausgegeben von Prof. Dr. Carl Ram- SF 
sauer, Berlin, und Prof. Dr. Hans Rukop, Berlin. | 
Redaktion: Dipl.-Ing. Prof. Dr. Wilhelm Hort, ee 
1935 erscheint der 16. Jahrgang. Jährlih 12 Hefe. 
Halbjährlich RM. 24.— 
Neben Originalarbeiten werden zusammenfassende Berichte ge- 
boten aus allen Gebieten der technischen Physik unter Berück- 
sichtigung der technischen Erfordernisse, der neuesten Errungen- 


schaften und Veröffentlichungen physikalischer und 
technischer Natur. 


Heft 11 und 12 enthalten den Bericht über den 


10. Deutschen 
Physiker- und Mathematikertag ag 
1934 Bad Pyrmont 


Es ist zu beachten, daß dieser Bericht nicht als Sonderausgabe erscheint, 
sondern vollständig nur im Rahmen der „Zeitschrift für tech- 
nische Physik“ veröffentlicht wird. Er umfahjt insgesamt 71 Vo- 
träge über folgende Stoffgebiete: : 


er 


Mechanik elastischer Körper - Mechanische Schwingungen - Gas- 
entladungen - Elektrizitätsdurchgang durch Festkörper und 
Flüssigkeiten - Magnetik - Elektrische Meßtechnik - Allgemeine 
Optik - Spektroskopie - Metallphysik - Stoffphysik - Photo- 
graphische Prozesse, Photochemie - Geophysik - Glasphysik - 


Radiologie - Allgemeine Meßtechnik - Allgemeine technische 


sik +» Physikalische Chemie - Aerodynamik - Strömungsphysik - 
Allgemeine Wärmelehre - Technische Wärmelehre - Wärmeüber- 
tragung, Wärme- und Kälteschutz - Kältetechnik - Allgemeine 
Elektrik + Elektrizitätsaustritt aus festen Körpern - Nach- 
richtenübermittlung durch Draht und Raum - BE 
Bildtelegraphie - Elektrische Schwingungen (Sende- und Emp- 
fangsapparate) - Lichttechnik - Allgemeine technische Akustik- 
Hörtechnik - Allgemeine Chemie. 


Umfang beider Hefte: 96 und 190 Seiten mit zahlreichen At a 
bildungen im Text und auf Tafeln, sowie Tabellen. RM, 16.— 
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Grundlagen der technischen 
Röntgendurchstrahlung 


Von Dr.-Ing. Rudolf Berthold, Berlin. VI, 109 Seiten mit 62 Fi- 
guren im Text. 1930. gr. 8°. RM. 8.10, geb. RM. 9.72 


Inhalt: Allgemeine Grundlagen: Nutzstrahlung. — Schwächung 
der Nutzstrahlung beim Durchgang durch Materie. — Photo- 
graphische Schicht ohne Verstärker. — Photographische Schicht 
mit Verstärker. — Leuchtschirm. — Physiologische Grundlagen. 
Die besonderen Grundlagen des Verfahrens: Belichtung. — Fehler- 
erkennbarkeit. — Grenzen des Verfahrens, Technische Hilfsmittel: 
Hochspannungserzeuger. — Röntgenröhre. — Sonstige Hilfsmittel, 
Röntgenanlagen: Schutzmaßnahmen, Wirtschaftlichkeit des Ver- 
fahrens: Kostenberechnung von Röntgenaufnahmen. 


Physikalische Zeitschrift: Das klar und ansprechend geschriebene Biichlein wendet sich in 
erster Linie an den Techniker. Der Techniker, aber auch der Physiker, der im Laboratorium 
mit Röntgenstrahlen zu arbeiten hat, wird viel Anregendes und Belehrendes in dem Büch- 
lein finden. L. Bewilogua, 


Einführung in die Physik der 


Réntgen- und Gammastrahlen 


Von Maurice de Broglie, Membre de |' Institut, Paris, Louis 
de Broglie, Docteur des Sciences, Paris, Nobelpreisträger für 
Physik 1929. 


Von den Verfassern ergänzte deutsche Ausgabe. Übersetzung 
von L.vonSeuffert, München. Mit Einführung von L. Graetz, 
München. 


VIil, 208 Seiten mit 33 Abbildungen im Text und 11 Tafeln. 
1930. gr. 8°. RM. 18.90, geb. RM. 20.70 


RA 


VDI-Zeitschrift: Das Buch behandelt in erschöpfender Weise alle Erscheinungen auf dem 
Gebiete der Röntgen- und Gammastrahlen. „Insbesondere bietet es auch sehr wertvolle Aus- 
blicke auf den Widerstreit der elektromagnetischen Theorie und der Quantentheorie der Strah- 
lung. In den letzten Abschnitten leiten die Verfasser in die wichtigsten Fragen der Wellen- 
mechanik über, die ja von L. de Broglie der Wissenschaft zuerst erschlossen wurde. Wegen 
der gut gelungenen Zusammenfassung aller Beobachtungen und versuchten Erklärungen auf 
diesem Gebiet, bei der auch immer auf die noch nicht gelösten Aufgaben hingewiesen wird, 
ist das Erscheinen einer deutschen Ausgabe des wertvollen Buches sicherlich besonders zu 
begrüßen. 


Verlangen Sie bitte mein ausführliches Fachverzeichnis ,,Naturwissenschaften” 
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